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Abstract 

From  October,  1982  to  October,  1983  a  current  meter  mooring  reaching 
from  the  bottom  into  the  thermocllne  was  deployed  for  the  first  time  In 
the  Gulf  Stream  at  68*W.  The  temperatures,  pressures,  and  velocities  at 
the  uppermost  instrument  Indicate  the  Gulf  Stream  moved  back  and  forth 
across  the  mooring  site,  so  that  the  entire  Stream  was  sampled  In  time; 
hence  the  data  may  be  used  to  examine  horizontal  as  well  as  vertical 
structure  of  the  Stream.  The  two  key  points  to  the  success  of  the 
analysis  are:  1 ) the  well-defined  relationship  between  temperature  and 
cross-stream  distance  In  the  thermocllne,  enabling  the  use  of  the  former 
as  a  horizontal  coordinate;  and  2)a  daily-changing  definition  of  Gulf 
Stream  flow  direction  based  on  the  shear  between  the  thermocllne  and 


2000  m  depth.  Time-series  of  daily-rotated  velocities  may  be  used  to 
calculate  empirical  orthogonal  functions  for  the  long-  and  cross-stream 
vertical  structures,  which  are  decoupled  and  are  respectively  barocllnlc 
and  barotropic.  Using  the  Inferred  horizontal  coordinate  one  can 
estimate  mass,  momentum  and  kinetic  energy  fluxes  for  four  individual 
events  when  the  entire  Stream  swept  by  the  mooring.  The  results  agree 
well  with  historical  data.  Bryden's  (1980)  method  has  been  used  to 
calculate  vertical  velocities  from  the  temperature  equation;  the 
resulting  time-series  of  w  are  visually  coherent  throughout  the  water 
column  and  their  vertical  amplitude  structure  Is  reminiscent  of  that  for 
a  two- layer  system.  The  rms  vertical  velocities  are  large  (0( . 05  cm/s)), 
and  these  as  well  as  other  estimates  have  been  used  to  explore  the 
validity  of  the  quasi -geostrophic  approximation  at  the  mooring  site.  The 
Rossby  number  for  the  thermocllne  flow  Is  about  0.3,  and  for  the  deep 
flow  Is  _<  0.1. 

The  entire  data  set  may  also  be  used  to  construct  a  horizontal  and 
vertical  profile  of  velocity  in  the  Gulf  Stream,  from  which  a 
cross-section  of  the  mean  potential  vorticity  can  be  produced.  The 
latter  shares  many  common  feature  with  cross-sections  from  past  work  for 
a  nearby  site,  as  well  as  analogous  data  from  a  three-layer  numerical 
model,  thus  suggesting  that  they  are  robust  features  of  Gulf  Stream-like 
currents.  These  features  are.  In  particular,  a  strong  jump  from  low  to 
high  values  crossing  the  Stream  from  south  to  north;  and  a  change  In  the 
sign  of  the  potential  vorticity  gradient  on  Isothermal  surfaces  for 
T  >  12*C. 


To  complement  the  analysis  of  the  observational  data,  a  set  of 
diagnostic  calculations  has  been  performed  on  an  eddy- resol ving  general 
circulation  model,  to  provide  a  complete  picture  of  the  kinetic  energy 
budgets  of  the  free  jet  and  its  environs.  It  is  found  that  the 
downstream  convergence  of  kinetic  energy  in  the  decelerating  jet  is 
balanced  primarily  by  an  ageostrophic  flow  against  the  pressure  gradient, 
which  in  turn  implies  some  conversion  of  kinetic  to  available  potential 
energy  in  the  region.  Energetic  analysis  of  the  observations  as  well  as 
the  numerical  data  suggests  barotropic  and  barocllnic  instabilities  may 
be  equally  important  to  the  kinetic  energy  budgets  in  the  Stream. 

Because  there  is  but  one  mooring,  the  dynamics  governing  the 
fluctuations  remain  elusive.  Nonetheless,  a  kinematic  framework  is 
proposed,  which  is  consistent  with  the  data  and  accounts  for  a  variety  of 
unusual  features  that  arise  in  the  original  analysis  (for  example, 
distinct  asyimte tries  in  the  four  Gulf  Stream  crossings,  and  the  rather 
large  vertical  velocities).  It  is  suggested  that  the  data  we  are  now 
capable  of  collecting  is  proffering  fundamentally  new  attributes  of  the 
Gulf  Stream,  which  must  be  included  and  accounted  for  in  future 
theoretical  work. 
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Chapter  1.  Introduction 


Perhaps  the  most  remarkable  feature  of  the  general  circulation  of 
the  world  oceans  is  the  persistence  of  strong  western  boundary  currents 
closing  the  subtropical  gyre  circulations  of  each  ocean  to  the  west.  The 
Gulf  Stream  in  the  North  Atlantic,  the  Brazil  Current  in  the  South 
Atlantic,  the  Kuroshio  in  the  North  Pacific  and  the  Agulhas  In  the  Indian 
Ocean  are  all  narrow,  energetic  swiftly  flowing  currents  carrying  large 
mass  transports.  Indeed,  these  boundary  currents  are  responsible  for 
transporting  interior,  wind-driven  equatorward  flow  poleward  to  maintain 
mass  continuity;  presumably  a  large  fraction  of  the  energy  and  vorticlty 
input  by  wind  must  be  dissipated  in  these  currents  as  well.  They  are 
clearly  an  Integral  part  of  the  general  circulation,  and  until  we  fully 
understand  what  governs  their  behavior,  we  cannot  claim  to  understand  the 
ocean's  general  circulation  fully. 

The  Gulf  Stream  was  first  described  over  450  years  ago  (Stommel , 
1950,  traces  the  development  of  ideas  and  observations  in  the  Stream), 
and  since  then  has  been  observed  and  monitored  In  a  number  of  ways. 

While  it  Is  probably  the  most  heavily  documented  feature  in  the  deep 
North  Atlantic,  it  is  also  terribly  complicated,  and  continues  to  defy 
complete  understanding.  Part  of  the  problem  lies  In  the  type  of 
observations  that  have  been  feasible  In  the  Gulf  Stream.  The  earliest 
data  came  from  ship  drift  measurements,  which  served  to  outline  the 
general  currents  associated  with  the  surface  flow  along  the  coast  and 
eastward  into  the  North  Atlantic.  Once  it  was  recognized  that  a 
temperature  front  was  a  feature  of  the  Stream,  temperature  measurements 


were  used  to  define  the  surface  path  of  the  current.  Early  versions  of 
floats  and  drifters  were  employed  as  well.  Around  1900  the  geostrophlc 
relationship  was  employed  to  determine  oceanic  currents  (Stommel,  1950). 
In  more  recent  years  the  surface  thermal  front  has  been  monitored  by 
satellite,  which  has  the  advantage  of  producing  continuous  long 
time-series  of  front  position,  but  which  has  the  disadvantage  that  clouds 
often  obscure  the  front,  severely  reducing  the  amount  of  data  collected, 
particularly  in  winter. 

To  understand  the  dynamics  of  this  energetic  current,  however,  one 
must  know  something  of  the  three-dimensional  structure  of  Its  velocity 
field,  and  how  it  varies  with  time.  Hydrographic  surveys  such  as  Gulf 
Stream  '60  (Fugllster,  1963)  have  served  to  define  the  barocllnlc 
structure  rather  well,  but  since  that  study  had  but  a  few  float 
trajectories  to  aid  In  choice  of  a  reference  level  for  geostrophlc 
velocities,  the  barotropic  aspects  of  the  flow  remained  largely  unknown. 
Moreover,  a  hydrographic  survey  Is  at  best  a  snapshot  of  the  flow  In 
time.  The  deployment  of  large  numbers  of  floats  over  the  past  twenty 
years  has  helped  In  describing  gross  features  of  the  Gulf  Stream's  time 
variability,  --  e.g.,  how  meanders  and  eddies  affect  the  eddy  kinetic 
energy  patterns  In  the  North  Atlantic  (Richardson,  1983)  --  but  because 
there  Is  little  control  over  the  floats  after  deployment,  it  has  been 
virtually  Impossible  to  learn  how  the  barocllnlc  structure  Itself  varies 
with  time.  Clearly,  this  aspect  must  be  monitored  by  long-term  fixed 
arrays  with  Instruments  throughout  the  water  column,  but  until  recently, 
mooring  technology  was  Incapable  of  successfully  deploying  such  moorings 


in  a  current  as  strong  as  the  Gulf  Stream.  Now  that  the  technology  is 
available,  we  need  to  proceed  carefully  in  interpreting  the  data  obtained 
from  such  moorings,  particularly  when  there  is  data  from  but  a  single 
mooring,  as  is  the  case  here. 

Development  of  theoretical  Gulf  Stream  investigations  began  with 
explanations  of  a  steady  structure,  primarily  determined  by  either 
friction  (Stommel ,  1948;  Munk,  1950)  or  inertial  effects  (Fofonoff, 

1954).  Attempts  to  predict  the  time-varying  path  of  the  free  eastward 
Gulf  Stream  from  upstream  inlet  conditions  began  with  Warren  (1963),  and 
later  continued  with  a  model  by  Robinson  and  Niiler  (1967).  Various 
types  of  instability  models  have  been  used  to  explain  the 
time-variabil ity  of  the  Stream.  Orlanski  (1969)  used  a  two-layer  linear 
model  with  different  bottom  topographies  to  predict  time  and  space  scales 
for  instabl il i ties  of  the  Stream  on  the  continental  shelf  and  in  deeper 
water  over  the  continental  rise.  Luyten  and  Robinson  (1974)  and 
Robinson,  Luyten,  and  Flierl  (1975)  have  discussed  the  long-wavelength 
instabilities  of  a  thin,  quasi-geostrophic  meandering  jet  (assumed  to 
move  coherently  from  top  to  bottom).  Talley  (1982)  has  used  a 
horizontally  unbounded,  two-layer,  linear  model  to  examine  the  radiation 
of  energy  by  instabilities  away  from  various  jet  configurations  meant  to 
approximate  the  Gulf  Stream.  In  the  past  decade,  great  progress  has  been 
made  numerically,  as  well.  Holland  and  Lin  (1975)  were  the  first  to  run 
a  numerical  model  of  a  baroclinic  ocean  in  which  the  horizontal 
resolution  was  small  enough  (on  the  order  of  the  deformation  radius),  and 
the  viscosity  was  low  enough,  to  demonstrate  that  mesoscale  eddies  are 
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generated  spontaneously  in  a  steadily-forced  model.  The  numerical 
simplification  involved  in  introducing  quasi-geostrophic  eddy-resolving 
general  circulation  models  (Holland,  1978)  led  to  extensive  parameter 
studies  of  the  models  (Holland  and  Haidvogel ,  1980)  as  well  as 
investigations  of  the  stability  properties  of  the  model  jets  (Haidvogel 
and  Holland,  1978) . 

The  interpretaion  that  will  unfold  in  succeeding  chapters  is  seated 
in  the  body  of  knowledge  accumulated  from  these  theories  and 
observations.  Two  approaches  may  be  adopted  for  interpreting  the  data: 
in  the  first,  the  mean  flow  at  the  mooring  site  is  assumed  to  "define" 
the  Gulf  Stream  there.  In  the  second,  the  well  documented  baroclinic 
structure  of  the  Stream  is  used  to  infer  a  daily-changing  "Streamwise 
direction  of  flow."  That  is,  the  Stream  is  recognized  as  a  permanent 
front  that  may  change  its  orientation,  with  quasi -permanent  attributes 
such  as  strong  vertical  shear  of  the  flow  velocity.  Recent  work  by  Johns 
(1985)  demonstrates  that  the  primary  mode  of  displacement  for  the  Gulf 
Stream  thermocline  is  the  simple  translation  of  a  coherent  feature;  this 
result  is  crucial  in  determining  some  horizontal  structure  from  a  single 
mooring.  The  basic  tenet  of  this  thesis  is  that  such  a  well-defined 
feature  exists  and  that  its  description,  in  a  time-averaged  sense, 
differs  from  the  mean  measured  flow  at  the  mooring  site.  In  particular, 
the  analysis  concentrates  on  developing  a  description  of  this  feature  and 
enumerating  the  distinctions  between  the  two  possible  interpretations  of 
the  Gulf  Stream. 

Chapter  2  first  presents  the  vertical  structure  of  the  flow  at  the 
mooring  site  and  places  it  in  the  context  of  past  work.  Next  the 
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machinery  is  developed  to  extract  horizontal  information  from  the  data 
set,  allowing:  1)  estimates  of  mass,  momentum,  and  kinetic  energy 
transports;  2}  calculation  of  vertical  velocities  from  the  temperature 
equation;  and  3)  an  assessment  of  the  quasi-geostrophic  approximation  at 
the  mooring  site. 

Little  has  been  said  thus  far  of  the  downstream  changes  in  the  Gulf 
Stream,  but  that  is  an  area  in  which  little  progress  has  been  made.  The 
average  downstream  changes  in  potential  vorticity  structure  of  the 
Stream,  for  example,  are  so  gradual  that  they  are  swamped  by  local 
variability.  Fofonoff  and  Hall  (1983)  documented  downstream  changes  in 
Gulf  Stream  transports  of  mass,  momentum  and  kinetic  energy  ,  but  had 
difficulty  in  accounting  for  the  changes.  Long  time-series  at  several 
downstream  positions  are  necessary  to  assess  energetic  balances  in  the 
vicinity  of  the  Stream,  and  at  this  point  such  data  are  available  only 
from  numerical  simulations  of  oceanic  circulation.  In  Chapter  3  the 
energety  budgets  are  discussed  for  the  numerical  Gulf  Stream  analog  and 
its  environs  in  a  Holland  (1978)  model,  assessing  the  relative  importance 
of  various  mechanisms  in  accelerating  and  decelerating  the  jet,  as  well 
as  their  roles  in  the  inertial  recirculations  and  the  primarily 
wind-driven  portions  of  the  domain.  It  should  be  pointed  out  that  this 
analysis  resorts  to  the  first  rather  than  second  Gulf  Stream  flow 
definition  described  above:  that  is,  the  time-averaged  Eulerian  flow  in 
the  jet  region  is  assumed  to  be  the  time-averaged  jet. 

Chapter  4  is  an  extension  of  the  classical  development  of  the  data, 
as  a  cross-section  of  potential  vorticity  in  the  Stream  is  constructed 


and  described.  As  previously  mentioned,  the  structure  is  far  from 
simple.  However,  comparisons  with  other  data  and  with  numerical  model 
data  point  up  the  salient  features  in  cross-sections  of  potential 
vorticity  in  the  observed  Gulf  Stream  and  its  analog  in  the  numerical 
model.  These  are  evidently  robust  features,  which  have  immediate 
implications  for  the  dynamics  of  the  flow. 

In  Chapter  5,  an  effort  is  made  to  analyze  the  energetic  budgets  at 
the  mooring  site  for  comparison  with  past  work  and  with  the  numerical 
analysis  of  Chapter  3.  Then  a  kinematic  interpretation  of  the  flow 
regime  is  explored,  which  is  consistent  with  the  data  and  accounts  for 
some  of  the  more  unusual  results  described  in  Chapter  2.  Finally,  some 
suggestions  are  made  for  the  future  directions  Gulf  Stream  observational 
work  might  profitably  take. 


Chapter  2.  Profiling  the  Gulf  Stream  with  a  current  meter  mooring 

2.1  Introduction 

The  GUSTO  mooring,  deployed  for  one  year  at  37*37'  N,  68*00'  W,  the 
mean  position  of  the  Stream  at  that  longitude  (Halliwell  and  Mooers, 
1983),  represents  the  first  opportunity  to  examine  long  time  series  of 
current  measurements  throughout  the  water  column  in  the  Gulf  Stream. 
Currents,  temperature,  and  pressure  were  recorded  at  the  nominal 
deployment  depths  of  400,  700,  and  1000m,  while  current  and  temperature 
only  were  measured  at  2000  and  4000  m.  The  only  missing  data  is  at  the 
middle  instrument  (1000  m) ,  where  the  VACM  stopped  working  after  64  days, 
but  resumed  after  56  more  days.  Although  only  one  mooring  was  deployed, 
yielding  only  vertical  resolution  of  the  flow,  a  remarkable  amount  of 
horizontal  structure  can  be  inferred:  as  the  Stream  meanders  back  and 
forth  across  the  mooring,  the  temperature  and  pressure  measurements  at 
the  uppermost  instruments  can  be  used,  in  conjunction  with  historical 
data,  to  determine  how  far  north  or  south  of  the  current  axis  the  mooring 
is.  Thus,  the  data  profiles  the  horizontal  structure  of  the  Stream  in 
time.  Inspection  of  the  National  Weather  Service  analyses  of  satellite 
data,  indicating  the  approximate  surface  expression  of  the  Gulf  Stream, 
shows  that  the  mooring  was  in  the  Stream  58  percent  of  the  time,  in  the 
Slope  Water  to  the  north  12  percent  of  the  time,  and  Sargasso  water  to 
the  south  30  percent  of  the  time. 

Because  the  strong  currents  tilted  the  mooring  --  average  pressure 
at  the  top  instrument  was  498  <ft>ar,  with  a  minimum  of  433  dbar  and  a 


maximum  of  680  dbar  --  temperatures  and  velocities  were  interpolated  or 
extrapolated  to  intermediate  standard  pressures  of  575  dbar,  875  dbar, 
and  1175  dbar  to  create  time  series  that  could  be  used  in  a  consistent 
manner  (Raymer,  Spencer,  and  Bryden,  1984).  Even  a  cursory  inspection  of 
the  records  reveals  that  while  the  mean  velocity  vector  at  575  dbar  was 
directed  nearly  due  east,  instantaneously  the  flow  was  usually  north  or 
south  of  east.  It  has  been  pointed  out  that  there  are  at  least  two  ways 
to  discuss  a  "mean"  Gulf  Stream,  either  as  the  Eulerian  average  flow  in  a 
particular  region  or  as  a  discrete  feature  that  may  change  its  position 
or  orientation  while  retaining  certain  fundamental  characteristics.  The 
former  is  a  more  traditional  approach  and  is  discussed  first,  for 
comparison  with  past  work.  The  second  approach  is  more  fruitful  for 
describing  a  meandering,  frontal  jet  like  the  Gulf  Stream. 

2.2  Mean  statistics  of  the  flow 

All  the  data  from  the  mooring  was  low-pass  filtered  with  a  24-hour 
Gaussian  filter,  then  subsampled  daily,  to  provide  time  series  of  360 
daily  values.  Table  2-1  shows  the  record-length  mean  east  and  north 
velocities  and  temperatures  at  the  five  standard  depths.  Also  shown  are 
the  variances  of  these  quantities.  In  spite  of  the  often  large  deviation 
from  East  of  the  current  direction  at  575  dbar,  the  mean  velocity  there 
is  directed  essentially  due  East.  It  is,  however,  much  smaller  than  the 
maximum  speeds  recorded  there,  which  are  well  over  100  cm7s.  Mean 
northward  velocities  throughout  the  water  column  are  extremely  small  and 
are  nearly  barotropic,  though  on  a  daily  basis  there  may  be  considerable 


shear  in  that  direction  (Figure  2.1a).  Since  v  is  barotropic  while  u  is 
baroclinic,  the  mean  velocity  vectors  turn  cyclonically  with  increasing 
depth,  implying  a  mean  downward  vertical  velocity  (Bryden,  1980).  The 
mean  velocity  at  4000  m  is  directed  nearly  along- isobath,  which  is  about 
7o°  true  (bottom  depth  is  4688  m) ,  and  the  average  zonal  velocity  there 
is  westward  rather  than  eastward.  This  deep  flow  is  reminiscent  of  what 
Luyten  (1977)  found  in  his  "upper  rise"  regime  at  70°W,  with 
along-isobath  flow  directed  mostly  westward.  The  mean  temperature  at  575 
dbar  is  very  close  to  that  associated  with  the  Gulf  Stream  axis,  as 
discussed  in  Section  2.3,  suggesting  that  this  was  indeed  the  average 
location  of  the  Stream  for  the  year. 

The  eddy  kinetic  energy  is  surface  intensified,  as  Richardson  (1983) 
found  in  constructing  a  vertical  section  of  EKE  (eddy  kinetic  £nergy) 
along  55°W  from  drifter,  float,  and  current  meter  data.  In  that  section, 

however,  the  values  at  575  db  in  the  Gulf  Stream  region  are  roughly  500 

2-2  2-2 
cm  s  ,  only  about  half  that  of  around  1050  cm  s  from  Table 

2-1.  The  surface  intensification  appearing  in  both  velocity  components 

indicates  that  barocl  ini  city  associated  with  the  Gulf  Stream  appears  in 

north  as  well  as  east  velocities.  Below  1000  m  or  so,  the  EKE  values 

decay  very  little  with  depth  and  are  very  nearly  equi partitioned  between 

the  two  velocity  components.  Richardson's  55  W  values  do  not  decay  as 

2  -2 

rapidly  with  depth:  at  2000  m  he  displays  a  value  of  136  cm  s  , 

2  -2 

compared  to  108  cm  s  from  Table  2-1,  but  his  values  below  1500  m  or 

so  under  the  Stream  axis  remain  fairly  constant  at  around  130-140 
2  -2 

cm  s  .  Schmitz  (1984)  found  similar  results  for  abyssal  eddy 
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Figure  2.1.  a)  Time  series  of  east  (top)  and  north  (bottom)  velocities 
at  five  standard  depths;  b)  Flow  direction,  defined  as  direction  of  shear 
between  400  and  2000  m  instruments;  c)  Time  series  of  along-  and 
cross- stream  velocities  (see  text)  at  five  standard  depths.  At  top  is 
time-series  of  T575,  temperature  corrected  to  575  <tt>ar  (see  text). 
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kinetic  energy  in  particular,  that  its  variation  with  depth  is  much  less 

in  deep  water  than  shallow;  and  that  abyssal  (4000  m)  EKE  displays 

similar  scales  of  variation  in  the  zonal  and  meridional  directions.  He 

a  7  7 

found  actual  values  near  the  Gulf  Stream  at  70  W  of  104  cm  s  at 

o  -  ? 

4000  m,  about  twice  the  value  of  58  cm  s  form  Table  2-1,  but  less 

7  -7 

than  Richardson's  value  of  138  cm  s’  .  Schmitz  also  points  out  that 
there  is  a  gap  in  data  coverage  between  55  W  and  70  W,  but  that  a  maximum 
in  eddy  kinetic  energies  might  be  expected  there.  In  summary,  the  values 
in  Table  2-1  seem  to  fit  in  well  with  other  documented  values  in  this 
region,  both  qualitatively  and  quantitatively. 

2.3  Vertical  structure  of  the  "average  Gulf  Stream" 

In  order  to  obtain  a  description  of  the  Gulf  Stream  as  a  discrete 
feature  from  the  flow  at  the  current  meter  site,  it  is  necessary  to 
define  what  is  meant  in  referring  to  the  Stream.  As  the  Gulf  Stream 
meanders  and  changes  direction,  a  significant  part  of  the  "along-stream" 
flow  may  be  contained  in  the  northward  component  of  velocity;  hence  a 
definition  of  what  direction  the  Stream  is  flowing  at  any  given  time  is 
required,  to  determine  how  large  the  along-stream  flow  is.  The  direction 
of  the  shear  between  the  measured  current  at  the  uppermost  (400  m) 
Instrument  and  that  at  the  2000  m  instrument  has  been  chosen  as  the 
definition  of  the  along-stream  flow  direction  for  several  reasons.  While 
the  velocity  vector  may  rotate  with  depth  in  the  current  records,  that 
can  be  due  to  a  small  barotropic  flow  superimposed  on  the  primarily 
baroclinic  jet:  strong  vertical  shear  has  long  been  recognized  as  a 


signature  of  Gulf  Stream  flow.  The  time  series  of  east  and  north 
currents  at  the  five  Instruments  demonstrate  that  when  the  Gulf  Stream  is 
present  (as  evidenced  by  the  current  speeds  and  the  temperature  at  575 
dbar),  it  appears  to  penetrate  all  the  way  to  the  bottom  instrument 
(Figure  2.1a);  at  other  times,  for  example  the  beginning  of  April,  there 
is  a  bottom-intensified  westward  flow  in  the  deep  water.  To  avoid  any 
problems  with  such  reversals  at  4000  m,  the  direction  of  shear  between 
that  and  the  400  m  instrument  was  not  chosen.  Furthermore,  the  choice  of 
2000  m  is  a  classical  reference  level  of  no  motion  used  in  the  North 
Atlantic.  With  such  a  definition,  the  along-stream  direction  changes  day 
by  day  as  shown  in  Figure  2.1b. 

The  time  series  of  along-  and  cross-stream  velocities  are  shown  in 
Figure  2.1c,  along  with  T575  (temperature  corrected  to  575  db).  A 
monotonic  increase  or  decrease  in  Tg75  signals  the  passage  of  the  Gulf 
Stream  across  the  mooring  site:  four  clear  examples  occur  in  March, 

June,  late  August/early  September,  and  later  September.  The  Gulf  Stream 
also  occupied  the  mooring  site  for  long  periods  when  it  did  not  sweep 
completely  across,  such  as  from  November,  1982  to  March,  1983. 

Comparison  of  Figures  2.1a  and  c  shows  that  the  third  event  is  an 
excellent  argument  for  defining  a  daily  along-stream  direction:  the  flow 
then  was  to  the  northwest  (Figure  2.1b),  yet  according  to  T575  and 
corroborating  evidence  from  the  NOAA  satellite  pictures,  this  flow  was 
indeed  the  Gulf  Stream.  In  fact,  the  top-to-bottom  coherence  of  Gulf 
Stream  flow  is  more  apparent  in  Figure  2.1c  than  2.1a.  The  baroclinicity 
of  the  along-stream  flow  is  obscured  due  to  the  changing  velocity  scales 


with  depth,  which  were  chosen  so  that  the  4000  m  velocities,  e.g.,  would 
in  fact  be  discernible.  However,  note  that  all  the  scales  on  the 
cross- stream  velocities  are  the  same,  displaying  the  remarkably 
barotropic  character  of  the  cross-stream  velocities:  this  feature  would 
never  have  been  isolated  in  strictly  east/north  coordinates. 

To  quantify  the  baroclinic  and  barotropic  structures  of  the  along- 
and  cross-stream  velocities,  empirical  orthogonal  functions  (EOF's)  for 
their  vertical  structures  were  computed.  Time  series  from  the  standard 
depths  were  used,  and  EOF's  were  computed  using  data  from  the  four 
complete  (360-day)  records.  EOF's  may  be  computed  for  the  separate 
velocity  components  as  well  as  for  the  vector  velocity  as  a  function  of 
depth.  If  there  is  a  strong  coupling  between  the  two  components,  the 
structure  of  the  latter  may  be  considerably  different  thatn  the  structure 
obtained  by  adding  the  EOF  structures  for  the  individual  components. 

EOF's  were  computed  both  ways  and  compared,  but  no  such  differences  were 
found,  indicating  that  there  is  little  coupling  between  along-  and 
cross-stream  velocity  components.  The  correlation  coefficient  for  the 
two  amplitude  time  series  attains  a  maximum  value  of  C  =  0.1,  when  the 
cross-stream  series  is  not  lagged  at  all.  For  an  estimated  30  degrees  of 
freedom,  C  must  be  _>  0.3  for  a  significant  correlation  at  the  95  percent 
confidence  level.  Only  calculations  for  the  separate  components  are 
described  here.  Note  that  when  data  is  available  at  N  points  in  the 
vertical,  N  EOF's  will  be  computed. 

When  data  from  the  575,  875,  2000,  and  4000  dbar  records  are  used, 
92.4  percent  of  the  variance  in  cross-stream  velocity  is  accounted  for  by 
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the  first  mode,  with  the  vertical  structure  shown  in  Figure  2.2a.  Notice 
how  nearly  barotropic  the  mode  is.  The  second,  third,  and  fourth  modes 
contain,  respectively,  5.9,  1.7,  and  0.03  percent  of  the  variance  and  are 
thus  of  little  importance  in  comparison  to  the  first  mode.  In  the 
along-stream  direction,  the  first  calculated  EOF  accounts  for  89.8 
percent  of  the  variance,  and  has  the  strongly  baroclinic  vertical 
structure  shown  in  Figure  2.2b.  There  is  no  reversal  of  the  amplitude  at 
4000  m,  indicating  that  the  deep  flow  in  these  coordinates  is 
statistically  in  the  same  direction  as  the  thermocline  flow.  The  second, 
third,  and  fourth  modes  contain  6.2,  3.7,  and  0.3  percent  of  the  energy, 
respecti vely. 

Using  hat  notation  to  denote  rotated  velocities,  the  along-  and 
cross-stream  velocities  may  be  expressed  in  terms  of  east  and  north 
velocities,  and  the  direction  a  of  the  flow  (see  Figure  2.6  for 
schematic) : 

^  A  A 

U  =  UCOSa  +  VSina  U  =  UCOSa  -  VSina 

V  =  VCOSa  -  usina  V  =  VCOSa  +  usina  (2-1). 

Since  the  time  mean  (denoted  by  (  ))  is  first  removed  when  calculating 


EOF's,  the  first  EOF  actually  yields  time  series  for 
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(2-2) 


where  m2+  m2+  m2  +  mjj  =  1,  n2+  n2+  n2+  nj  =  1,  and  Ja2(t)dt  = 

jb2(t)dt  =  nj,  the  first  eigenvalues.  Thus,  in  terms  of  east  and 


Figure  2.2.  a)Vertica1  structure  of  the  first  EOF  calculated  for  cross 


stream  velocity  (see  text  for  definition),  using  complete 
(360-day)  records  at  4  depths;  b)same,  but  for  long-stream 
velocity . 


north  velocities,  what  the  EOF's  actually  yield  is  the  following 
descri ption: 
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The  two  decoupled  modes  of  the  rotated  frame  combine  to  yield  the  more 
complicated  description  of  the  data  in  terms  of  east  and  north 
components.  The  actual  time  series  of  (u,v)  and  the  derived  series  from 
(2-3a,b)  compare  extremely  well. 


Depth 

(db) 

u  (cm/s) 

v(cm/s) 

(F7!— ) 

s2 

-k+ 2  ,cmz . 
v  z  ( - ) 

s2 

575 

51.003 

-.600 

918.60 

45.01 

875 

25.020 

-.448 

248.67 

56.55 

1175 

11.942 

-.027 

77.17 

45.09 

2000 

8.852 

-.638 

36.48 

44.67 

4000 

2.118 

-.147 

37.72 

75.15 

Table  2-2.  Record-length  statistics  for  velocities  in  rotated  co¬ 
ordinate  system  at  standard  depths. 
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Table  2-2  shows  the  mean  statistics  for  the  velocities  in  the 
rotated  coordinate  system.  The  mean  along-stream  velocities  are  stronger 
than  the  mean  east  velocities  throughout  the  water  column,  and  are  In  the 
same  direction  top  to  bottom.  The  cross-stream  velocity  variance  Is  very 
barotroplc,  and  the  total  EKE  throughout  the  water  column  is  less  for 
these  coordinates,  indicating  that  much  of  the  energy  calculated  in  Table 
2-1  is  associated  with  the  meandering  of  the  Gulf  Stream  rather  than 
changes  in  its  inherent  structure. 

In  summary,  the  successful  definition  of  an  along-stream  flow 
direction  leads  to  a  clean  description  of  the  flow's  vertical  structure 
at  the  mooring  site.  In  the  rotated  coordinates,  the  fluctuations,  as 
well  as  the  mean,  break  down  neatly  into  a  baroclinlc  along-stream  mode 
and  a  barotroplc  cross- stream  mode,  which  are  decoupled,  and  each  of 
which  contains  about  90  percent  of  the  variance.  Furthermore,  the 
structure  of  both  the  mean  and  fluctuating  along-stream  mode  show  no  flow 
reversals  at  depth,  Indicating  that  the  Gulf  Stream  as  defined  here  does 
Indeed  penetrate  to  the  bottom. 

2.4  Horizontal  structure  of  the  Gulf  Stream 

The  meandering  of  the  Gulf  Stream  back  and  forth  past  the  mooring 
site  suggests  that  a  horizontal  description  of  the  Stream  might  be 
deduced  from  the  mooring  data,  to  complement  the  vertical 
description.  A  scatter  plot  of  along-stream  velocity  versus 
temperature  at  575  dbar  (Fig.  2.4)  shows  that  the  former  Is  a  strong 
function  of  the  latter.  If  the  cross-stream  temperature  structure  at 


575  dbar  remains  fairly  constant  in  time  {though  it  may  meander 
about),  then  at  any  given  time  the  corrected  temperature  at  575  dbar 
ought  to  indicate  how  far  In  the  cross-stream  direction  the  mooring  is 
located  relative  to  the  axis  of  the  long-stream  flow.  Johns  (1985) 
has  recently  shown  that  the  primary  mode  of  displacement  of  the 
thermocline  in  the  Gulf  Stream  is  a  strict  translation  of  the 
isotherms,  rather  than  either  a  tilting  or  squeezing,  for  example. 
Thus,  this  assumption  of  fixed  cross-stream  structure  is  justified. 

The  much  longer  scales  of  variation  In  the  downstream  direction 
prevent  a  similar  approach  for  a  description  of  long-stream  flow 
structure,  which  will  henceforth  rarely  be  discussed.  The  question  is 
now  this:  how  can  this  notion  of  temperature  and  cross-stream 
distance  being  functions  of  one  another  be  quantified? 

Consider  Figure  2.3,  which  shows  a  scatter  plot  of  aT/ay  as  a 
function  of  T  at  575  db,  where  the  A  notation  refers  to  rotated 
coordinates.  With  thermal  wind  (pQfau/az  =  gap/ay)  and  the 
assumption  that  salinity  Is  a  function  of  temperature  alone  such  that 
dp  =  -aQdT  where  -aQ  =  3p/aT  +  <  a p/aS) (dS/dT ) ,  the  cross-stream 
temperature  gradient  at  575  dbar  may  be  obtined  from  the  measured 
velocities  and  pressures  at  the  400  and  700  m  Instruments.  Since 
rotated  velocities  were  used,  technically  (f  *  3a/at)  ought  to  be  used 
In  place  of  f  In  the  geostrophic  relation;  however  at  Is  generally 
several  orders  of  magnitude  less  than  f,  and  so  it  can  be  Ignored. 

The  points  were  divided  Into  two  categories,  those  for  which  Tg^g  > 
13*C  or  Tg7g  <  13*C,  and  least  squares  linear  fits  were  obtained 


dT/d 


for  T*,  which  are  also  shown.  The  exact  forms  for  the  two  lines  are 
T  <  13°C.:  il  =  0.08849-^  -  0.02438  km-1x  T  (2-4a); 


>  13*C.:  4-  =  -1 .0043— r~  +  0.( 


.-1. 


T  >  13'C.:  -4-  =  -1.0043— r~  +  0.06025  km_1x  T  (2-4b) 
ay 

If  y  =  0  for  T  =  13*C,  (2-4  a,b)  may  be  integrated  in  either  direction 
to  obtain  T  =  T(y),  and  inverted  to  obtain  y  =  y ^575)  •  Tlie  results 
are: 

T  <  13° C.:  T  =  3.6296°C.  ♦  9.3704°C.  e“-02438y/km  {2-5a ) 


-km  1 _.08849°C./km  -  .02438km_1T 

_  -  1  p - 


(2-6a) 


.02438  -0.22845  *C/km 

T  >  13°C.:  T  =  16.670"C.  -  3.6695°C.  e*06025y/kni  (2-5b ) 


y  =  ■■  -k-m-  1  n( 4.5427  -  ,27252(*C. )-1T)  (2-6b) 

.06025 

It  should  be  noted  that  at  depths  greater  that  575  dbar,  there  is  not 
such  a  strong  and  clear  dependence  of  the  shear  (and  hence  T^)  on 
temperature,  making  any  other  level  less  suitable  for  this  procedure 
than  is  the  575  dbar  level. 

The  scatter  plot  of  lig^  (along-stream  velocity  at  575  dbar) 
vs.  Tg7g  (Figure  2.4)  again  suggests  that  a  linear  least  squares  fit 
is  a  reasonable  description  of  the  data--  note  particularly  the  break 
in  functional  form  at  T  =  13*C.  The  results  of  the  linear  fit  are: 

T  <  13*C. : 

ug75  =  -40 .9 5 lcm/ s  ♦  11 .332cm/s/*C.  x  T  (2-7a) 

T  >  13°C. : 

Ug7g  =  479.35  cm/s  -  27.463cm/s/*C.  x  T 


(2-7b ) . 
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Combining  equations  (2-5)  with  (2-7)  gives  u  =  u(y ) : 

T575<  13  C: 

u(y)  =  u(y  =  0)  -  106.19  H  (e'*02438  y/km-  1) 


—  =  -.2589  x  10"4s_1  e'*02438  y/km 

ay 
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(2-8a ) 
(2-9a) 

(2-8b) 
(2-9b ) 


Unfortunately,  these  descriptions  yield  discontinuities  in  both  u  and 
au/ay;  however,  they  are  used  primarily  for  qualitative  reasons  and  the 
two  regimes  are  often  considered  separately.  An  immediate  observation 
from  these  results  is  that  the  warm,  anticyclonic  side  of  the  Stream  is 
actually  sharper  than  the  cold,  cyclonic  side  in  the  sense 

that  —  >  —  ,  in  fact  nearly  three  times  as 

3y  r575>13  C  l *T575<  13  C 

sharp.  This  result  contradicts  the  notion  presented  in  Fofonoff  and  Hall 
(1983)  that  the  northern  edge  of  the  Stream,  in  analogy  with  the 
two-layer  inertial  jet  model,  should  exhibit  a  much  more  sudden  jump  in 
velocity  than  the  warm  side.  It  may  be  that  the  structure  proposed  by 
Fofonoff  and  Hall  is  found  higher  up  in  the  water  column.  The  North  Wall 
of  the  Gulf  Stream,  defined  as  the  place  where  the  15*C  isotherm  is  at 
200  m  depth  (Fuglister,  1963)  has  long  been  used  as  an  indicator  of  Gulf 
Stream  posftion,  but  clearly  that  is  well  above  the  depth  range 
detectable  by  the  mooring.  Analysis  of  the  Pegasus  sections  at  73*  W  by 


Halkin  and  Rossby  (1985)  shows  that  instantaneously  as  well  as  on  average 
the  horizontal  velocity  shear  is  much  stronger  on  the  cold  than  on  the 
warm  side  in  the  upper  few  hundred  meters;  even  at  600  m  depth,  though 
less  pronounced,  the  same  trend  is  definitely  there.  The  velocity 
section  constructed  by  Warren  and  Volkmann  (1969),  on  the  other  hand, 
gives  a  general  impression  that  au/sy  is  stronger  on  the  cold  side,  but 
at  600  m  depth  the  shear  appears  to  be  nearly  symmetric  about  the  maximum 
along-stream  velocity  at  that  depth.  Thus,  more  information  in  the  upper 
500  m  of  the  water  column,  taken  concurrently  with  deeper  information, 
will  be  required  to  resolve  this  question. 

2.4.1  Horizontally  integrated  fluxes 

In  Hall  and  Bryden  (1985),  an  average  velocity  profile  of  the  Gulf 
Stream  at  the  mooring  site  has  been  constructed,  based  on  the  assumption 
that  the  cross-stream  temperature  structure  is  constant  (Figure  2.5). 

From  the  profile,  an  average  transport  value  of  103  Sv  can  be 
calculated.  However,  it  is  also  instructive  to  examine  individually  and 
compare  the  four  events  when  the  Gulf  Stream  swept  by  the  mooring,  for 
the  differences  in  the  events  are  illuminating. 

Time  may  be  used  parametrically ,  with  u  =  u(t)  and  y  =  y(T(t)),  to 
integrate  streamwise  fluxes  horizontally  as  well  as  vertically.  The 
primary  difficulty  in  calculating  integrated  transports  is  defining  the 
"edges"  of  the  Stream.  This  has  long  been  a  troublesome  problem.  One 
example  of  the  type  of  problem  that  can  arise  was  discussed  in  Fofonoff 
and  Hall's  (1983)  treatment  of  the  Gulf  Stream  '60  data  (Fuglister, 

1963).  The  geostrophically  calculated  velocities  at  Section  I  exhibit 


Along- Stream  Velocity  (cms-’J 


profile  was  obtained  In  Hall  and  Bryden  (1985). 
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closed  circulations  north  and  south  of  the  main  front:  should  the 
contributions  from  these  circulations  be  included  in  the  transport  One 
way  of  delineating  the  edges  is  to  look  for  the  point  at  which  either  the 
transport  per  unit  width,  or  the  velocity  anywhere  in  the  water  column, 
goes  to  zero.  While  that  approach  may  be  successful  for  the  north  edge, 
velocities  at  the  southern  edge  tend  to  decay  very  slowly.  Moreover, 
Schmitz(1980)  has  identified  what  is  apparently  a  weakly  depth-dependent 
recirculation  regime  about  200-300  km  south  of  the  Gulf  Stream  axis  at 
55°W,  and  has  suggested  that  it  might  be  a  recirculation  increasing  the 
total  Gulf  Stream  transport.  In  this  analysis,  the  problem  manifests 
itself  in  that  while  the  relationship  between  aT/ay  and  T  is  less  tight 
for  temperatures  far  from  the  axis  of  the  flow,  so  that  a  small  change  in 
T575  leads  to  a  large  change  in  y,  velocities  there  may  be  too  strong 
to  be  excluded  from  the  integration.  It  is  best  to  approach  each  event 
individually,  with  the  restriction  that  only  data  in  the  temperature 
range  5°C  _<  T575  <_  16.5°C.  are  to  be  considered.  Since  none  of  the 
passages  spans  this  full  range,  temperature  and/or  velocity  extrema  are 
used  to  define  the  integration  limits  for  each  event. 

Since  the  data  only  extended  upward  as  far  as  575  dbar,  and  since 
much  of  the  transport  is  known  to  occur  above  that  depth,  some 
extrapolation  scheme  was  required  to  extend  velocities  to  the  surface. 
Different  schemes  were  tested  on  hydrographic  data  from  Gulf  Stream  '60 
at  68.5°W  to  determine  which  one  best  reproduced  the  actual  shear  in  the 
upper  water  column.  If  (u^,P4)  are  the  measured  along-stream 
velocities  and  pressures  at  the  400  and  700  m  instruments  respectively,  a 
value  of  u^  =  u^  +  1/2  (u^  -  u7)  was  assigned  to  the  level  P^  =  2P^  - 
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P i  (linear  extrapolation),  and  this  value  was  considered  to  be  constant 

to  the  surface.  At  the  bottom,  the  velocity  was  assumed  to  decrease 

linearly  from  its  value  at  4000  m  to  zero  at  the  bottom  depth  of  4688  m. 

The  actual  measured  velocities  and  pressures  were  used  since  these  are 

depth-integrations,  not  requiring  standard  levels. 

Table  2-3  gives  the  pertinent  information  for  the  four  events,  which 

divide  themselves  more  or  less  into  two  categories.  In  March  and  early 

September,  T575  was  decreasing,  indicating  that  the  translational 

cross-stream  velocity  of  the  flow  was  negative  (southward-like);  yet  the 

average  measured  cross-stream  velocities  for  those  two  events  were 

positive  (see  Fig.  2.1c,  bottom).  In  June  and  late  September, 

translational  velocities  were  positive,  and  so  were  the  average  measured 

cross-stream  velocities.  The  maximum  velocities  at  575  db  for  the  former 

two  events  are  considerably  lower  than  for  the  latter  two  (about  20-30 

cm/s),  yet  mass  transport  is  correspondingly  lower  only  for  the  March 

event,  if  the  narrower  width  of  the  early  September  passage  is  taken  into 

account.  A  glance  at  Figure  2.1c  shows  that  strong  along-stream 

velocities  at  4000  m  in  early  September  are  responsible  for  the 

difference.  The  values  in  Table  2-3  are  similar  to  other  estimates  near 

this  same  longitude.  Fuglister  (1963)  used  hydrographic  data  and  an 

assumption  of  zero  bottom  velocity  to  calculate  a  Gulf  Stream  transport 
«  6  3 

at  68.5  W  of  136  x  10  m  /s.  Using  velocity  measurements  from  floats 

along  with  hydrographic  data  at  38°N,  69°W,  Warren  and  Volkmann  (1968) 

6  3 

estimated  a  Gulf  Stream  transport  there  of  101  x  10  m  /s,  with  an 
average  bottom  velocity  of  1.5  cm/s  in  the  same  direction  as  the  surface 
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Integrated  momentum  and  kinetic  energy  fluxes  have  also  been 
calculated  and  are  listed  in  Table  2-3,  For  comparison,  values  from 
Section  I  of  Gulf  Stream  *60  (Fofonoff  and  Hall,  1983)  are  also 
presented.  The  June  event  most  resembles  the  flow  that  was  surveyed  In 
Gulf  Stream  60,  In  terms  of  total  fluxes.  The  extrapolation  scheme  to 
the  surface  for  the  GUSTO  data  may  be  somewhat  conservative,  yet  it  Is  In 
the  upper  water  that  the  greatest  contributions  to  momentum  and  kinetic 
energy  fluxes  arise.  That  may  explain  why  the  mass  transport  Is 
estimated  reasonably  well,  while  the  momentum  and  kinetic  energy 
transports  seem  to  be  underestimated.  In  conclusion,  the  transport 
values  obtained  here  are  reasonable  for  data  from  a  single  mooring,  but 
It  would  have  been  helpful  to  have  more  Information  In  the  upper  part  of 
the  water  column . 


2.4.2  Potential  vortlclty  at  575  dbar 

Ertel's  theorem  (see  Pedlosky,  1979,  for  a  good  discussion)  shows 
that  water  parcels  conserve  their  potential  vortlclty  In  the  absence  of 
forcing,  dissipation,  and  mixing,  where  potential  vortlclty  q  Is  defined 
as 


q  =  (  2fl+  ur)*Tx 


(2-10). 


Here  2  Is  planetary  vortlclty;  ur  =  Vxu  Is  relative  vortlclty;  p  Is 
density;  and  x  Is  any  conservative  scalar  property  of  the  fluid.  Taking 
x  to  be  potential  temperature  and  p  approximately  constant  yields: 

q  =  (f  +  vx  -  uy)ez  ♦  (h  +  uz  -  wx)Ty  ♦  (wy  -  vz)Tx 
where  (f,h)  =  (vertical,  horizontal)  components  of  planetary  vortlclty. 
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Since  terras  are  Implicitly  considered  In  the  rotated  frame.  In  the  Gulf 
Stream,  the  following  Inequalities  may  be  applied: 

Og  «  u«;  Ta  «  Ta;  h  «  uz; 

vz  «  uz;  wa  «  uz;  wa  <<  vz 

to  get  a  good  approximation  of  q  as 


q  =  (f  -  up)ez  +  uzT*  =  (f  -  fi*)9z  -  ^o  t| 


using  thermal  wind  to  evaluate  T*.  The  analytic  forms  for  all  terms  at 
575  dbar  as  functions  of  y  or  Tgyg  may  be  used  to  obtain  a  qualitative 
picture  of  the  cross- stream  potential  vortlclty  structure  at  575  dbar 
across  the  Gulf  Stream.  Notice  that  in  quasi -geostrophic  theory 
the  term  -UaTz  would  be  order  Rossby  number  smaller  than  f  x  ©z,  while 
uzT»  would  be  order  Rossby  number  squared  smaller  than  f©z.  The  latter  would 
certainly  be  Ignored  in  an  estimate  of  q,  and  In  areas  of  large  scale 
flow,  so  Is  the  former.  Table  2-4  gives  values  of  the  various  components 
across  the  Stream,  estimated  from  the  analytic  forms  for  u»,  «z,  and 
T*  at  575  dbar.  Although  fez  clearly  dominates  the  q  values,  the 
remaining  terms  make  important  contributions  near  the  Stream  axis,  where 
horizontal  and  vertical  gradients  are  strongest.  Other  Investigators 
(McDowell,  Keffer  and  Rhlnes,  1982;  McCartney,  1982)  have  pointed  out 
that  the  relative  vortlclty  contribution  ought  to  be  Included  in  an 
evaluation  of  potential  vortlclty  In  the  vicinity  of  strong  boundary 
currents  like  the  Gulf  Stream.  The  results  In  Table  2-4  suggest  that 
UaTz  and  uzT*  ought  to  be  retained  when  constructing  a  section  of 
potential  vortlclty  In  the  Gulf  Stream  (Chapter  4). 
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2.5  Governing  equations  in  rotated  coordinate  system 

Up  to  now,  the  notion  of  rotated  velocities  has  been  used  only 
mechanistically.  The  momentum  equations  for  the  rotated  system  are 
actually  quite  complicated  for  this  rotating  system,  even  with  the 
assumption  that  the  origins  (x,y)  =  (0,0)  and  (x,y)  =  (0,0)  remain  always 
at  the  mooring  site,  which  in  turn  implies  that  as  the  Stream  rotates,  it 
always  does  so  about  an  axis  fixed  at  the  site  (Fig.  2.6).  The  change  in 
the  Stream  orientation  from  day  to  day  is  generally  small,  but  away  from 
the  origin  the  rotation  leads  to  centrifugal  terms  which  make  this  system 
a  questionable  one  for  investigating  dynamics  away  from  the  mooring  site 
(difficult  with  just  one  mooring  in  any  case'.).  However,  the  rotation 
has  been  used  primarily  to  construct  a  picture  of  the  Gulf  Stream  from  a 
number  of  snapshots  taken  at  different  angles.  The  momentum  equations  at 
x=y  =  x=  y  =  0  are  presented  nonetheless.  The  transformations  Involve 
terms  like  u  =  u(x(x,y,t),y(x,y,t),t)  whence: 


-|^(ucosa  +  vsina)  +  -^-(ucosa  +  vsinc*)-|^  + 

A 

— (ucosa  +  vsina)!^- 

^  d  L 

3y 


(2-11). 


•  ••.'  -.1 


But  3x/3t  =  (-xsina  +  ycosa)at  =  0  at  x  =  y  =  0,  so  that  these  terms  do 
not  enter  into  the  momentum  equations  at  the  site.  The  resulting 


momentum  equations  (valid  only  at  the  site)  are: 

Ut  +  UU»  +  VU»  +  wGz  -  (f  +  at)v  =  -P* 


( 2- 12a ) 


Vt  +  UVJ  +  VV*  +  WVZ  +  (f  +  at)u  =  -Pp  (2- 12b ) . 

Thus,  thinking  in  terms  of  rotated  coordinates  at  the  site  adds  only  one 


V  A  S  \  ’  • 


•-V. 


term  to  each  equation,  and  {a^/fKO.08  during  the  whole  time  series. 
Notice  also  that  possible  curvature  of  the  Stream  path  has  been  excluded 
essentially  by  taking  a  4  a(x,y). 

The  vorticity  equation  at  the  origin  may  be  derived  as  follows:  at 
the  current  meter  site,  the  total  derivative  following  a  fluid  parcel 
retains  its  form,  i.e., 


i  +  ui  +  v  —  -  i-  +  i ;  i-  ♦  c  i- 

3t  u  ax  v  ay  at  3* 


(2-13)  . 


It  is  also  true  that  £=  vy  "  uw  =  v5  "  uO*  l’hus»  the  vorticity  equation 

x  y  x  y 


becomes 


(3  ♦  u  1-  ♦  V  i_)  (™  -  i“)  ♦  BV  =  fVZ 
aT  ax  ay'  ax  ay'  az 


(1_  +  5  1_  ♦  C  t*l  .  2±\  ♦  bv  -  f 

'at  u  >  v  a1  v  *  »'  p  az 

ax  ay  ax  ay 


(2- 14a) 


(2- 14b) . 


at  the  mooring  site.  In  equation  (2-14b)  £  can  be  approximated  by  -u», 
as  noted  before. 

Finally,  the  equation  for  heat  conservation,  using  standard 
notation,  is: 


il  ♦  u  il  ♦  v  il  ♦  w  2®  .  RHS 

at  ax  ay  az  '  KM:> 


(2-15) 


where  the  RHS  may  include  source  and  sink  terms  which  will  be  assumed 
negligible.  This  equation  assumes  the  same  form  in  rotated  coordinates, 
i.e.,  T^  +  uT*  +  vTa  +  w«z  =  0. 

With  data  from  just  one  mooring,  it  is  difficult  to  evaluate  most  of 


the  terms  in  these  equations,  in  order  to  determine  what  balances 
obtain.  Certain  attributes  of  the  flow  may  be  tested,  however.  It  is  a 


straightforward  matter  to  compare  the  local  time  derivatives  with  the 
Coriolis  terms  in  the  momentum  equations  to  obtain  some  idea  of  the 
strength  of  geostrophic  balance  at  the  mooring  site.  That  has  been  done 
for  both  the  east/north  and  along/cross- stream  velocities,  by  comparing 
the  rms  values  of  the  local  accelerations  with  those  of  the  associated 
Coriolis  accelerations,  and  assuming  that  the  advective  terms  are  of 
comparable  magnitude  to  (or  smaller  than)  the  local  accelerations.  Table 
2-5  shows  the  results  for  575  and  875  dbar.  There  is  little  question 
that  in  either  the  north  or  cross-stream  direction,  geostrophic  balance 
is  quite  strong,  when  data  from  the  whole  time  series  are  used,  the 
ratio  r  of  the  terms  at  575  (875)  dbar  is  .115  (.058)  for  the 
along-stream  direction  and  .032  (.032)  for  the  east  momentum  equation. 
When  data  from  a  single  strong  Gulf  Stream  event  are  used  (the  June 
event),  r  increases  to  .186  (.090)  for  along-stream  and  .057  (.058)  for 
east;  thus,  using  the  geostrophic  approximation  to  evaluate  temperature 
gradients  from  velocity  shear  should  introduce  errors  no  larger  than  20 
percent. 


2.5.1  Vertical  velocities  at  the  mooring  site 

Now  consider  the  balance  of  terms  in  the  heat  conservation 
equation.  Given  that  the  geostrophic  approximation  is  valid  to  within  20 
percent  at  the  mooring  site,  thermal  wind  yields: 


C  3U  „  3p  „  3T  . 

P0^  3Z  =  9  A="9ao  A» 

ay  ay 


f  3V  „  3p  _  3T 

pof  3Z  =  -9  a  9ao  a 

3X  3X 


(2-16) . 


The  quantity  aQ,  defined  before,  can  be  evaluated  from  historical  data 


-4  3o 

and  is  approximately  10  gm/cm  /  C,  with  variations  of  about  25 


Depth  Pi  recti  on 


Ratio 


Depth 

Di recti  on 

^ut  Vms 

(cm/s2) 

(fvVns 

575 

East 

6  .98x10*5 

2.19xl0'3 

575 

L- Stream 

6.88xl0‘5 

5.99xl0'4 

975 

East 

3.90xl0'5 

1.24X10'3 

875 

L-Stream 

3.90xl0'5 

6.70xl0'4 

Depth 

Direction 

(vtVms 

(cm/s2) 

^uVms 

575 

North 

6.45xl0‘5 

4.84xl0'3 

575 

X- Stream 

2.15xl0"5 

5.28xl0"3 

875 

North 

4.62xl0'5 

2.42xl0'3 

875 

X-Stream 

3.67xl0'5 

2 .63xl0”3 

.032 

.115 

.032 

.058 

Ratio 

.013 

.004 

.019 

.014 


Table  2-5.  Comparison  of  rms  local  and  Coriolis  accelerations  for 
momentum  equations  in  east/north  and  long/cross-stream  coordinates,  at 
575  and  875  dbar,  for  year-long  time  series. 
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percent  through  the  water  column.  Furthermore,  if  the  velocities  are 
expressed  as  u  =  R  cos/6,  v  =  R  sinj6.  where  R  is  the  speed  and 
i  =  tan'*(v/u)  is  the  angle  the  velocity  vector  makes  with  east,  the 
temperature  equation  (2-15)  may  be  rewritten: 


&  .  1°.  r2  „  |S  . 

3t  3  Z  3  Z 

Under  an  assumption  of  negligible  mixing  (i.e.,  RHS  is  small),  vertical 
velocities  may  be  calculated  from 


»  ■  -<S  ‘  r2  -n  >/i#  <2-17> 

J  o 

as  originally  pointed  out  by  Bryden  (1976).  Equivalently, 
w  »  -<Tt  +  GTX  +  °Ty)/ez 

If  a  balance  between  local  time  changes  and  horizontal  advection  is 

not  achieved  at  the  mooring  site,  then  vertical  velocities  must  exist. 

That  is  in  fact  the  case  for  all  depths  on  the  mooring.  At  575  dbar, 

uT»  has  generally  small  amplitude  and  is  very  noisy;  also,  Tt  and 

vT^  do  not  show  any  significant  negative  correlation,  which  would  be 

the  case  if  T  ~  -vT».  Although  conventional  scaling  arguments 

^  y 

(e.g.,  Pedlosky,  1979)  suggest  that  uT*  ~  vT',  the  rms  value  of  the 

a  y 

latter  at  575  dbar  is  about  three  times  as  great  as  that  of  the  former, 
and  the  rms  value  of  aT/at  is  about  2  1/2  times  that  of  uT*.  These 
results  suggested  calculating  w  at  575  dbar  from  w  =  -(T^  + 
vT-)/oz  while  at  875,  1175,  and  2000  dbar  (2-17)  was  used. 

To  avoid  introducing  yet  more  noise  by  calculating  «z  from  the 
recorded  data,  the  analytic  fits  used  in  the  temperature  correction 


schemes  at  575,  875,  and  1175  dbar  were  used.  At  2000  dbar,  a  constant 
value  of  .0971  °C./l00m  (determined  from  historical  data)  was  used.  At 
575  dbar,  aT/at  was  calculated  from  the  corrected  temperature  time  series 
at  that  level,  v  was  also  the  corrected  value,  and  au/az  (used  to  obtain 
Tp)  was  determined  from  the  recorded  velocity  and  pressure  time  series 
at  400  and  700  m.  Similar  schemes  were  used  at  the  other  depths: 
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The  latter  forms  of  0Z  are  used  when  data  from  the  1000  m  instrument  is 
lacking,  and  they  are  based  on  a  comparison  of  the  time  series  of  rfz 
calculated  in  the  two  different  ways.  In  the  calculations  for  2000  m, 
^3000  1S  the  direction  the  velocity  linearly  interpolated  between 
2000  and  4000  m.  At  4000  m,  there  is  a  lack  of  vertical  information  to 
compute  w  using  the  temperature  equation;  however,  since  this  current 
meter  was  only  688  m  off  the  bottom,  a  fair  idea  of  the  deep  vertical 
velocities  may  be  obtained  by  calculating  w^qqq  =  u.tfh,  where  Vh  is  the 
slope  of  the  bottom  topography.  At  the  time  of  mooring  recovery,  a 
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Table  2-6,  Top:  Record-lenqth  mean  and  rms  velocities  for 
vertical  velocities  calculated  as  described  in  text.  Value 
at  1175  db  is  based  on  240  days  of  data,  the  rest  on  360  days. 
Bottom:  Cross-products  of  different  vertical  velocity  pairs, 
and  their  correlation  co-ef f icients , 

bathymetric  survey  was  conducted  about  the  mooring  site;  the  survey 
indicated  a  slope  of  about  4  x  10"  ,  with  alongslope  directed  along  70° 
true.  It  should  be  emphasized  that  the  calculation  of  deep  vertical 
velocity  is  absolutely  independent  of  all  those  calculations  made  higher 
in  the  water  column. 

The  resulting  time-series  of  w  at  the  five  depths  are  displayed  in 
Fig.  2.7,  after  being  smoothed  with  a  five-day  running  mean  (except  at 
4000  n,  where  the  time-series  was  already  fairly  smooth).  The  vertical 
velocities  exhibit  good  visual  coherence  throughout  the  water  column,  and 
calculating  correlation  coefficients  for  pairs  of  series  bears  out  what 
is  evident  from  the  plots,  that  vertical  velocities  are  indeed  well 
correlated  throughout  the  water  column.  Table  2-6  lists  the  correlation 
coefficients,  all  of  which  are  significant  at  the  95  percent  confidence 
level  for  an  estimated  30  degrees  of  freedom.  Table  2-6  also  shows  the 
mean  and  rms  values  of  w  at  the  five  levels.  The  mean  values  of  the 
vertical  velocities  are  very  small,  an  order  of  magnitude  smaller  than 
the  rms  values,  and  probably  not  significantly  different  from  zero.  The 
maximum  rms  amplitude  is  attained  mid-depth,  at  875  db;  a  close  look  at 
the  time  series  reveal  that  instantaneously  wmax  occurs  at  thermocline 
levels,  which  vary  from  875  dbar  on  the  warm  side  of  the  Stream  to  575 
dbar  on  the  cold  side.  The  monotonic  decay  above  and  below  the 
thermocline  maximum  appears  in  the  first  EOF  for  the  vertical  structure 
of  vertical  velocity,  which  contains  81.9  percent  of  the  variance  and 
looks  like  the  vertical  velocity  structure  associated  with  the  first 
baroclinic  mode  (Fig.  2.8).  This  result,  taken  with  the  neat  decoupling 


Figure  2.7.  Vertical  velocities  at  designated  levels,  calculated  as 

described  in  text.  From  top  to  bottom  are  vertical  velocities 
at  575,  875,  1175,  2000,  and  4000  dbar.  Units  are  10'3cm/s. 
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Figure  2.8.  Vertical  structure  of  the  first  EOF  calculated  for  vertical 


velocity  from  complete  (360-day)  records  at  4  depths.  The 
dashed  line  goes  to  w  =  0  at  the  surface,  which  Is  the  boundary 
condition  there. 


Eigenfunction  Amplitude 


of  horizontal  velocities  into  just  two  modes,  suggests  that  a  two-layer 
model  (which  reproduces  a  barotropic  and  baroclinic  mode)  can  adequately 
describe  the  data. 

2.5.2  Vorticity  balances  at  the  mooring  site 

A  discussion  of  the  vorticity  equation  is  limited  by  the  lack  of 
horizontal  resolution  afforded  by  a  single  current  meter,  due  to  the 
multiple  derivatives  involved,  for  example,  in  the  advective  terms. 
However,  some  general  features  should  be  recognized  before  trying  to 
determine  the  primary  balances:  l)the  bottom  slope,  with  a  value  of 
about  4  x  10  ,  produces  a  B-effect  of  f  |Vh(/H  =  (0.89  x  10”*  s"1) 
x(4  x  l0"3)/4000  m  =  9  x  10~ilm”*s~1 ),  4  to  5  times  the  planetary 
B-effect;  2)on  the  anticyclonic  side  of  the  Stream,  the  analytic  form 
for  u  gives  a  maximum  velocity  curvature  at  575  db  of  (u»*)max  = 

Q  11  a 

3.60  x  10”  m”  s”  ,  two  orders  of  magnitude  greater  than  b,  and  b  -  u»*  <0 

from  the  axis  to  some  87  km  south  of  the  axis;  and  3)  calculating  fwz 

from  the  time  series  of  w  at  various  depths  turns  out  to  be  fairly  noisy; 

but  rough  estimates  show  that  it  is  generally  an  order  of  magnitude 

larger  than  bv.  Suppose,  *or  example,  there  is  a  change  in  w  of  50  x 
-3 

10  cm/s  from  875  dbar  to  4000  dbar.  The  v  required  to  balance  fwz 
is  then 

«  .  f  =  (0.8892  x  10-4)(5D  »  W3)  cm/$  .  „  cm/ 

«  3125  x  1.9  x  W11 

which  is  greater  than  the  maximum  value  of  v  in  that  part  of  the  water 
column  by  a  factor  of  about  three.  Moreover,  the  above  is  a  conservative 


estimate  of  f  wz> 


Although  these  two  terms,  the  only  two  in  (2- 14b)  that  are 
straightforward  to  evaluate,  are  consistently  of  different  magnitude, 
they  are  nevertheless  correlated  strongly  enough  to  suggest  that  there  is 
an  important  dynamical  connection  between  them.  There  is  even  a  stronger 
correlation  between  v  and  w  at  any  particular  depth  than  between  v  and 
w2,  a  result  that  may  be  interpreted  in  two  ways.  If  the  system  is 
indeed  behaving  as  the  first  baroclinic  mode  suggested  by  Figure  2.8, 
then  w  and  wz  should  display  the  same  behavior  in  time,  but  wz  is  in 
general  a  much  noisier  term--  hence  the  higher  correlation  between  v  and 
w.  The  second  interpretation  depends  on  water  parcels  conserving  their 
temperature:  to  the  extent  that 

-vTy  Tv 

w  ~  and  ~  constant,  w  and  v  ought  to  be  correlated.  This 

z  z 

notion  is  elaborated  on  in  Chapter  5. 

The  analytic  values  for  terms  like  u»*  have  been  used  to  make 
daily  estimates  of  u^t  and  vu^;  however,  the  results  are  very  noisy, 
and  no  systematic  balances  emerge  from  these  calculations.  It  may  be 
fairly  concluded  that  data  from  the  mooring  alone  is  inadequate  to  make  a 
complete  dynamical  investigation.  The  curvature  of  the  Stream  may  be 
important  to  the  vorticity  balance,  as  suggested  by  the  satellite 
composites  for  the  year. 

2.6  Discussion  of  quasi-geostrophic  approximation 

In  quasi-geostrophic  theory,  the  Rossby  number  e  is  defined  as  U/fL 
(where  U  is  a  velocity  scale  and  L  a  horizontal  length  scale),  a  measure 


of  tiie  strength  of  advection  relative  to  the  Coriolis  terms  in  the 


momentum  equations,  or  of  relative  to  planetary  vorticity.  For  a 
strongly  non-isotropic  flow  like  the  Gulf  Stream,  there  will  be  a 
different  estimate  for  each  momentum  equation:  for  the  long-stream 
direction  e  ~  U/fLu  and  for  the  cross-stream  eu  ~  V/fL.  If  U  »  V  and 
Lx  >>  Ly,  clearly  ex  >>  ey.  As  long  as  ex  <<  1,  however, 
quasi -geos trophic  theory  still  applies,  and  variables  can  be  expanded  in 
a  power  series  of  ex>  In  what  follows,  all  estimates  are  implicitly  of 
ex,  since  what  is  important  is  the  upper  bound  on  the  Rossby  number. 

The  subscript  x  is  dropped  for  convenience.  The  Rossby  number  crops  up 
repeatedly  as  a  limit  on  the  relative  size  of  other  flow  attributes,  such 
as  vertical  velocities  or  Isotherm  slopes.  From  the  data  and  analytic 
fits,  it  is  possible  to  estimate  the  Rossby  number  from  these  other 
quantities  as  well  as  from  the  original  definition.  When  the  mooring  is 
in  the  Gulf  Stream,  velocities  range  well  over  100  cm/s  at  575  dbar  (see 
Fig.  2.1c);  the  scale  of  cross-stream  variation  is  about  50  km,  from  the 
fits  (2-7a,b).  Thus,  the  definition  e  =  U/fLy  gives 
100  cm/s  oo 

e  .  =  .ci, 

.89  x  10  s  x  50  km 

It  was  pointed  out  in  Section  2.4,2  that: 
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Using  the  analytic  forms  only,  it  was  found  that  (Table  2-4) 
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where  the  values  of  e  fall  off  rapidly  away  from  the  Stream  axis.  These 
estimates  are  probably  unrealistically  large.  One  e-folding  distance 
south  of  the  maximum  in  -u*,  (2-9c)  gives  u*/f  ~  .24,  more  in  line 

with  the  above  estimate.  The  average  velocity  section  constructed  from 
the  data  may  be  used  to  estimate  a  maximum  ratio  for  -u»/f  of  0.31. 

When  the  vertical  velocity  calculations  were  presented,  no  mention 
was  made  of  their  absolute  magnitude;  but  anyone  familiar  with  estimates 
of  w  for  other  flow  regimes  might  be  uncomfortable  with  rms  magnitudes  of 
0.05  cm/s  (although  Johns  and  Watts  (1985)  recently  estimated  similar 
values  of  vertical  velocity  in  the  Gulf  Stream  100-200  km  northeast  of 
Cape  ’iatteras).  Consider  a  scale  analysis  of  the  continuity  equation. 
With  (U,L  ) ,  ( V,L  )  denoting  long-  and  cross-stream  (velocity, 
length)  scales,  and  (W,H)  the  corresponding  vertical  scales,  then 
according  to  quasigeostrophic  theory. 
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so  that  the  horizontal  velocity  can  remain  nondivergent  to  lowest  order 

in  e.  For  the  Gulf  Stream,  typical  scales  are  U  =  100  cm/s,  V  =  5  cm/s, 

_2 

l  =  1000  km,  L  =  50  km,  H  =  1  km,  and  w  =  5  x  10  cm/s.  Then  e  ~ 
x  y 

WLVUH  -  WL  /VH  =  (0.05  cm/s)(100  km)/(10cm/s)(l  km)  =  0.5.  This 
*  y 

estimate  is  essentially  an  indication  of  the  strong  isotherm  slopes  in 
the  Gulf  Stream.  Quasi-geostrophlc  scaling  of  the  temperature  equation, 
assuming  local  time  derivatives  do  not  dominate,  gives: 
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ay 
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In  these  relations,  H  i 


meant  to  be  the  vertical  scale  of  variation,  e.g.,  for  the  horizontal 
velocities.  Quasi -geostrophy  requires  that  an  isotherm  not  vary  Its 
depth  an  0(1)  amount  within  the  length  scales  under  consideration.  How¬ 
ever,  an  inspection  of  historical  hydrographic  sections  shows  that 


while  with  H  =  1  km,  and  Ly  =  50  km,  e  ~  0.5,  just  the  same  estimate  as 

that  associated  with  vertical  velocities. 

Most  of  the  above  Rossby  number  estimates  apply  only  to  the  flow  at 
575  dbar,  that  is,  the  thermocllne  Gulf  Stream,  and  exclusive  of  the 
estimates  depending  on  the  strong  Isotherm  slopes,  they  do  not  exceed 
values  of  .3.  Deeper  in  the  water  column,  similar  estimates  yield 
substantially  lower  Rossby  numbers  --  all,  that  Is,  except  for  those 
derived  from  the  vertical  velocities,  which  are  substantial  throughout 
the  water  column  (see  Table  2-6).  Below  the  575  dbar  level,  thermal  wind 
balance  holds  to  within  10  percent  In  both  the  along-  and  cross-stream 
directions.  At  575  dbar,  using  the  geostrophlc  approximation  probably 
introduces  errors  of  about  20-30  percent,  and  certainly  the  errors  would 
increase  at  still  shallower  levels.  Moreover,  a  Rossby  number  of  0.3 
suggests  that  a  quasl-geostrophlc  dynamical  model  may  be  Inadequate  for 
explaining  the  time-evolution  of  the  Gulf  Stream. 


Chapter  3.  A  diagnostic  investigation  of  energy  budgets  in  a  numerical 
model 

3.1  Introduction 

An  obvious  feature  of  the  Gulf  Stream  as  it  flows  eastward  into  the 
North  Atlantic  is  its  eventual  weakening:  mass  transport  decreases,  and 
mean  and  eddy  energies  fall  off  from  their  strong  maxima  farther  to  the 
west.  The  mechanisms  behind  this  weakening  are  not  thoroughly  understood, 
but  radiating  instabilities,  topographic  features,  and  conversions  of  mean 
kinetic  to  potential  energy  are  all  likely  candidates.  Fofonoff  and  Hall 
(1983)  have  documented  downstream  changes  in  mass,  momentum,  and  kinetic 
energy  fluxes  as  evaluated  from  hydrographic  data  spanning  sixteen  degrees 
of  longitude,  or  roughly  1500  km.  While  mass  transport  decreases  to  the 
east  only  mildly  over  the  survey,  kinetic  energy  flux  falls  from  its 
maximum  at  the  westernmost  section  to  ten  percent  of  that  value  at  the 
easternmost.  Dewar  and  Bane  (1984)  have  estimated  a  corresponding 
acceleration  of  the  Stream  in  the  South  Atlantic  Bight.  Such  dramatic 
change  in  the  overall  structure  of  this  strong  current  has  implications  for 
the  energetics  throughout  the  subtropical  gyre. 

Different  approaches  have  been  taken  in  examining  this  problem. 

Talley  (1982)  discussed  the  possible  radiation  of  energy  by  barotropic, 
baroclinic  and  mixed  instabilities  for  a  number  of  configurations  (meant  to 
approximate  a  Gulf  Stream-like  current)  in  a  linear,  two  layer,  unbounded 
model.  A  number  of  observational  studies  (Hansen,  1970;  Halliwell  and 
Mooers,  1983;  Watts  and  Johns,  1983)  have  sought  to  document  the 
predominant  time  and  space  scales  of  fluctuations  in  various  parts  of  the 
Gulf  Stream,  but  such  studies  do  not  establish  the  dynamic  reason  for  the 
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meant  to  be  the  vertical  scale  of  variation,  e.g.,  for  the  horizontal 
velocities.  Quasi -geostrophy  requires  that  an  isotherm  not  vary  its 
depth  an  0(1)  amount  within  the  length  scales  under  consideration.  How¬ 
ever,  an  inspection  of  historical  hydrographic  sections  shows  that 
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while  with  M  =  1  km,  and  =  50  km,  e  ~  0.5,  just  the  same  estimate  as 

that  associated  with  vertical  velocities. 

Most  of  the  above  Rossby  number  estimates  apply  only  to  the  flow  at 
575  dbar,  that  is,  the  thermocline  Gulf  Stream,  and  exclusive  of  the 
estimates  depending  on  the  strong  isotherm  slopes,  they  do  not  exceed 
values  of  .3.  Deeper  in  the  water  column,  similar  estimates  yield 
substantially  lower  Rossby  numbers  --  all,  that  is,  except  for  those 
derived  from  the  vertical  velocities,  which  are  substantial  throughout 
the  water  column  (see  Table  2-6).  Below  the  575  dbar  level,  thermal  wind 
balance  holds  to  within  10  percent  in  both  the  along-  and  cross- stream 
directions.  At  575  dbar,  using  the  geostrophic  approximation  probably 
introduces  errors  of  about  20-30  percent,  and  certainly  the  errors  would 
increase  at  still  shallower  levels.  Moreover,  a  Rossby  number  of  0.3 
suggests  that  a  quasi-geostrophic  dynamical  model  may  be  inadequate  for 
explaining  the  time-evolution  of  the  Gulf  Stream. 
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Chapter  3.  A  diagnostic  Investigation  of  energy  budgets  in  a  numerical 
model 

3.1  Introduction 

An  obvious  feature  of  the  Gulf  Stream  as  it  flows  eastward  into  the 
North  Atlantic  is  its  eventual  weakening:  mass  transport  decreases,  and 
mean  and  eddy  energies  fall  off  from  their  strong  maxima  farther  to  the 
west.  The  mechanisms  behind  this  weakening  are  not  thoroughly  understood, 
but  radiating  instabilities,  topographic  features,  and  conversions  of  mean 
kinetic  to  potential  energy  are  all  likely  candidates.  Fofonoff  and  Hall 
(1983)  have  documented  downstream  changes  in  mass,  momentum,  and  kinetic 
energy  fluxes  as  evaluated  from  hydrographic  data  spanning  sixteen  degrees 
of  longitude,  or  roughly  1500  km.  While  mass  transport  decreases  to  the 
east  only  mildly  over  the  survey,  kinetic  energy  flux  falls  from  its 
maximum  at  the  westernmost  section  to  ten  percent  of  that  value  at  the 
easternmost.  Dewar  and  Bane  (1984)  have  estimated  a  corresponding 
acceleration  of  the  Stream  in  the  South  Atlantic  Bight.  Such  dramatic 
change  in  the  overall  structure  of  this  strong  current  has  Implications  for 
the  energetics  throughout  the  subtropical  gyre. 

Different  approaches  have  been  taken  in  examining  this  problem. 

Talley  (1982)  discussed  the  possible  radiation  of  energy  by  barotropic, 
baroclinic  and  mixed  Instabilities  for  a  number  of  configurations  (meant  to 
approximate  a  Gulf  Stream-like  current)  in  a  linear,  two  layer,  unbounded 
model.  A  number  of  observational  studies  (Hansen,  1970;  Halliwell  and 
Mooers,  1983;  Watts  and  Johns,  1983)  have  sought  to  document  the 
predominant  time  and  space  scales  of  fluctuations  in  various  parts  of  the 
Gulf  Stream,  but  such  studies  do  not  establish  the  dynamic  reason  for  the 


existence  of  the  fluctuations.  Bryden  (1979)  and  Wright  (1981)  have 
addressed  analogous  questions  in  the  Antarctic  Circumpolar  Current  by 
modelling  the  flow  as  strictly  zonal  and  confined  to  a  zonal  channel,  then 
investigating  the  normal  mode  linear  instabilities  that  can  arise. 

Haidvogel  and  Holland  (1978)  have  used  a  similar  model  on  both 
instantaneous  and  mean  profiles  of  the  jet  from  a  two-layer 
quasi -geos trophic  numerical  model,  and  related  the  results  to  observed 
scales  of  variability  in  the  model.  Johns  (1985)  recently  used  a  linear 
baroclinic  instability  model  successfully  to  predict  time  and  space  scales 
of  fluctuations  100-200  km  northeast  of  Cape  Hatteras,  where  linear  theory 
might  still  be  expected  to  apply.  Pedlosky  (1970,  1982)  has  investigated 
the  effects  of  allowing  the  instabilities  to  attain  finite  (but  still 
small)  amplitude. 

There  are  philosophical  as  well  as  technical  problems  in  extending 
such  investigations  to  the  Gulf  Stream  where  it  becomes  a  free 
eastward- flowing  current.  Observations  have  long  suggested  that  the  Gulf 
Stream  is  strongly  non-linear,  so  that  Pedlosky's  weakly  nonlinear  theory 
has  dubious  application.  Furthermore,  while  channel  models  may  be 
geometrically  realistic  for  the  Antarctic  Circumpolar  Current,  the  Gulf 
Stream  is  certainly  not  confined  to  a  channel  but  is  free  to  meander  north 
and  south.  While  observations  serve  to  describe  the  fluctuations,  they  can 
do  little  to  explain  them:  because  a  wave  radiates  energy  primarily 
through  correlations  of  velocity  and  pressure,  one  would  like  to  measure 
these  properties  to  evaluate  the  radiation.  However,  for  the  energy  budget 
of  an  enclosed  area,  it  is  the  divergence  of  the  horizontal  pressure  work 


that  contributes,  and  in  a  prominently  geostrophic  regime,  this  quantity  is 
degenerate  at  lowest  order.  Measurement  techniques  are  not  yet  refined 
enough  to  give  reliable  estimates  of  correlations  between  pressure  and 
small,  ageostrophic  velocities. 

On  the  other  hand,  numerical  eddy- resolving  general  circulation 
models  (EGCM's)  can  provide  not  only  the  degree  of  non-linearity  actually 
observed  in  the  ocean,  but  also  time-series  of  "data"  as  extensive  and 
detailed  in  space  and  time  as  desired.  W.  R.  Holland  of  NCAR  has  created 
several  "generations"  of  quasi-geostrophic  EGCM's,  and  has  recently 
reviewed  (Holland,  1984)  the  successes  and  compromises  involved  in 
comparing  such  models  to  the  real  ocean.  In  all  the  models,  there  is  a 
strong  western  boundary  current  that  turns  eastward  at  some  latitude  to 
become  either  a  northern  boundary  current  (in  some  single-gyre  models)  or  a 
free  jet  flowing  into  the  interior  of  the  basin.  In  a  number  of  early 
models,  the  free  jet  did  not  penetrate  realistically  far  into  the  model 
basin  before  dying  out  completely.  Holland  and  Schmitz  (1984)  discuss  the 
factors  controlling  the  penetration  scale,  based  on  a  study  of  many  two-, 
three-,  and  eight-layer  models.  They  find  that  the  strength  of  friction 
and  topography,  as  well  as  the  vertical  resolution  of  the  model,  affect  the 
penetration  scale. 

The  model  chosen  for  analysis  here  was  recommended  by  Holland 
(personal  communication)  as  being  the  most  realistic  model  for  such  an 
investigation,  and  is  referred  to  by  Holland  as  3L-4;  its  features  are 
summarized  in  Table  3-1.  Figure  3.1  shows  the  time-mean  streamfunction  in 
each  layer.  The  free  jet  penetrates  about  halfway  across  the  4000  km 


Layer  depths: 


300  m 


Interfaces: 

Rotation: 

Wind  stress 

Friction: 

Horizontal 


H2  =  700  m 
H3  =  4000  m 

9  3/2  =  m  s 

g ' 5/ 2  =  .0162  m  s“2 

f  =  9.3  x  10  J  s_1 
6  =  2  x  10"11  m' V1 

:  ZQ  =  0.1  N/m2 

-ZQ  cos  (  40^0  |c"m  ^  ®  i.  y  5.  ^000  k  m 

9  -4-2 

lateral:  biharmonic,  A  =  8  x  10  m  s 
bottom:  linear,  y  *  1  x  10“7  $"* 

dimensions:  4000  km  x  4000  km 


Table  3-1.  Summary  of  features  of  Holland  numerical  model  3L-4 


Figure  3.1.  Mean  streamf unction  for  Holland  numerical  model  3L-4 : 
allayer  1;  contour  from  -1.5xl05  to  1.5xl05  m2/s;  Cl  =  104  m2/s;  b)layer 
2:  contour  from  -5.4xl04  to  5.4xl04  m2/s;  Cl  =  6xl03  m2/s;  c)layer  3: 
contour  from  -1.8xl04  to  1.8xl04  m2/s;  Cl  =  2xl03  m2/s. 


square  basin.  Three  layers  appear  to  be  sufficient  for  capturing  the 
proper  degree  of  non-linearity  in  the  free  jet;  moreover,  in  a  three-layer 
as  opposed  to  the  two-layer  models,  there  is  a  layer  shielded  from  both 
surface  and  bottom  effects,  which  displays  vast  regions  of  homogenized 
potential  vorticity.  Since  such  regions  are  predicted  by  theory  (Rhines 
and  Young,  1982)  and  found  observationally  (McDowell  et^h ,  1982),  it 
seems  important  to  work  with  a  model  reproducing  this  feature.  On  the 
other  hand,  adding  still  more  layers  appears  to  do  little  more  than  enrich 
the  vertical  resolution  of  structures  already  observed  In  three-layer 
model s. 

In  the  remainder  of  the  chapter,  first  the  equations  for  mean  and 
eddy  kinetic  energy  are  derived  for  a  three-layer  quasi -geostrophic  system, 
starting  either  with  the  momentum  equations  or  the  potential  vorticity 
equation.  Next  the  calculations  from  the  numerical  data  are  described 
briefly,  finally,  the  results  are  presented,  in  which  essentially  six 
different  energetic  regions  are  identified  and  described. 

3.2  Equations 

In  formally  deriving  the  equations  governing  a  quasi-geostrophic  flow 
regime,  it  is  useful  to  scale  the  dimensions  out  of  most  of  the  variables; 
then  one  or  more  nondimensional  parameters  govern  the  relative  importance 
of  terms  in  any  equation,  and  consistency  Is  achieved  without  confusion. 
However,  when  evaluating  terms  of  an  equation  from  a  data  set,  it  is  less 
awkward  and  physically  more  intuitive  to  work  with  a  dimensional  set  of 
equations,  applying  the  results  of  quasi-geostrophic  theory  directly. 
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Thus,  in  this  text,  such  results  will  be  cited  frequently  without  rigorous 
derivation;  anyone  unfamiliar  with  quasi -geostrophic  theory  is  referred  to 
Pedlosky  (1979).  Notation  is  kept  as  much  as  possible  like  Holland's 
(1984)  notation,  but  in  addition,  expansion  of  variables  in  terms  of  a 
small  parameter  --  usually  the  Rossby  number  e  --  will  be  allowed  for. 

Then  subscripts  refer  to  layers,  superscripts  to  the  order  of  a  variable  in 
its  Rossby  number  expansion.  Table  3-2  gives  other  notation  and 
conventions. 

In  a  three-layer  system,  there  are  five  variables  to  solve  for,  u, 
v,  w,  p,  and  h;  the  governing  equations  are  as  follows: 

a)  In  each  layer,  there  are  two  horizontal  momentum  equations  and  a 
continuity  equation: 

,u-  ,uk  *  v„  3v  ...  -1  . 


‘<7 £*" 


k  —  k  rr]  -f\  *  ’^3x  8(k'1)  fyk-s<k-3>”“k  -  A-’4  “k 
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k  =  1,  2,  3 


(3-lc 


b)  At  each  interface,  a  "thermodynamic’1  equation  for  evaluating  w, 
and  a  relation  between  interface  height  and  pressure: 


ah 


ah. 
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k+S 
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EKE 

MKE 

RHS 

LHS 


(east,  north)  coordinates 
(east,  north)  velocity  in  layer  k 

deviation  of  interface  between  layers  k,  k+  1  from  rest 
state;  h  >  0  for  upward  displacement 

(const amt)  thickness  of  layer  k  in  rest  state 

vertical  velocity  of  interface  h^^ 

basic  density  of  fluid 

actual  density  of  fluid  in  layer  k 


Pk+1  “  Pk 


deviation  of  pressure  from  basic  state,  where  basic  state 
is  in  hydrostatic  balance 


k+*i 


g  *  9.81  m/s 


linear  (bottom)  friction  coefficient 
biharmonic  (lateral)  friction  coefficient 


KE^  =  kinetic  energy  in  layer  k 

*  +  ±1' 

>s  H,  (u,/2  +  vT^) 


VHk(uk2  +  vk2) 


available  potential  energy 
eddy  kinetic  energy 
mean  kinetic  energy 
right  hand  side 
left  hand  side 


Table  3-2.  Definitions  for  noatation  used  in  Chapter  3. 
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Boundary  conditions  at  top  and  bottom  are  that  w  =  0  there.  The 
Kronecker  delta  function  Is  used  to  Include  forcing  and  dissipation 
occurring  only  In  particular  layers.  Wind  stress  and  bottom  friction  act 
as  body  forces  on  a  layer.  Each  variable  will  be  expanded  In  terms  of  the 
Rossby  number  e,  which  must  be  <<  1  for  quasi -geostrophlc  theory  to  be 
valid.*  Then,  for  example, 

uk  =  u^°^+U||^e+u^^e^+. . .  Rarely  will  more  than  two  terms  of 

any  expansion  be  required  here. 

In  a  quasi-geostrophlc  regime,  the  flow  field  is  horizontally 
non-divergent  at  lowest  order  so  that  a  streamfunctlon  ^  can  be 


introduced  --  in  fact,  rk  =  — p-,  and 


•3fk  3*k  .  f0  ,4,  iL  S 

ay"’  k  =  3x“*  k+1/2  3  gj^J  lVl  '  rk> 


(3-2) 


Inserting  the  e-expanslon  for  each  variable,  cross-differentiating 
(3-la,b),  substituting  from  (3-ld)  and  (3-2),  and  retaining  only  lowest 
order  quantities,  one  can  derive  a  vortlclty  equation: 


(3-3) 


Tt^2\  +  J(Vf  *  =  (wk-l/2  "  Vl/2) 

-  8(k-D  1T17  -  «(k-3)  _  AV6t^ 

which  may  be  further  rearranged  using  (3-1)  and  (3-2)  to  obtain  the 

potential  vortlclty  equation  as  given  in  Holland  and  Schmitz  (1984): 


It  Is  also  required  that  forcing  and  dissipation  not  enter  the  lowest 
order  momentum  balances  (3-la,b). 


.  -  -  *  »  ■  v  ■  "  •»_ 
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dQ,  ov 

Ht“  =  VFk*  Qk=7  k+fo+By+fl7(hk+l/2"hk-l/2) 


(3-4) 


d  3 


=  +  ).  Here,  of  course,  Tk  = 


s{k-1)  Twrh"' 


while  Fk  =  -<5 (k-3)  yV2*k  -  AV6^. 


Now  consider  the  derivation  of  energy  equations.  Of  primary  interest 
here  is  the  kinetic  energy,  which  to  lowest  order  is  just 

i  u  fe  <-k(0’2  *  V0,2> 

KEk  ■  ?  «k  Kr-"^ — rJs — 

since  at  this  order  the  vertical  velocity  makes  no  contribution.  Although 
in  the  end  the  equations  are  integrated  over  depth,  separate  equations  for 
each  layer  are  derived  first  to  show  how  energy  is  exchanged  vertically 
between  layers.  An  energy  equation  for  each  layer  is  obtained  by 
multiplying  (3-4)  by  -^k  Hk,  then  rearranging  terms  to  obtain: 

hr-  (*ix+  ^iy,‘fo^i  ir^  _Hi  fx^^iV^Vtq’  h3/2)+*i^ixt} 

-H1  fy^lVV  >^h 3/^Vlyt)  =  *lfA  (3’5a) 
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+  ^2xt)  -  »2  |y{v2^(V2V  ^h5/rh3/2))^2f2yt]-.  -AH 2f/f2  (3-5b) 
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"H3  h5/2)  +  ^3yt}  =  ~  M3^6% 


In  Appendix  A.l  is  derived  the  useful  identity 
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ay  1  'k"k "H'k '  ax  '“k  Hk  '  ay 
In  layer  k,  the  total  derivative  following  the  fluid  is 


(3-6) 


it  =  at  +  uk  ax  +  vk  ay  =  1ST  +  J^k,  * 
while  at  an  interface  k  +  1/2  it  is 

dk+l/2  a  .  a  .  „  a  3  +  » 

it  a  TE  uk+l/2  ax  vk+l/2  ay  *  1TE  Vl/2*  ' 


where  the  streamf unction  at  the  interface  is  defined  following  Holland* s 

i/  ^k  ^k*l+  ^k+1  ^ k 

(1978)  convention  as  ^k+1/2  =  - q — TTl - -•  In  ^act,  as  shown  in 


Appendix  A. 2, 

J(Vl/2,hkn/2)  =  J(^k,hkn/2)  3  J(^k+l,hk*l/2) 

so  that  In  future  the  subscript  on  the  total  derivative  Is  omitted,  since 
it  Is  evident  what  is  meant  by 

Using  the  above,  equations  (3-5)  can  now  be  written  as 


( 3-  7  a ) 
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dhc 


dv. 


u  d  rl,*2  .  ^  2n  f  U  ■sB/2  +  u  3/L-^  +  f  A 
H3  Ht  Vr3x  V3  =  -fo  3  dt  H3  17^?® 3  3  ’3T) 


du. 


*H3  ly  <‘*3  TT>  *  *3  A  *  ^sA  l3-7c) 

Finally,  it  is  useful  to  break  down  the  first  term  on  the  RHS  of  (3-7a)  as 


„  dhi/?  ,  f_H  •  . 

fo\  tfE  =  fo^lw3/2  =  fo  ~ TT^fCT  ”  *3/2  *  HT^TT  w: 
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H1 

=  fo^3/2W3/2  "  ffJ+Hj  h3/2  W3/2  g'3/2 


(3-8) 


Similar  treatment  for  the  other  two  layers  gives  the  final  versions  of  these 
energy  equations: 

Hi  ,i  jo  dv,  .  ,  du. 


It  KE1  =  fo^3/2W3/2  ’  g,3/2h3/2w3/2+Hl  IT^lly* ^HT1 
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+H2  htbfyh  ^r-)+H 2  ly^'^2  Tr)+AH2V6,2 


(3-9b) 


^3  =  -f o^5/2 w5/2~ H^TT9 '5/2h5/2W5/ 2+H3  !x<^3+  ?  7T)+H3  ly*'^ 


^3V2V3  +  M&6% 


(3-9c) 


The  first  term  on  the  RHS  of  each  equation  Is  like  a  vertical  pressure  work 
term,  which  transfers  energy  between  the  layers.  The  second  term  is  the 
conversion  In  that  layer  of  available  potential  to  kinetic  energy,  where 
the  total  available  potential  energy  Is 

D  1  fo  x  2J  *o  Mk-V',)2  1„.  u2  A  1„,  „2 

P  =  7  <TyY{2  '  V  7  cTTJ*  3  2  =  793/2h3/2+  Z*5/2h5/2  • 
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The  total  conversion,  summed  over  the  three  layers,  Is 

•g,3/2W3/2h3/2  -  g,5/2w5/2h5/2*  To  obta1n  *i%  result 
directly,  multiply  (3- Id)  by  9'|l+i/2hk+i/2  and  sum  over  tbe  two 
Interfaces.  The  result  is 

■ffE  ■  g*3/2w3/2h3/2  +  g5/2w5/2h5/2* 

The  Interpretation  of  the  other  terms  on  the  RHS  of  (3-9a-c)  remains 
ambiguous,  so  consider  now  the  horizontal  momentum  equations  (3-la,b).  To 
form  a  kinetic  energy  equation,  multiply  (3-la)  byfH^  +  ~  hk+l/2^uk» 

(3-lb)  by  (Hk  +  \_y2  -  hk+i/2^  vk*  and  add*  Note  that  v*hereas  before 
It  was  understood  that  u^,  v^  were  the  lowest  order  portions  In  the 
e-expanslons,  now  one  must  be  careful  to  Include  enough  orders  in  each 
variable's  expansion  so  as  not  to  miss  any  part  of  the  balance  (hence  the 
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retention  of  h^^  the  layer  depths).  Omitting  the  forcing  terms 
in  layer  1,  which  are  straightforward  to  Include,  noting  that 

U  3?  =  and  writ^n9  h3/2  =  eHln3/2  *n3/2*s 


and  non-dimensional)  yields  for  that  layer: 
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The  lowest  order  balance  is  degenerate: 

pi0)  p[°> 

0  ■  'Hlu[0)  lx  "o 'Hri0)  ly  “o 


since  u[°^pj°^x  +  vi°^Pi°^y  =  °*  The  next  order  yields  the  equation 


for  the  lowest  order  kinetic  energy: 
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AH1(u[o)V4u|o)+v}o)V4v{o)) 

when  forcing  is  included.  The  first  term  on  the  WiS  can  be  written 
(p{0,+ep{1,){u{0)*eu{1))  +7T(Pl0)+ePl1),,-(ul0,+  eUl11)  = 


H  w(1)  u 

n  1  /  r\  \  /r>  -H  i 


poe  po  1  "l  poe  o  1  3/2 

Similar  treatment  of  all  three  layers  yields: 

h  y, 
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d  ,  -^3 

dt(KE3)a"V3W5/2  V^.(P3^3)-AH3(f3y74lf3y^3x74f3x)-^3(^+f2)  (3.12c)< 

The  first  term  on  the  RHS  of  each  equation  is  the  same  one  that  was 
rewritten  in  (3-9)  as  the  sum  of  a  vertical  transfer  of  energy  and  a 
conversion  between  potential  and  kinetic  energies.  The  second  term  is  the 
divergence  of  the  horizontal  pressure  work,  a  term  that  ultimately  will  be 
calculated  directly.  The  remaining  terms  are  input  and  extraction  of 
energy  by  wind  and  lateral  and  bottom  friction. 

The  next  step  is  to  note  the  following  identities: 


Wlx^x1  *ly 

V\  -  &*Ax>  *  -  V\x-  %'% 

If  these  are  put  into  equations  (3-9)  then  a  comparison  of  terms  between 
the  resultant  forms  and  equations  (3-12)  shows,  finally: 

7^  ^2e  V  \  ’ar+A W+Hi  V4V 


^.(p2&2)=h2  \  ^r+AV4^2x)+H2  t^-%. 


f*(p3&3)=H3  IT+AV4^3x+t^3x)+H3 


(3-13 


Thus,  the  divergence  of  horizontal  pressure  work  is  calculable  entirely  in 
terms  accessible  from  the  numerical  model  data. 

It  is  clear  that  in  the  mean  the  Gulf  Stream  weakens  as  it  flows 


eastward  into  the  North  Atlantic  (Fofonoff  and  Hall,  1983);  the  same  is 
true  of  its  analog,  the  free  jet,  in  the  numerical  model.  Yet  the 
possibility  must  be  anticipated  that  variations  from  the  mean  are  important 
in  determining  the  mean  distribution  of  energy  in  the  ocean  basin.  Thus, 
all  quantities  are  separated  into  time-averaged  and  time-varying  components 


( 
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and  T  is  the  total  length  of  the  time  series  being  analyzed.  Thus,  by 
definition  =  0.  Time-varying  or  "eddy"  terms  are  not  restricted  in 

magnitude  relative  to  the  mean.  It  is  evident  from  Fig.  3.1  that  because 
the  model  forcing  and  equations  are  symmetric  about  the  mid-latitude  of  the 
model  domain,  the  time-mean  jet  flows  directly  eastward.  Notice  that  this 
energetic  analysis  treats  the  jet  as  the  Eulerian  average  flow  of  a 
particular  region,  not  as  a  feature  changing  position  and  orientation 
constantly.  There  are  two  reasons  for  having  taken  this  approach:  first, 
it  was  not  technically  feasible  to  take  the  latter  approach  due  to 
computing  limitations;  second,  energetic  analyses  have  traditionally  been 
conducted  in  an  Eulerian  frame,  so  for  comparison  it  makes  sense  to  use 
that  approach  here. 

Bryden  (1982)  discusses  the  derivation  of  equations  describing  the 
change  of  kinetic  energy  of  the  mean  flow, 

if?2  +  ^)>  and  of  the  mean  value  of  the  kinetic  energy  of  the  eddies, 

i(Y2  +  f'^) .  Inteyrated  over  all  three  layers,  the  equations  are 
u  x  y 


respectively. 
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( 3- 14b ) . 


Several  comments  are  in  order  here.  First,  superscripts  denoting  the 
e-expansion  order  of  variables  have  been  dropped:  in  all  cases  save  the 
pressure  work  term,  only  lowest  order  quantities  are  required.  Since  the 

calculation  of  -V.(pti)  actually  is  made  from  (3-13),  which  involves  only 
lowest  order  variables,  the  superscripts  are  omitted  for  neatness.  Second, 
it  should  be  noted  that  the  time- varying  field  obeys  quasi-geostrophic 
dynamics  as  well,  whence  u'^  =  vk=^  kx'  As  a  resu^t*  any 

term  of  the  sort  D.V[(  )]  can  be  written  7 . C ti(  )]  regardless  of 

how  time-averaging  enters  the  term.  This  last  point  will  be  useful  in 
interpreting  the  equations  in  terms  of  kinetic  energy  fluxes. 

The  equation  for  the  time  average  of  the  total  kinetic  energy  is 


obtained  either  by  adding  ( 3- 14a  and  b)  or  directly  by  time-averaging  the 
sum  of  equations  (3-12).  The  result  is: 


(3-15) 
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where  K  =  ^-(u^+v^lH,  K'^u  ,2+v|2)H.  In  (3- 14a),  the  mean  advection  of  mean 

kinetic  energy  is  balanced  by  the  pressure  work  term  due  to  the  mean  field, 
conversions  to  and  from  mean  potential  energy,  input  of  energy  by  the  wind, 
removal  of  energy  by  bottom  and  lateral  stresses  acting  on  the  mean  field, 
and  a  term  which  represents  the  conversion  of  eddy  to  mean  kinetic  energy. 

In  ( 3- 14b ) ,  the  mean  plus  eddy  advection  of  eddy  kinetic  energy  is  balanced 
by  the  horizontal  pressure  work  due  to  the  eddy  field,  conversions  to  and 
from  eddy  potential  energy,  losses  to  friction  due  to  the  eddy  field,  and  a 
term  giving  the  conversion  of  mean  to  eddy  kinetic  energy.  When  the  two 
equations  are  added,  as  Bryden  (1982)  has  pointed  out,  the  conversion  terms 
between  mean  and  eddy  kinetic  energy  do  not  cancel,  but  give  rise  to  the 
pecul iar-lookinq  part  of  the  divergence  on  the  LHS  of  (3-15),  i.e., 

I-tCu  u'  +  7  u’V)  +  -f-(u  u'v'  +7  V2) .  The  mystery  of  this  term  dis- 

a  X  Oj 

appears  when  it  is  observed  that  the  LHS  of  (3-15)  is  just  the  total 
divergence  of  the  total  kinetic  energy  flux,  time-averaged: 


K  K 
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LHS  of  (3-15)  >  ^  (|j  okHk  ^  *  |y  «kHk^r!L  >• 


The  kinetic  energy  budget  will  be  considered  over  volumes  of  ocean 
with  open  side  boundaries  at  x  =  (xw,x^-)  to  west  and  east,  y  = 

(y$,yN)  to  south  and  north.  Integrating  (3-15)  over  the  area  yields: 
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The  following  points  are  made  regarding  (3-16): 
|XE 


16) 


1)  (*) 


Is  the  value  of  {*)  at  xc  minus  its  value  at  x..; 


2)  the  divergence  of  pressure  work  has  been  replaced  by  the  terms 
used  to  calculate  it,  but  the  individual  fl uxes  across  the 
boundaries  are  not  necessarily  equal  to  the  fluxes  of  pu: 
only  the  net  divergence  over  the  volume  has  any  meaning  in  this 
context; 

3)  those  terms  for  which  separate  mean  and  eddy  components  have 
been  calculated  are  broken  down  into  those  components  in  (3-16) 

Before  the  kinetic  energy  budgets  for  the  model  are  discussed,  the  model 
data  and  calculations  are  briefly  described. 

3.3  Calculations 

Important  attributes  of  numerical  model  3L-4  analyzed  here  are  given 
in  Table  3-1.  Little  will  be  said  concerning  the  actual  numerics  used  to 
run  the  model,  as  Chow  and  Holland  (1985)  describe  this  in  detail.  The 
wind  stress  in  3L-4  is  steady  and  drives  two  gyres  with  opposite  vorticity 
input.  The  circulation  is  spun  up  from  a  rest  state,  then  run  for  a  number 
of  years  after  reaching  statistical  equilibrium.  Although  the  model  time 
step  is  a  fraction  of  a  day,  instantaneous  values  of  the  streamfunction  are 
stored  only  every  two  days.  The  data  treated  here  use  streamfunction 
values  for  every  fourth  day,  extending  over  about  four  years  of  model  time, 
resulting  in  a  total  of  360  instantaneous  values  of  streamfunction. 

Due  to  limited  resources,  it  was  necessary  to  restrict  the 
calculations  in  several  ways.  First,  only  a  subset  of  the  entire  domain 
was  investigated,  extending  from  820  km  south  of  mid-basin  to  780  km  north, 
and  from  the  western  boundary  to  the  middle  of  the  basin,  2000  km 


eastward.  (It  was  originally  intended  to  go  from  800  km  south  to  800  km 
north  of  mid-basin,  but  the  data  were  accidentally  offset  by  20  km.) 
Equation  (3-16)  was  evaluated  by  first  calculating  all  the  necessary 
time-averages  at  every  grid-point,  then  using  the  resulting  data  sets  to 
examine  budgets  of  various  volumes.  While  it  seemed  useful  to  break  down 
some  terms  into  mean  and  eddy  contributions,  there  would  simply 

“ ,,  dvk 

be  too  many  moments  to  evaluate  in  a  term  like  -g-, 

so  that  only  the  overall  time-average  has  been  calculated.  Also,  initial 
calculations  of  terms  involving  the  biharmonic  friction  proved  very  noisy 
and  very  small:  indeed,  in  models  of  this  type  (see,  .e.g.,  Holland,  1978) 
this  term  contributes  negligibly  to  the  dissipation  of  energy  (but  is 
important  in  dissipating  enstrophy).  Hence  it  is  not  included  In  the 
budgets. 

Evaluation  of  terms  in  the  Instantaneous  vorticity  equation  using  the 
model  data  yields  relatively  large  imbalances.  Evidently  that  is  due  to 
the  poor  estimate  of  local  time  change  of  vorticity  afforded  by  the 
sampling  interval  of  four  days  (Holland,  personal  communication). 

Similarly,  calculations  of  "instantaneous”  kinetic  energy  budgets  show  that 
on  a  short  time  scale  there  are  apparent  imbalances,  but  again, 

|^KE  would  be  poorly  estimated.  On  the  other  hand,  once  the  averages  of 

terms  have  been  formed,  the  well-defined  meridional  structure  of  individual 
quantities  suggests  the  time-averaged  budgets  give  a  reliable  assessment  of 
the  energetics  in  the  domain.  Also,  due  to  the  strength  of  flow  right  at 
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the  western  boundary,  the  instantaneous  values  may  be  extremely  noisy;  thus 
for  examining  the  budgets,  the  western  wall  itself  has  not  been  used,  but 
rather  the  data  from  the  first  grid-point  in  (20  km  east  of  the  boundary). 

First  and  second  derivatives  have  been  calculated  using  ordinary 
centered  finite-difference  schemes.  Any  term  involving  a  Jacobean  has  been 
evaluated  using  an  Arakawan  scheme.  Time-differences  were  also  evaluated 
using  a  centered  scheme,  and  in  terms  where  they  are  part  of  an  average, 

such  as  in  the  poor  estimate  of  local  time  change  may  be  partly  res¬ 

ponsible  for  the  residuals  that  are  found  in  the  energy  budgets. 
Straightforward  trapezoidal  integrations  have  been  used  in  the  horizontal 
in  these  calculations. 

3.4  Results 

Once  the  tools  are  at  hand  to  evaluate  kinetic  energy  budgets  for  any 
volume  of  the  domain,  it  is  important  to  consider  carefully  the  choice  of 
volumes  to  be  investigated.  Sometimes  one  is  interested  in  the  exchanges 
of  energy  between  the  mean  and  eddy  flows;  Harrison  (1979),  for  example, 
investigated  the  interaction  between  the  mean  and  eddy  kinetic  energy 
fields  of  a  variety  of  numerical  experiments  by  integrating  over  the  entire 
basin.  He  points  out  that  local  budgets  of  open-bounded  regions  may  be 
quite  different  from  the  basin-averaged  budgets,  however.  Harrison  and 
Robinson  (1978)  have  in  fact  done  someanalysls  of  this  sort,  which  will  be 
discussed  later.  Here  the  point  is  to  Investigate  the  very  strong 
convergence  or  divergence  of  kinetic  energy  in  the  vicinity  of  a  free  jet 


like  the  Gulf  Stream,  and  to  determine  which,  if  any,  of  the  mechanisms  on 
the  RHS  of  (3-16)  is  dominant  in  producing  a  balance. 

First  a  number  of  volumes  were  examined,  chosen  according  to 
experience  with  the  Gulf  Stream  '60  data.  These  volumes  were  of  small 
enough  size  to  prevent  extensive  overlap  of  energetically  different 
regions,  until  the  gross  energetic  patterns  emerged.  A  comparison  of  the 
areas  suggested  by  this  search  with  the  mean  streamlines  in  all  three 
layers  helped  to  refine  the  horizontal  boundaries  of  the  volumes  and 
offered  rigorous  criteria  for  selecting  them.  Essentially  six  types  of 
energetic  regimes  have  been  identified,  but  since  the  model  is  symmetric 
about  mid-latitude,  four  regimes  appear  on  each  side  of  the  jet,  so  that 
there  are  actually  ten  different  physical  volumes  for  which  budgets  have 
been  calculated.  The  criteria  for  determining  the  latitudinal  and 
meridional  boundaries  of  these  volumes  are  as  follows: 

1)  Of  paramount  importance  is  delimiting  the  north-south  extent  of  the 
jet.  Throughout  the  jet,  the  deep-layer  zonal  velocities  nearly 
always  change  sign  between  the  grid-points  corresponding  to 

y  a  140  km  and  y  =  160  km  north  and  south  of  the  jet  axis  at  y  = 

0.  This  also  turns  out  to  be  the  latitude  at  which  upper  layer  zonal 

p 

velocity  has  dropped  off  to  approximately  1/e  times  its  maximum 
value  for  any  given  section.  Thus,  the  jet's  boundaries  are  chosen 
to  be  y  =  140  km. 

2)  Although  the  motivation  for  this  work  was  the  weakening  of  the  free 
jet,  there  is  a  significant  accelerating  portion  starting  at  the 
western  boundary.  The  jet  is  considered  to  change  from  accelerating 


to  decelerating  where  the  mean  zonal  flux  of  mean  kinetic  energy 
attains  a  maximum,  660  km  from  the  western  boundary. 

3)  From  Fig.  3.1,  it  is  clear  that  in  the  top  and  middle  layers  there  is 
both  a  broad  Sverdrup  circulation  and  a  tighter  Inertial 
recirculation,  while  in  the  deep  layer  only  the  latter  is  visible. 
Holland  and  Rhines  (1980)  have  found  that  the  exact  Sverdrup  balance 
ev  =  fwz  obtains  over  a  relatively  small  area  of  numerical  models, 
but  in  that  work  it  was  found  that  it  did  describe  the  region  in  the 
top  layer  that  appears  in  Fig.  3.1a  as  broadly  spaced  streamlines 
with  a  primarily  meridional  orientation.  The  latitude  at  which  the 
deep- layer  mean  streamf uncti on  goes  to  zero,  defining  the  meridional 
extent  of  the  recirculation,  is  very  close  to  the  latitude  where  the 
mean  vertical  heat  flux  gwp  ,  Integrated  over  the  eastern  part  of  the 
domain,  changes  sign.  The  latitude  midway  between  the  two  is  chosen 
as  the  boundary  for  this  regime.  The  last  region,  referred  to  as  the 
Sverdrup  regime,  extends  from  the  edge  of  the  recirculation  to  the 
open  boundaries  of  the  domain,  at  y  =  -820  km  and  y  =  +780  km. 

4)  For  each  of  the  latitude  bands  defined  In  3),  an  east-west  division 
is  also  made  at  x  =  660  km,  to  correspond  to  the  division  made  in 
the  jet  regions.  As  will  be  seen,  the  eastern  and  western  portions 
of  the  recirculation  and  Sverdrupian  regimes  are  characterized  by 
different  energetic  budgets. 

In  Fig.  3.2,  the  obvious  feature  common  to  all  six  volumes  Is  that 

there  is  a  dominant  balance  between  just  two  of  the  many  terms  included  in 


figure  3.2.  Su weary  of  depth-integrated  energy  budgets  for  six  regions 
(described  In  text).  Arrows  show  direction  of  fluxes  Into  or  out  of  boxes 
nuafcers  are  teres  Integrated  over  the  volume.  Pattern  Is  given  In  upper 
left  box  and  Is  same  for  others.  Sign  Is  appropriate  to  LHS  of  energy 
equation  for  all  teres.  Key:  KE  .  total  kinetic  energy  flux  and  Its 
divergence  over  the  voluae;  PH  .  pressure  work  teres;  HF  ■  heat  flux,  or 
conversions  between  potential  and  kinetic  energy;  WW  .  work  done  on  volue 
by  wind  stress;  FR  ■  dissipation  due  to  bottom  friction.  Residual 
lafcalance  Is  shown  In  box. 


the  budget.  (Notice  that  the  residuals  never  exceed  about  1/3  the 
magnitude  of  the  dominant  terms.)  In  the  accelerating  portion  of  the  jet, 
the  total  flux  of  kinetic  energy  out  of  the  region  is  nearly  double  that 
coming  in.  This  divergence  appears  to  be  balanced  by  a  net  amount  of  work 
being  provided  by  pressure-velocity  correlations;  of  secondary  importance 
is  the  fact  that  the  wind  stress  tends  to  accelerate  the  flow  throughout 
the  jet  region.  In  the  decelerating  jet,  the  roles  of  kinetic  energy  flux 
and  pressure  work  are  reversed:  the  former  provides  a  tremendous 
convergence  of  energy  which  is  radiated  mostly  north  and  south  via  the 
pressure  work.  The  secondary  terms  in  this  region,  which  are  a  third  the 
size  of  the  dominant  ones,  are  a  net  input  of  energy  by  the  wind  stress  and 
a  net  conversion  of  kinetic  to  available  potential  energy.  Because  they 
are  so  similar  in  magnitude,  it  is  tempting  to  suggest  that  all  the  wind 
work  is  going  directly  into  potential  energy,  but  of  course  there  is  no 
real  basis  for  drawing  this  conclusion.  It  is  possible  that  the  overall 
conversions  in  this  region  are  quite  complicated. 

The  balances  in  the  recirculation  regions  are  of  a  different  nature, 
for  here  conversions  to  and  from  potential  energy  play  a  dominant  role.  In 
the  eastern  portion  of  the  recirculation  region,  south  of  the  decelerating 
jet,  there  is  a  substantial  amount  of  pressure  work  exerted  on  the  region, 
which  is  primarily  balanced  by  a  sizable  conversion  of  kinetic  to  potential 
energy.  Of  secondary  importance  in  this  region  is  a  convergence  of  kinetic 
energy  flux,  due  primarily  to  a  large  flux  coming  across  the  northern 
boundary  from  the  jet.  In  the  western  recirculation  pressure  work 
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is  unimportant ,  and  there  is  a  conversion  from  potential  to  kinetic  energy  , 
which  is  returned  to  the  accelerating  jet  region  as  It  fluxes  across  the 
northern  boundary  of  this  area.  In  all  the  regimes,  of  course,  bottom 
friction  extracts  energy  from  the  system;  everywhere  except  in  the  jet 
regions,  wind  stress  acts  to  remove  energy  as  well,  usually  on  a  small 
scale  comparable  to  or  less  important  than  friction. 

The  Sverdrup  regimes  depart  from  strict  two-term  balances.  In  the 
western  portion,  kinetic  energy  is  fed  into  the  jet  to  the  west  and  north, 
resulting  in  a  divergence  that  is  balanced  half  by  the  pressure  work 
exerted  on  the  region  and  half  by  a  conversion  of  potential  to  kinetic 
energy.  The  eastern  Sverdrup  region  also  sees  a  conversion  of  potential  to 
kinetic  energy,  which  it  appears  to  export  to  surrounding  areas  via 
pressure  work.  The  secondary  terms  in  this  region,  however ,  are  nearly  50 
percent  as  large  as  the  dominant  ones:  there  is  a  convergence  of  kinetic 
energy  roughly  equal  in  magnitude  to  the  net  loss  of  energy  to  bottom 
friction. 

An  examination  of  the  Individual  energy  fluxes  across  the  boundaries 
of  all  the  regions  (Figure  3.3)  reveals  several  interesting  asymmetries, 
for  example,  that  the  accelerating  and  decelerating  jet  regimes  are  not 
"mirror  Images"  in  an  energetic  sense.  In  the  west,  mean  flux  of  mean 
kinetic  energy  accounts  for  only  about  half  the  net  divergence,  the  deficit 

being  made  up  mostly  by  the  terms  |j(u  u"*^  +  v  u'v')  +  ~(u  u'v'  +  v  v'^). 

Because  of  the  role  these  terms  play  In  the  individual  mean  and  eddy 
equations,  they  will  be  referred  to  as  the  mean/eddy  exchange  terms.  In 


Figure  3.3.  Summary  of  Independent  components  of  depth-integrated  kinetic 
energy  flux.  Arrows  show  directions  of  fluxes  Into  or  out  of  boxes.  Key: 
MM  •  mean  flux  of  MK£;  ME  »  mean  flux  of  EKE;  EE  *  eddy  flux  of  EKE;  EX 
«  flux  of  mean/eddy  exchange  terms  (see  text).  Divergence  of  various 
fluxes  are  listed,  and  total  divergence  Is  In  heavy  type.  Left  to  right, 
top  to  bottom,  are:  accelerating  and  decelerating  Jet;  western  and  eastern 
recirculation;  western  and  eastern  Sverdrup. 
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the  east,  the  convergence  of  kinetic  energy  flux  is  95  percent  due  to  the 
mean  fluxes.  The  eddy  flux  of  EKE  is  divergent  in  the  decelerating  as  well 
as  accelerating  part  of  the  jet,  so  that  it  acts  as  another,  smaller  scale 
"radiation  agent"  In  the  eastern  portion,  while  the  mean  flux  of  EKE  is 
everywhere  convergent  in  the  jet.  This  opposition  of  the  tendencies  for 
mean  and  eddy  fluxes  of  EKE  extends  to  other  parts  of  the  domain  and  is 
interesting  because  the  two  terms  are  similar  in  magnitude,  while  any 
linear  or  small-amp! itude  non-linear  theory  would  completely  ignore  the 
eddy  fluxes  of  EKE. 

Still  more  insight  may  be  gained  by  examining  separately  the  budgets 
for  mean  and  eddy  kinetic  energies.  In  this  context,  one  must  be  careful 
in  the  interpretation  of  the  mean/eddy  exchange  terms.  Harrison  and 
Robinson  (1978)  have  discussed  at  some  length  the  interpretation  of  these 
terms  in  open  regions.  They  classify  three  types  of  regions,  which  may  be 
summarized  as  follows:  1)  those  regions  in  which  the  conversion  of  mean  to 
eddy  kinetic  energy  implied  by  the  equation  for  K  is  approximately  equal 
and  opposite  to  the  conversion  of  eddy  to  mean  kinetic  energy  implied  by 
the  equation  for  K':  in  other  words,  the  divergence  of  the  mean/eddy 
exchange  terms  over  the  volume  is  nearly  zero;  2)  regions  in  which  the  two 
conversions  are  small  compared  to  other  processes  and  may  be  Ignored  in  the 
energy  budgets;  3)  regions  In  which  the  two  conversions  are  large,  but 
equal  and  opposite:  In  other  words,  the  devergence  of  the  mean/eddy 
exchange  terms  over  the  volume  is  considerable.  Only  in  regions  of  type  1 
or  2  Is  there  a  clear  interpretation  of  the  energy  budget  for  the  region  in 
terms  of  mean  to  eddy  conversion  processes.  Here  a  different  interpre- 


tation  is  adopted.  The  budget  for  EKE  over  a  region,  for  example,  may 
indicate  that  eddy  energy  is  increasing  at  the  expense  of  MKE;  yet  over  the 
same  region,  the  budget  for  MKE  may  suggest  that  there  is  a  conversion  of 
eddy  to  mean  kinetic  energy.  The  region  is  of  type  3  according  to  Harrison 
and  Robinson.  Indeed,  the  exchange  terms  as  they  appear  in  the  two 
equations  are  not  equal  and  opposite,  but  differ  by  a  divergence  of  the 
mean/eddy  exchange  fluxes  across  the  boundary,  which  redistribute  eddy  and 
mean  kinetic  energy  regionally,  where  it  may  then  appear  as  a  conversion 
from  one  type  to  the  other.  Because  the  major  thrust  of  this  work  is  to 
examine  the  divergence  of  fluxes  across  the  open  boundaries,  rather  than  to 
evaluate  the  conversions  between  eddy  and  mean  kinetic  energies  as  in 
Harrison  and  Robinson  (1978),  a  type  3  region  is  not  considered  a  "problem" 
in  the  present  work. 

Figures  3.4  through  3.6  display  the  interdependent  budgets  for  mean 

and  eddy  kinetic  energy  in  the  six  regimes,  and  are  to  be  interpreted  as 

follows.  Consider  Fig.  3.4a  as  an  example.  Here,  in  the  accelerating  jet, 

g 

the  net  divergence  of  MKE  fluxes  over  the  area,  equal  to  48.2  x  10  J/s, 

implies  that  mean  kinetic  energy  must  be  supplied  to  the  region  at  this 

rate  by  other  mechanisms  in  order  to  maintain  the  budget.  Moreover,  there 

g 

is  a  constant  conversion  to  EKE  of  32.0  x  10  J/s,  which  is  interpreted 

as  being  fluxed  out  of  the  region,  and  a  net  loss  to  dissipation  of  1.4  x 

10  J/s.  On  the  supply  side,  mean  potential  energy  is  being  converted  to 

g 

MKE  at  the  rate  of  5.5  x  10  J/s  and  wind  is  supplying  eneregy  at  16.6  x 
109  J/s.  There  is  still  a  net  deficit  of  59.5  x  109  J/s,  which  must  be 
supplied  by  the  pressure  work  divergence,  the  only  term  remaining  in  the 


Figure  3.4.  Schematic  of  Interconnected  mean  and  eddy  energy  budgets  In 
a)western;  b)eastern  Jet  region.  See  text  for  explanation  of  how  to 
Interpret.  Key:  K  *  MKE;  K'  «  EKE;  P  «  mean  potential  energy;  P'  « 
eddy  potential  energy;  f >  work  done  by  wind  stress;  5  *  dissipation  due 
to  bottom  friction  acting  on  mean  flow;  D'  »  dissipation  due  to  bottom 
friction  acting  on  eddy  terms;  PW  »  pressure  work. 


Western  Jet 


AD-A154  753  HORIZONTAL  AND  VERTICAL  STRUCTURE  OF  VELOCITY  POTENTIAL 
^  VORTICITY  AND  ENE. . (U>  HOODS  HOLE  OCEANOGRAPHIC 

INSTITUTION  HA  H  H  HALL  FEB  85  MHOI-85-16 
UNCLASSIFIED  N00014-82-C-0019  F/G  8/3 


2, 

NL 


Same  as  Figure  3.4,  but  for  recirculation  region 


» 


94 


budget  for  MKE.  Similar  analysis  of  the  eddy  budget  suggests  the  pressure 

9  9 

work  is  supplying  17  .8  x  10  J/s,  for  a  total  of  77.3  x  10  J/s.  In 

g 

parentheses  is  the  calculated  pressure  work  divergence  of  89.7  x  10  J/s, 
and  the  two  values  differ  by  the  residual  Imbalance  shown  in  Fig.  3.2. 

Still  concentrating  on  the  jet  region,  observe  that  the  fluxes  and 
exchanges  Involved  In  maintaining  the  MKE  budget  are  nearly  an  order  of 
magnitude  larger  than  those  In  the  EKE  budget.  In  particular,  one  Is 
forced  to  interpret  most  of  the  pressure  work  divergence  as  being  due  to 
the  correlations  between  the  mean  pressure  field  and  a  mean  ageostrophlc 
velocity  field.  That  interpretation  runs  counter  to  the  notion  that  the 
pressure  work  Is  radiating  energy  via  wavelike  structures,  through  terms 
like  p'u'.  However,  recall  that  the  mean  ageostrophlc  velocity  field 
depends  on  time-dependent  quantities  arising  from  Inertial  accelerations  In 
the  momentum  equations.  In  Hie  decelerating  jet  (Fig.  3.4b),  going  beyond 
the  dominant  balance,  there  are  two  more  interesting  features.  First, 
there  is  a  sizable  conversion  of  mean  kinetic  to  potential  energy,  the 
mechanism  proffered  by  Fofonoff  and  Hall  (1983)  to  account  for  downstream 
decreases  In  kinetic  energy  flux.  Second,  In  the  eddy  budget,  there  Is  an 
energy  pathway  suggestive  of  barotroplc  Instability,  as  MKE,  Interpreted  as 

9 

being  fluxed  into  the  region.  Is  converted  to  EKE  at  a  rate  of  13.1  x  10 
J/s.  Meanwhile,  the  direction  of  energy  flow  between  eddy  kinetic  and 
potential  energies  Is  opposite  of  what  would  be  anticipated  were  baroclinlc 
instability  operative  In  this  region. 

The  story  told  by  the  budgets  in  the  recirculation  regime  Is  quite 
different  (Figure  3.5).  In  the  first  place,  mean  and  eddy  fluxes  and 
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exchanges  are  of  like  magnitude.  In  the  east  there  Is  a  continuation  of 

the  conversion  of  MKE  to  mean  potential  energy.  However,  eddy  budgets 

throughout  the  recirculation  now  display  an  energy  pathway  Implying 

barocllnlc  Instability,  where  eddy  potential  energy  Is  converted  to  EKE  (at 
g 

rates  of  3-4  x  10  J/s),  and  EKE  Is  converted  to  MKE  (at  rates  of  3-3.5  x 
g 

10  J/s).  Curiously,  throughout  the  recirculation,  the  eddy  budget 
implies  that  pressure  work  due  to  eddy  terms  radiates  energy  into,  rather 
than  out  of,  the  region.  The  mean  budget  in  the  west  shows  the  other  half 
of  the  inertial  recirculation  of  energy  as  potential  energy  --  that  Is, 
there  is  now  a  conversion  of  mean  potential  to  mean  kinetic  energy,  which 
helps  to  accelerate  the  flow. 

Figure  3.6,  showing  the  Sverdrup  regions,  reveals  little  more  than 
has  already  been  learned,  save  that  the  energy  budgets  here  are  dominated 
by  mean  fluxes  and  exchanges.  In  the  west,  there  Is  again  acceleration  of 
the  mean  flow  as  mean  potential  energy  Is  converted  to  MKE. 

Comparison  of  these  results  with  Harrison  and  Robinson's  (1978) 
results  Is  limited,  for  several  reasons.  In  the  first  place,  they  analyzed 
a  single  gyre  model  with  the  eastward  jet  flowing  along  the  northern  wall 
rather  than  freely  in  mid-ocean,  and  the  basin  was  only  2000  km  by  2000 
km.  Second,  they  did  not  examine  the  accelerating  and  decelerating 
portions  of  the  eastward  jet  separately,  so  the  very  large  fluxes  of  mean 
kinetic  energy  in  the  jet  are  of  no  Importance  In  their  budgets. 
Volume-averaged  magnitudes  of  the  mean  to  eddy  conversion  are  roughly 
comparable  for  the  two  investigations ,  but  because  of  the  importance  of  the 
mean  flow  terms  In  the  budgets  presented  here,  most  of  the  regions  In  this 


analysis  would  be  considered  type  2  (with  the  exception  of  the  western 
recirculation  and  possibly  the  western  jet).  Notice  that  if  the  regions 
had  not  been  divided  into  their  eastern  and  western  portions,  conversions 
in  the  mean  energy  equations  would  tend  to  cancel  out,  and  the  mean/eddy 
conversions  would  become  relatively  more  Important;  then  most  of  the 
regions  would  be  type  3. 

Finally,  it  is  worth  considering  the  energy  picture  for  the  entire 

domain,  shown  in  Figure  3.7.  Although  there  may  still  be  fluxes  across  the 

open  boundaries,  the  very  strong  Internal  fluxes  and  exchanges  tend  to 

cancel  out,  and  forcing  and  dissipation  assume  much  more  important  roles. 

In  the  complete  absence  of  fluxes  across  the  boundaries  of  the  domain,  the 

two  would  necessarily  be  found  equal  and  opposite.  Even  in  this  open 

domain,  that  balance  is  the  most  obvious  part  of  the  budget,  as  bottom 

g 

friction  removes  energy  at  a  rate  of  27.1  x  10  J/s  while  wind  puts 

g 

energy  in  at  34.9  x  10  J/s.  Also  note  that  most  of  the  energy  is 
dissipated  through  eddy  terms:  In  spite  of  the  dominance  of  mean  terms  in 
the  localized  budgets,  ultimately  the  total  system  depends  on  the  eddies  in 
a  direct  way  to  balance  the  energy  budget.  The  fact  that  the  system  is 
entirely  forced  by  the  wind  stress,  which  directly  affects  only  the  mean 
kinetic  energy  budget,  means  that  there  is  an  Inescapable  coupling  between 
mean  and  eddy  fields,  as  Harrison  (1979)  concluded. 

Individually,  the  mean  and  eddy  kinetic  energy  fluxes  over  the  domain 
appear  to  be  convergent.  However,  taking  into  account  the  mean/ eddy 
exchange  terms,  it  is  found  that  the  total  energy  fluxes  are  in  fact 
divergent;  there  is  also  a  net  conversion  of  kinetic  to  potential  energy. 
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Although  the  excess  of  forcing  over  dissipation  makes  up  part  of  the 
deficit,  notice  that  energy  must  still  be  Imported  via  pressure  work  terms 
--  curiously,  this  occurs  in  the  eddy  field,  so  that  waves  from  regions 
outside  this  domain  evidently  radiate  energy  into  itJ 

This  discussion  of  energy  budgets  in  the  numerical  model  must  be 
concluded  with  several  warnings.  As  pointed  out  at  the  beginning  of  this 
chapter,  the  analysis  has  treated  the  time-averaged  flows  in  the  vicinity 
of  an  intense  jet  which  meanders  strongly  on  an  instantaneous  basis.  Thus, 
this  analysis  addresses  not  so  much  the  energetics  of  the  jet  as  the 
energetics  of  that  region  in  which  the  jet  is  usually  found.  The  jet 
itself  may  have  some  average  "structure”  that  it  carries  about  as  it 
meanders.  The  deviations  from  such  an  average  are  lumped  together  with  the 
meandering  to  comprise  the  eddy  portion  in  an  Euler ian  time-average  of  the 
flow  field.  Without  separating  various  effects,  one  cannot  carry  out  a 
phenomenological  investigation  of  the  energetics,  identifying  time  and 
space  scales,  energy  sources,  and  propagation  characteristics  of  those 
features  most  important  to  maintaining  the  energy  budgets. 

Rather,  a  zeroth  order  evaluation  of  the  budgets  has  been  presented, 
with  the  idea  that  it  can  provide  guidelines  for  future  work  of  this  type 
either  on  numerical  models  or  in  the  ocean.  Of  the  six  energetic  regimes 
that  have  been  identified,  balances  in  both  the  accelerating  and 
decelerating  jet  regions  are  dominated  by  mean  quantities.  In  particular 
the  downstream  convergence  of  kinetic  energy  in  the  decelerating  jet  is 
balanced  primarily  by  a  mean  ageostrophic  flow  against  the  pressure 
gradient,  which  in  turn  implies  some  conversion  of  kinetic  to  available 


potential  energy.  Mean  and  eddy  quantities  are  of  equal  Importance  in  the 
recirculation  regimes,  but  the  mean  flow  again  dominates  energetic  balances 
in  the  Sverdrup  regions.  That  so  much  energy  is  exported  from  the  jet 
region  implies  that  the  jet  is  not  energetically  isolated  from  its 
environs,  so  channel  models  probably  are  not  a  good  means  for  investigating 
the  jet's  dynamics  and  energetics.  Indeed,  Harrison  (1979)  has  pointed  out 
the  necessity  for  local  energetic  analyses  of  regions  with  open  boundaries 
in  numerical  models.  The  conversions  between  kinetic  and  potential  energy 
in  both  the  jet  and  recirculation  regions  are  consistent  with  Fofonoff  and 
Hall's  (1983)  conclusions  regarding  energy  conversions  evaluated  from  the 
Gulf  Stream  '60  data.  Finally,  the  energetic  budgets  suggest  that 
barotropic  instability  may  be  occurring  in  the  jet  itself,  while  baroclinic 
instability  energy  pathways  appear  only  in  the  recirculation  regime. 
Energetically,  the  jet  and  recirculations  are  evidently  parts  of  an 
inseparable  whole;  further  observational  and  numerical  investigations 
should  be  designed  with  that  result  in  mind. 


Chapter  4.  Cross-sections  of  potential  vorticity  in  the  Gulf  Stream 

4.1  Introduction 

Potential  vorticity  of  a  fluid  is  a  dynamically  important  quantity 
because  it  is  nearly  conserved  following  fluid  parcels  if  forcing, 
dissipation,  and  mixing  are  sufficiently  small.  Fofonoff  (1962)  and 
Stonriel  (1965),  among  others,  have  suggested  modelling  the  Gulf  Stream  as  a 
layer  of  water  with  uniform  potential  vorticity  ,  which  corresponds  to  q 
being  constant  on  temperature  (or  density)  surfaces,  and  simplifies  the 
mathematics  of  the  problem.  On  the  other  hand.  If  potential  vorticity  is 
not  constant  on  surfaces  of  constant  temperature,  then  if  temperature  is 
also  conserved  following  fluid  parcels,  paths  of  flow  may  be  determined 
from  the  intersection  of  surfaces  of  constant  potential  vorticity  and 
constant  temperature.  Thus  it  is  important  to  determine  how  potential 
vorticity  is  distributed  on  isothermal  surfaces.  Recent  investigators 
mapping  large  scale  potential  vorticity  fields  in  the  oceans  (McDowell  et. 
al . ,  1982;  McCartney,  1982)  have  recognized  this  point,  but  their  work  has 
been  restricted  to  regions  of  relatively  quiet  flow,  where  q  can  be 
approximated  by  fez.  From  the  GUSTO  data  set  it  is  possible  also  to 
calculate  the  relative  vorticity  contributions  -uyez  and  uzTy.  The 
resulting  potential  vorticity  section,  when  compared  with  the  temperature 
structure  in  the  Stream,  should  offer  some  insight  as  to  whether  a  uniform 
potential  vorticity  model  of  the  Gulf  Stream  is  indeed  appropriate. 

The  mean  potential  vorticity  field  of  a  quasi-geostrophic  flow  regime 
also  determines  the  instability  properties  of  the  flow.  In  particular,  the 


following  constitute  a  necessary  condition  for  Instability  of  the  flow  to 
infinitessimal  perturbations  (see  Pedlosky,  1979,  for  a  good  discussion): 
In  a  zonal ly  bounded  flat-bottomed  channel  In  which  a  basic  flow 


U  =  U(y,z)  only  exists,  then  If: 

1) The  potential  vortlcity  gradient  qy  changes  sign  within  the 
flow;  or 

2) qy  is  somewhere  the  opposite  sign  of  Uz  at  the  surface,  or  the 
same  sign  as  U2  at  the  bottom;  or 

3) UZ  at  the  surface  Is  the  same  sign  as  at  the  bottom,  then  the 
flow  may  be  unstable. 

The  easiest  of  these  to  test  from  data  is  whether  or  not  qy  changes  sign 
somewhere  in  the  flow;  if  It  does,  then  the  flow  under  consideration  Is 
potentially  unstable.  8ecause  isothermal  surfaces  depart  markedly  from 
horizontal  surfaces,  the  gradient  of  q  on  both  will  be  examined,  and  the 
results  compared. 

One  test  for  the  validity  of  a  numerical  model  is  how  well  It 
reproduces  the  observed  potential  vortlcity  distribution  of  the  real 
ocean.  Since  the  model  jet  is  highly  energetic.  It  is  especially  important 
that  this  region  be  reproduced  realistically.  Thus,  It  Is  useful  to 
compare  the  observed  potential  vortlcity  distribution  with  that  found  in 
the  numerical  model  jet.  Because  the  layers  In  the  model  are  Isopycnal,  q 
in  each  layer  should  again  be  compared  with  q  along  Isothermal  surfaces  In 
the  data. 


4.2  Potential  vorticity  section  from  the  mooring  data 


It  was  established  in  Chapter  2  that  the  appropriate  form  for 
potential  vorticity  at  the  mooring  site  is 


q  =  (f  -  uy)  oz  +  u2  Ty. 

To  construct  an  average  potential  vorticity  section,  the  mooring  data  were 
first  divided  into  temperature  bins  determined  by  Tgyg.  From  Tg^g  = 

4.5*C  to  7.5°C,  and  from  Tgyg  =  13.5*C  to  17°C,  the  bins  are  .5#C  wide; 
from  Tgyg  =  7.5°C  to  13.5*C,  they  are  1*C  wide,  giving  a  total  of  19 
bins.  The  average  along  stream  velocity  shear  between  the  400  and  700  m 
instruments  was  then  calculated  for  each  bin,  and  from  these  values,  an 
average  Ty  was  obtained  for  each  bin.  (Since  only  rotated  velocities 
will  be  discussed,  hat  notation  will  be  dropped  for  this  chapter.)  This 
average  function  Ty(T)  at  575  db  was  then  integrated  across  all  values 
of  Tg75  to  obtain  y  as  a  function  of  Tg^g.  This  procedure  was 
completely  analogous  to  that  used  in  Chapter  2  to  obtain  a  continuous 
function  y(T),  only  there  the  function  Ty(T)  was  determined  continuously, 
as  two  best  fit  lines,  rather  than  discretely.  The  values  obtained  both 
ways  are  compared  in  Table  4-1,  where  the  origin  for  y  is  at  13.0°C. 
Except  at  the  very  coldest  estimate,  the  comparison  is  excellent.  Next 
along-stream  velocities  at  all  depths  were  sorted  into  temperature  bins  and 
an  average  velocity  for  each  bin  at  each  depth  was  obtained:  these  values 
were  assigned  to  the  mid-points  of  the  bin.  In  order  to  assure  that 
along-stream  velocity  could  go  to  zero  at  the  edges  of  the  average  profile, 
average  eastward  (rather  than  along-stream)  velocity  was  calculated  in  the 
warmest  and  coldest  bins.  The  along-stream  velocities  were  center 


y (analytical) 
(km) 


y(from  average) 
(km) 


T  (“C) 
575'  J 


4.5 

97.5 

123.5 

5,0 

78.9 

82.2 

5.5 

66.1 

65.3 

6.0 

56.4 

55.0 

6.5 

48.5 

46.6 

7.0 

41.9 

40.9 

7.5 

36.3 

35.7 

8.5 

26.8 

27.5 

9.5 

19.2 

20.5 

10.5 

12.7 

13.8 

11.5 

7.2 

8.3 

12.5 

2.3 

2.6 

13.5 

-2.4 

-2.6 

14.0 

-5.3 

-5.4 

14.5 

-8.7 

-8.9 

15.0 

-13.1 

-13.0 

15.5 

-19.0 

-19.8 

16.0 

-28.3 

-30.7 

16.5 

-51.1 

-50.8 

Table  4-1. 

Intercomparison  of  horizontal 

coordinate 

values  obtained  analytically  from  equations  2-5,  2-6 
with  values  obtained  from  averaging  procedure  described 
in  Chapter  4.  In  both  cases,  y  is  set  equal  to  zero 
at  T  =  13°C. 

D  /-> 


T  Bin  Y  _  0y(575)  0y<875)  0y(1175  0y(2OOO)  Ov(4000) 

575  center 


(°C) 

(km) 

(All 

shears 

have  units  of  10  5  s 

*) 

4.593 

4.997 

5.340 

5.899 

6.400 

6.982 

7.658 

8.463 

9.441 

10.405 

11.475 

12.410 

13.276 

14.013 

14.590 

15.142 

15.663 

16.113 

16.374 

85.9 

-.688 

-.  398 

-.391 

-.403 

-.387 

79.1 

-.814 

-.459 

-.415 

-.455 

-.469 

62.2 

-.623 

-.199 

-.147 

-.235 

-.262 

51.9 

-.351 

.447 

.064 

.296 

.427 

44.5 

-1.138 

-.145 

-.402 

-.143 

-.046 

37.7 

-1.181 

-.292 

-.384 

-.325 

-.241 

31.1 

-1.737 

-.626 

-.250 

-.174 

.015 

24.4 

-1.394 

-.783 

-.307 

-.109 

.271 

17.5 

-1.196 

-.228 

.111 

.270 

.324 

11.3 

-1.917 

■1.100 

-.191 

-.140 

-.218 

5.5 

-1.615 

•1.303 

-.414 

-.215 

-.149 

0.7 

-.799 

-.826 

-.231 

.016 

.069 

-3.6 

1.950 

-.261 

.140 

-.207 

-.314 

-7.5 

2.726 

.117 

.547 

.065 

-.214 

-11.8 

2.720 

.170 

.720 

.442 

.371 

-18.1 

2.700 

.952 

.852 

.507 

.374 

-28.6 

1.312 

.647 

.264 

.156 

-.039 

-41.9 

1.383 

.927 

.209 

.200 

.003 

16.775 

-49.6 

1.958 

1.496 

.452 

.388 

.801 

Table  4-2. 

Temperature  bins 

arrived 

at  for 

averaging  procedure 

described  in  text,  along  with  the  value  of  y  halfway  between  the  end 
point  values  for  each  bin.  Horizontal  shear  of  long-stream  velocity 
is  from  center-differencing  average  u  values  from  Hall  and  Bryden 
(1984) ,  except  first  and  last  values,  which  could  not  be  centered. 
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flnite-differenced  to  obtain  estimates  of  au/ay;  because  the  calculated 
values  of  u  are  not  evenly  spaced  in  y,  this  step  required  calculating 
the  new  values  of  y  (and  hence  T575)  to  which  au/ay  values  applied. 
Correspondingly,  new  temperature  bins  were  also  defined,  as  shown  in  Table 
4-2,  along  with  values  of  au/ay.  In  the  new  temperature  bins,  average 
values  of  vertical  temperature  gradient  and  vertical  shear  of  along-stream 
velocity  were  calculated  from  the  mooring  data  for  the  top  three  standard 
levels.  Taken  all  together,  it  is  then  possible  to  obtain  estimates  of  q 
at  three  points  in  the  vertical  and  19  In  the  horizontal,  from 

Simple  linear  interpolation  or  extrapolation  was  used  to  obtain  continuous 
values  of  q  in  the  vertical  which  were  then  contoured  to  provide  the  map 
of  dashed  lines  shown  in  Figure  4.1. 

The  solid  lines  in  Figure  4.1  are  isotherms.  This  average 
temperature  cross  section  was  obtained  as  follows.  At  each  y,  values  of 
T575  are  already  known;  average  values  of  Tg j5  and  T^5  were 
calculated  for  the  corresponding  temperature  bins.  Then  the  analytic  fits 
for  3©/az  at  each  level  (Raymer,  Spencer  and  Bryden,  1984)  were  used  to 
integrate  up  or  down  to  obtain  nearly  continuous  vertical  profiles  at  each 
y.  (The  vertical  gradient  functions  change  form  at  700  and  1000  m.)  The 
temperature  values  were  then  contoured  to  yield  the  map  of  solid  lines  In 
Figure  4.1. 

Figure  4.1  requires  careful  consideration.  Recall  that  the  map  of 
q  was  constructed  from  only  3  points  in  the  vertical.  Thus,  for  example. 
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the  mid-depth  maximum  at  the  southern  edge  is  defined  by  nine  points:  at 
the  three  southernmost  values  of  y,  estimates  of  q  at  875  db  are  larger 
than  those  above  and  below.  The  strong  central  maximum  in  the  upper  layer 
results  from  estimates  only  at  575  db,  but  it  is  a  well-defined  trend  at 
this  level.  Comparison  with  a  section  from  Watts  (1983)  as  well  as  with 
sections  from  the  numerical  model  will  show  that  certainly  the  latter 
feature  and  probably  the  former  feature  are  both  real . 

On  the  cold  side  of  the  Stream,  for  Tg^g  <  11°C,  isotherms  tend 
to  parallel  isostrophes  (lines  of  constant  potential  vorticity).  Isotherms 
for  T  >  9°C  heading  northward  into  the  core  of  the  Stream  must  cross 
isostrophes  almost  perpendicularly,  going  from  lower  values  in  the  south  to 
higher  values  in  the  north.  In  addition,  most  of  these  isotherms  pass 
through  the  weak  relative  minimum  south  of  the  core.  Isotherms  for 
T  <  8°C,  on  the  other  hand,  tend  to  lie  on  the  same  isostrophe  at  both  the 
southern  and  northern  edges  of  the  section;  that  they  do  not  exactly 
parallel  isostrophes  in  between  the  endpoints  may  be  due  to  the 
uncertainties  involved  in  creating  the  cross-section. 

To  complement  Figure  4.1,  q  is  shown  as  a  function  of  y  for 
selected  isotherms  in  Figure  4.2,  which  is  derived  simply  by  reading  points 
off  Figure  4.1.  These  curves  may  be  compared  with  the  values  derived  from 
the  analytic  expressions  at  575  db  and  tabulated  in  Table  2-4:  both 
results  display  a  strong  "wall"  in  q  near  the  jet  axis,  a  minimum  just 
south  of  that,  and  an  increase  and  subsequent  leveling  off  as  we  head  into 
warmer  waters,  on  isotherms  T  =  ll'C  and  14*C.  The  e-effect  has  not  been 
explicitly  included  in  the  estimates  of  q,  and  it  is  now  shown  that  the 


omission  introduces  very  small  errors.  In  the  first  place,  By  «  f  and 
even  u^:  values  of  8  =  2  x  10"^m'*s~*  and  y  =  100  km  yield  maximum  values 
of  sy  =  2  x  10~®s_*.  The  e-effect  could  still  be  important  to  the 
gradient  of  q.  If  e  were  included  in  the  definition  of  q,  two  extra 
terms  eTz  and  ey  Tzy  would  appear  in  qy.  Estimates  of  qy 
from  Figure  4.2  are  compared  with  estimates  of  eTz  in  Table  4-3  (scale 
analysis  suggests  that  By  Tzy  would  be  the  same  order  as  bTz).  The 
latter  tends  to  be  one  to  two  orders  of  magnitude  smaller  than  aq/ay,  and 
clearly  would  not  alter  the  major  features  displayed  in  Figure  4.1. 

Values  of  q  on  the  remaining  Isotherms  are  less  easily 
interpreted.  Endpoints  of  the  5,  7,  and  9*C  isotherms  tend  to  lie  on  equal 
values  of  q,  but  all  three  exhibit  a  minimum  near  y  =  0.  Inspection  of 
Fig.  4.1  shows  that  while  the  relative  minimum  between  y  =  0  and  -25  km  is 
well-defined  above  800  m  depth  or  so,  it  decays  considerably  below  that,  so 
that  the  appearance  of  a  minimum  on  the  5*C  and  7*C  isotherms  in  Fig.  4.2 
may  be  an  artifact  of  the  contouring  in  Fig.  4.1. 

4.3  Comparison  with  past  results 

Watts  (1983)  discusses  the  potential  vorticity  distribution  across 
the  Gulf  Stream,  its  relation  to  the  general  circulation,  and  Its 
implications  for  instability.  In  particular,  he  points  out  the  following 
Important  contrast:  on  the  one  hand,  the  uniformity  of  q  In  the  gyre 
interior  taken  with  conservation  of  q  along  streamlines,  suggests  that 
potential  vorticity  in  the  Gulf  Stream  ought  to  be  uniform  as  well,  since 
streamlines  from  the  gyre  feed  the  Stream.  On  the  other  hand.  Watts  says 
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-5  y  2T  Tz  q  <*y  STz 

:)  (km)  (m)  (°C/m)  do”7  °C/m/s)  (°C/m2/s)  (°C/m2/s) 


4.96  85.9  405  1.62  x  10_2 
8.19  31.1  672  1.45  x  io'2 

12.86  0.7  659  2.51  x  io'2 
14.30  -7.5  715  2.30  x  io"2 
16.47  -49.6  880  2.15  x  io'2 

12.86  0.7  532  2.33  x  io"2 
14.30  -7.5  597  1.99  x  io'2 
15.90  -28.6  698  1.90  x  10~  2 
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10.17 
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9.86 

15.91 
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14.48 


-6.8  x  10 


•12 


3.1  x  io 


4.4  x  io 


-1.4  x  io 


•l  l 


1 1 


8.2  x  io"11 


-1.1  x  io 


11 


2.8  x  10 


1  3 


‘1  3 


4.8  x  lo“13 


4.4  x  io 


“1  3 


4.1  x  io 


,“13 


4.4  x  io 


“IS 


3.8  x  10 


1  3 


3.6  x  io 


“IS 


Table  4-3.  Comparison  of  S  effect  with  total  potential  vorticity 

gradient.  Values  of  q  for  isotherms  come  essentially  from 

Figure  4.1.  Values  of  T  are  from  analytic  fits  of  Raymer  et^  al . 

2  ^ 

(1985).  8  is  taken  as  1.9  x  10  11  n  1  s  ’. 


that  "for  all  isopycnal  layers  In  and  above  the  main  thermocline, 
q  ~  f  Tz  Is  several  times  higher  in  the  Slope  Water  than  in  the  Sargasso 
Sea;"  hence,  a  strong  jump  in  potential  vorticity  must  occur  across  the 
Stream.  The  same  features  were  noted  much  earlier  by  Stommel  (1965),  who 
demonstrated  the  uniformity  of  potential  vorticity  across  the  Stream,  from 
the  Sargasso  Sea  to  the  Inshore  edge  of  the  Stream,  where  there  is  a  sharp 
discontinuity  in  q. 

Two  of  Watts'  (1983)  figures  may  be  compared  with  Figures  4.1  and 
4.2.  His  Figure  12  Is  a  section  of  potential  vorticity  across  the  Gulf 
Stream  near  73* W,  from  the  surface  to  800  m  depth,  with  Isotherms 
superimposed.  Curiously,  below  400  m  in  this  figure,  isotherms  tend  to 
parallel  isostrophes  consistently,  while  above  400  m  features  very  similar 
to  those  seen  in  Figure  4.1  may  be  found.  In  particular  there  is  a  strong 
maximum  in  q  lying  just  around  the  axis  of  the  Stream,  with  q 
increasing  upward.  This  feature  corresponds  to  the  high  q  values  lying 
between  y  =  0  and  25  km,  at  depths  of  400-700  m  In  Fig.  4.1.  Just  south  of 
the  maximum  is  an  Intermediate  minimum  reaching  down  from  above,  much  like 
the  minimum  in  Figure  4.1  lying  between  y  =  0  and  -25  km,  at  depths  of 
about  500-800  m.  In  general,  the  highest  values  Watts  finds  for  q  below 
400  m  are  in  the  range  2  -  3  x  10”^#C  m~*  s"\  slightly  higher  than  the 
values  calculated  here. 

Watts'  Figure  15,  which  shows  potential  vorticity  across  the  Stream 
for  the  12°-17°C  layer  near  69*W,  may  be  compared  with  the  14°C  Isotherm  in 
Figure  4.2.  The  latter  shows  this  Isotherm  rising  from  its  minimum 
value  of  q  =  1.2  x  10"®*C  m“*  s”*  to  its  maximum  value  (below  400  m)  of 


2  x  10  Cm  s  about  10  km  northward;  Watts'  Figure  15  shows  an 
increase  from  a  similar  minimum  of  1  -  1.2  x  10~6°C  m" *  s"1  to  a  maximum 
of  around  4  x  10"®8C  m"*  s'*  some  40  km  northward  (where  this  layer  lies 
well  above  400  m) .  Moving  southward,  both  figures  show  a  slight  increase 
and  subsequent  decrease  in  q  over  similar  horizontal  scales  (50-60  km). 
Watts'  figures  were  constructed  partly  from  recent  hydrographic  data  and 
partly  from  a  section  taken  by  Warren  and  Volkmann  (1968),  when  deep 
reference  velocities  were  measured  by  floats.  The  present  evaluation  has 
the  advantage  of  direct  velocity  measurements  throughout  the  water  column, 
as  well  as  greater  coverage  in  time,  so  that  Fig.  4.1  may  be  considered  an 
average  and  not  a  synoptic  section. 

Now  consider  the  issues  raised  in  the  introduction  to  this  chapter. 

A  two-layer  model  for  a  free  zonal  inertial  jet  with  uniform  potential 
vorticity  is  governed  by  the  following  equations: 


y  <  0 


h  =  0  at  y  =  0;  h  -->  hQ  as  y  -->  -  do 


h  =  u  =  0  y>0 

where  h  is  the  upper  layer  depth,  hQ  is  the  quiet  upper  layer  depth  far  to 


the  south,  and  the  lower  layer  is  at  rest.  These  may  be  solved  to  obtain: 
h  =  ho(l  -  exy),  u  =  (ghQ)1/2  exy,  x  =  ( fQ )/ ( ghQ ) 1/2 


Thus,  the  zonal  velocity  jumps  discontinuously  from  zero  to  its  maximum 
value  where  the  upper  layer  depth  goes  to  zero,  then  decreases 
exponentially  across  the  jet  towards  the  Interior,  and  the  flow  in  the  jet 


Is  everywhere  anti -cyclonic,  so  that  this  model  cannot  be  expected  to 
describe  the  cyclonic  portion  of  the  observed  jet.  In  isothermal  layers  at 
temperatures  lying  above  those  observed  at  the  mooring,  there  evidently  is 
a  tendency  towards  uniform  q  in  the  anti -cyclonic  portion  of  the  Stream 
(Watts,  1983;  Stommel ,  1965);  between  y  =  -25  and  -50  km  in  Fig.  4.1,  there 
is  a  suggestion  that  the  isotherms  for  T  =  10  to  13  ’C  may  be  starting  to 
parallel  isostrophes  as  the  interior  is  approached.  Thus  it  is  likely  that 
the  two-layer  inertial  jet  model  gives  a  reasonable  rendition  of  the  steady 
(average)  cross-stream  structure  of  velocity  in  that  portion  of  the  Stream. 

Investigation  of  the  potential  vorticity  gradient  qy  can  lead  to  a 
fuller  understanding  of  the  dynamics  of  the  Stream,  beyond  the  prediction 
of  average  cross- stream  structure.  Two  features  of  qy  are  of 
predominant  importance.  One  is  the  possibility  that  qy  changes  sign 
across  the  Stream;  the  other  is  the  strength  of  the  gradient  near  the 
Stream  axis.  Each  point  will  be  addressed  in  turn. 

It  has  been  noted  that  a  necessary  condition  for  the  baroclinic 
instability  of  a  zonal  flow  under  certain  conditions  is  that  qy  change 
sign  somewhere  in  the  basic  flow.  This  condition  can  be  easily  tested  with 
data,  but  one  must  be  cautious  in  applying  these  criteria  to  the  potential 
vorticity  and  velocity  fields  at  the  mooring  site,  for  various  reasons. 

The  criteria  are  meant  to  apply  to  quasl-geostrophlc  flow,  and  although  the 
Rossby  number  of  flow  at  the  mooring  site  is  generally  £0.3,  the 
non-quasl-geostrophic  term  uzTy  Is  as  much  as  25  percent  the  size  of 
q,  and  its  gradient  may  be  as  large  as  qy. 


Values  of  potential  vorticity  q  on  the  575  db  surface  as  a 
function  of  cross- stream  distance  y  at  the  mooring  site,  and 
individual  components  making  up  total  potential  vorticity. 

Key:  dashed  line  =  total  q;  solid  line  =  f«z;  dash-dotted 
line  =  relative  vorticity  (-u.oj;  dotted  line  =  u,T  . 
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Figure  4.3  shows  potential  vorticlty  and  its  components  at  a 
horizontal  level,  575  dbar.  There  is  a  pronounced  maximum  near  the  core  of 
the  Stream,  as  in  Fig.  4.1,  and  moving  southward,  a  rather  weak  relative 
minimum  between  y  =  -25  and  -50  km.  Notice  that  <  0  to  the  north  of 
the  maximum  is  misleading  in  this  figure;  a  glance  at  Figure  4.1  shows 
that  q  tends  to  become  uniform  on  Isotherms  in  this  region.  The  vertical 
and  horizontal  relative  vorticity  contributions  combine  to  affect  the 
potential  vorticity  profile  in  two  ways.  North  of  the  maximum,  they  very 
nearly  offset  each  other  so  that  q  ~  f  Tz.  However,  passing  through  the 
Stream  axis,  where  u^  changes  sign,  the  combined  strongly  negative 
values  of  -uy  Tz  and  uz  Ty  are  enough  to  offset  significantly  the 
maximum  in  q  from  that  in  f  Tz,  and  to  produce  the  weak  minimum  just 
south  of  there.  In  the  region  where  qy  fs  slightly  negative,  its 

1  p  _  O  I 

average  value  is  -4.5  x  10“  C  m  s  ;  over  the  same  region,  b  would 
contribute  to  the  gradient  a  term  of  the  size  (2.8-4. 5)  x  10“13°C  m“2  s"1, 
an  order  of  magnitude  smaller.  Moreover,  recall  that  isotherms  are  sloping 
downward  to  the  south,  thus  passing  through  the  minimum  and  into  higher 
values  of  q  than  appear  at  y  =  22  km  in  Figure  4.3  (see  Figure  4.1).  It 
would  appear  that  qy  changes  sign  on  a  horizontal  surface  as  well  as  on 
isothermal  surfaces.  Deeper  down,  qy  tends  to  be  of  one  sign 
(qy  >  0;  see  Figure  4.1),  so  it  is  only  in  the  upper  part  of  the  water 
column  (T  >  9#C)  south  of  the  core  that  qy  <  0  somewhere  in  the 
Stream.  If  the  Gulf  Stream  can  be  modeled  as  a  quasi -geostrophic  zonal 
flow  confined  to  a  channel  (or  at  least  restricted  from  exchanging  energy 
with  its  environs),  then  linear  baroclinic  Instabilities  could  arise  in  the 


current.  So  many  qual iflcations  are  required  in  drawing  this  conclusion 
that  it  should  be  considered  a  guideline  to  further  studies  rather  than  a 
concrete  result. 

Inspection  of  a  potential  vorticity  section  constructed  from  the 
GUSTO  velocity  and  temperature  data  establishes  that  isostrophes  generally 
parallel  isotherms  north  of  y  =  10  km;  south  of  that,  there  are  two 
distinct  regimes  according  as  T  is  greater  or  less  than  about  10°C.  For 
T  <  10°C,  isostrophes  and  isotherms  remain  parallel.  For  T  >  10°C,  there 
are  strong  changes  in  potential  vorticity  on  isotherms,  including  reversals 
in  the  sign  of  q^  when  T  _>  12*C,  or  on  horizontal  surfaces  above  700  m. 
Relative  vorticity  contributions  to  q  are  as  large  as  25  percent,  and 
their  gradients  may  be  as  large  as  q  itself,  which  has  typical  values  an 
order  of  magnitude  greater  than  the  6-effect.  Finally,  there  is  a  very 
striking  change  from  low  to  high  values  of  q  on  isotherms  T  _>  10°C  near  y  = 
0,  which  appears  as  a  "wall"  of  potential  vorticity  when  q(T  =  constant)  is 
plotted  as  a  function  of  y,  as  in  Fig.  4.2. 

This  strong  "wall"  in  potential  vorticity  across  the  Gulf  Stream  has 
immediate  implications.  The  sharp  transition  from  low  to  high  values 
suggests  that  the  thermocline  Gulf  Stream  is  a  potential  vorticity  front 
and  can  be  modeled  as  an  interface  between  two  types  of  water  of  different 
potential  vorticity.  Pratt  and  Stern  (1985),  using  such  a  model,  have  been 
examining  the  time  growth  of  large  amplitude  meanders  on  the  interface, 
including  the  "wave-breaking"  of  the  meanders.  The  study  should  be  more 
relevant  to  the  development  of  Gulf  Stream  rings  and  meanders  than 
linearized  instability  models  (Talley,  1982,  e.g.),  because  it  includes 


large  amplitude  non-linear  effects. 

The  other  inmediate  Implication  is  for  the  motion  of  water  parcels  in 

the  Gulf  Stream.  In  the  absence  of  strong  forcing,  dissipative  or  mixing 

mechanisms,  water  parcels  tend  to  conserve  their  temperature  (density)  and 

potential  vorticity.  It  is  anticipated  that  over  long  downstream 

distances,  the  potential  vorticity  of  a  water  parcel  may  be  modified  due  to 

relatively  weak  processes,  though  quantifying  these  changes  is  a  difficult 

problem.  On  the  other  hand,  it  is  unlikely  that  a  water  parcel  on  the  12° 

isotherm,  for  example,  encounters  processes  sufficiently  strong  to  change 

-7»  -1  -1 

its  potential  vorticity  from  12  to  20  x  10  Cm  s  ,  the  change  in 

potential  vorticity  on  that  surface  over  just  50  km.  (Even  if  ez  had 

values  as  strong  as  .02°C/m,  as  it  does  only  at  575  dbar  in  the  core  of  the 

-7  « 

Stream,  a  change  in  q  of  3  x  10  C/m/s  would  correspond  to  increasing 

-4-1  o 

f  by  0.4  x  1(1  s  --in  other  words,  a  water  parcel  from  37  N  would 

have  to  go  to  about  63 °N'.)  Thus,  for  all  water  with  T  >  9°C,  there  is  a 

strong  constraint  against  cross-stream  excursions  of  water  parcels,  due  to 

the  potential  vorticity  wall  there.  Deeper  down,  where  isotherms  and 

isostrophes  tend  to  be  parallel,  this  constraint  is  relaxed,  suggesting 

that  water  parcels  should  have  little  trouble  in  crossing  the  Stream.  This 

conclusion  is  entirely  in  accord  with  what  has  been  observed  numerous  times 

in  data  from  S0FAR  floats:  namely,  that  floats  tracked  at  700  m  always 

seem  to  get  caught  by  the  Gulf  Stream  once  they  encounter  it,  while  tracks 

from  floats  at  2000  m  hardly  reveal  the  existence  of  the  Stream  above,  and 

apparently  have  no  difficulty  in  crossing  from  one  side  to  the  other. 

Although  they  are  isobaric  rather  than  isopycnal  floats,  the  same  reasons 


for  their  behavior  can  be  expected  to  apply,  since  (as  we  have  seen)  the 
potential  vorticity  wall  exists  on  horizontal  surfaces  as  well. 

Owens  (1984)  has  recently  discussed  the  tracks  of  17  such  floats 
deployed  at  700  and  2000  m  in  or  near  the  Gulf  Stream.  He  concludes  that 
the  difference  in  float  trajectories  could  be  attributed  to  the  relation 
between  particle  speeds  and  the  propagation  speed  of  meanders:  "At  shallow 
depths,  the  particles  are  advected  along  the  meandering  flow  field  while  at 
depth,  parcels  of  water  cannot  move  sufficiently  fast  enough  to  stay  with 
the  meander  patterns."  He  further  suggests  that  the  difference  between 
isobaric  and  isopycnal  surfaces  could  induce  an  "artificial  dispersion" 
beneath  the  thermocline  Gulf  Stream.  The  results  presented  here,  however, 
suggest  that  an  isopycnal  float  could  behave  in  much  the  same  way,  as  there 
is  no  apparent  dynamical  constraint  against  cross-stream  movement. 

4.4  Comparison  with  numerical  model  results 

Lastly  cross-sections  of  potential  vorticity  in  the  jet  of  Holland's 
numerical  model  3L-4  are  examined.  Because  the  density  is  constant  in  each 
layer,  a  cross-section  of  potential  vorticity  in  a  model  layer  corresponds 
to  that  on  an  isothermal  surface  In  the  data.  Also,  since  the  numerical 
model  has  been  constructed  In  a  quasi-geostrophic  framework,  just  the 
quasl-geostrophic  potential  vorticity  has  been  examined.  Temperature  does 
not  appear  explicitly  as  a  variable  In  the  numerical  model  formulation,  but 
rather  Implicitly  in  the  variation  of  layer  thicknesses  from  their  rest 
depths.  Then,  the  quasi-geostrophic  potential  vorticity  for  layer  k 
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+  (order  Rossby  number  terms),  since  jj-  =  0(e). 

Then,  because  the  are  constant,  is  more  simply  defined  as  in 
(3-4),  repeated  here  for  individual  layers: 


\  =  V  By  +V2V'1+  =  fQ(l  ♦  ^1)  +  ey  (4-la) 

q2  =  f0+  By  +V2V  h3-/2)  =  fQ(l  ♦  -5/^h-2)^y^2^  (4-lb) 


V  V  -  A- 


(4-lc) 


The  units  of  potential  vorticity  in  this  system  are  not  the  same  as  shown 
in  Fig.  4.1,  but  the  terms  are  analogous.  (Multiplying  (4-1)  by  some  value 
oz  yields  the  same  units.)  The  first  term  of  the  far  RHS  will  be 
referred  to  as  the  stretching  term,  and  corresponds  to  foz  in  the  form 
used  for  the  mooring  data. 

The  limited  vertical  resolution  of  the  numerical  model  and  the 
mooring  data  demands  that  care  be  taken  in  comparing  the  two.  The  top 
layer  of  the  model  reaches  to  only  300  meters’  depth,  while  the  mooring 
data  extends  no  farther  up  than  400  m.  However,  it  will  be  found  that  the 
structure  of  potential  vorticity  in  the  top  layer  is  comparable  to  that  at 
thermocllne  levels  at  the  mooring  site.  In  the  downstream  direction,  there 
is  little  qualitative  change  In  the  model  q  profiles  over  many  hundreds 
of  kilometers;  the  section  examined  has  been  chosen  for  computational 
reasons,  and  is  20  km  downstream  from  the  boundary  between  accelerating  and 


decelerating  jet  flow. 

Figure  4.4  (a  through  c)  shows  the  potential  vorticity  profile  for 

each  of  the  three  layers,  with  f  subtracted  out,  and  with  the  various 

components  plotted  as  well.  (Subtracting  out  fQ  would  be  comparable  to 

subtracting  out  a  value  of  fez  from  the  form  used  on  the  data,  where  §z 

at  a  level  was  a  spatial  average  over  the  whole  domain.)  First  consider 

the  potential  vorticity  cross-section  for  the  top  layer.  Figure  4.4a. 

Throughout  the  domain,  variations  in  q  are  dominated  by  the  stretching 

term,  while  8  makes  a  very  small  contribution,  and  relative  vorticity  is 

important  only  in  the  jet.  There  is  a  strong  gradient  of  q  in  the  jet, 

with  a  weak  minimum  and  secondary  maximum  to  the  south.  Although  these 

latter  features  are  present  in  the  stretching  term  alone,  they  are 

accentuated  by  the  contribution  of  relative  vorticity  to  the  profile.  The 

comparison  with  Figure  4.3,  showing  q  at  575  dbar,  is  quite  remarkable. 

Absolute  magnitudes  of  potential  vorticity  may  be  compared  by  dividing  q 

in  Figure  4.3  by  a  typical  value  of  ez  at  575  dbar  of  .015  °C/m,  then 

-4  -1 

subtracting  out  f  =  .89  x  10  s  .  Then,  for  example,  the  peak 

7® 

value  of  q  in  Fig.  4.3,  q  =  18  x  10  C/m/s,  becomes  q'  =  (q/ez)  - 
-5  -1 

fQ  =  3.1  x  10  s  ,  which  is  somewhat  less  than  the  maximum  values  of 
q  in  Fig.  4.4a.  In  both  profiles  the  stretching  term  is  dominant;  the 
ratio  of  relative  vorticity  to  stretching  is  about  25  percent  where 
attains  its  maximum  value  in  both  figures;  and  the  B-effect  is  an  order  of 
magnitude  smaller  than  the  stretching  terms.  The  qualitative  resemblance 
between  Figure  4.4a  and  one  of  the  Isothermal  profiles  of  q  shown  In 
Figure  4.2,  such  as  T  =  14#C,  is  equally  Important:  there  is  the  potential 
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Figure  4.4.  Values  of  quasi-geostrophic  potential  vorticity  Q  from  Holland 
numerical  model  3L-4.  Sections  are  720  km  east  of  western  boundary. 
Individual  components  are  indicated  directly  on  figure.  Horizontal  lines 
in  lower  half  mark  off  boundaries  between  energetic  regions  referred  to  in 
Chapter  3.  a)Layer  1;  b)layer  2;  cllayer  3. 


O 
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vortlclty  wall  at  or  near  the  jet  axis,  a  minimum  to  the  south,  and  a 
secondary  maximum  south  of  that.  The  potential  vortlclty  to  the  north  of 
the  wall  in  Fig.  4.4a  does  not  decrease  as  does  q  at  575  dbar,  but  again 
that  is  due  to  the  difference  in  evaluating  q  on  isothermal  and 
horizontal  surfaces. 

Notice  in  Figure  4.4a  that  the  jet  region,  as  characterized  by  the 
energetic  analysis,  corresponds  almost  exactly  to  the  region  over  which  q 
forms  the  strong  wall:  in  other  words,  the  gradient  of  q  does  not  change 
sign  in  the  jet  in  this  layer,  though  Figures  4.4b  and  c  show  very  weak 
reversals  in  qy  deeper  down.  On  the  basis  of  these  features,  it  might 
be  anticipated  that  baroclinic  instability  would  not  emerge  as  a  major 
energetic  component  in  the  jet  region;  indeed  in  Chapter  3  that  was  found 
to  be  the  case.  On  the  other  hand,  the  recirculation  region  of  layer  1 
shows  qy  <  0  for  the  most  part,  so  that  if  the  jet  and  recirculation 
areas  are  considered  together  qy  Indisputably  changes  sign  just  within 
the  top  layer.  In  Chapter  3  it  was  found  that  all  the  regions  are  coupled 
energetically,  particularly  these  two,  so  that  trying  to  apply  a  channel 
model  to  the  model  jet  is  unrealistic;  clearly  the  recirculation  region 
needs  to  be  considered,  especially  if  one  anticipates  the  existence  of 
baroclinic  instability.  Talley  (1982)  drew  a  similar  conclusion  from  a 
study  of  two-layer  jets  flanked  by  "westward  recirculations." 

Now  consider  Figures  4.4b  and  c  more  closely.  The  abscissa  scale  on 
these  two  figures  Is  different  than  that  of  Figure  4.4a,  as  Is  evident  by 
the  rather  dominant  contribution  of  sy  in  each  case.  The  middle  layer 
profile  shows  nicely  the  homogenization  of  potential  vortlclty  mentioned  in 
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Chapter  3,  and  how  it  arises  for  the  various  regions.  In  the  Sverdrup 
region,  stretching  and  b  offset  one  another;  this  same  balance  between 
layer  thickness  and  the  planetary  vortlclty  tendency  has  been  documented 
numerous  times  for  mid-depth  subtropical  gyre  flow  in  the  North  Atlantic 
(see,  e.g.,  Luyten  et  a1_. ,  1983;  or  McDowell  et  aK,  1982).  In  the  jet 
region  the  primary  balance  is  between  stretching  and  relative  vorticity, 
with  a  slight  net  increase  in  q  going  northward  through  the  jet.  In  the 
recirculation  region,  the  layer  thickness  is  nearly  uniform,  and  it  is  the 
relative  vorticity  offsetting  b  that  produces  uniform  q.  It  is 
difficult  to  make  close  comparisons  between  the  middle  layer  profile  and 
mooring  data,  except  to  note  that  in  the  latter  there  was  a  tendency  for 
q  to  be  constant  on  isotherms  for  T  <  8°C,  that  is,  below  the 
thermocl  ine. 

Finally,  Fig.  4.4c  shows  the  deep  layer  profile,  though  comparison 
with  the  mooring  data  Is  tenuous  at  best.  Potential  vorticity  variations 
are  dominated  by  sy  in  that  layer,  except  in  the  jet  region  where  the 
profile  is  fairly  flat.  Not  only  is  there  no  suggestion  of  a  potential 
vorticity  "wall"  in  this  layer,  but  q^  Is  actually  slightly  negative 
underneath  the  jet.  The  5#C  Isotherm  in  Figure  4.3  shows  that  while  q 
tends  to  be  constant  on  deep  Isotherms,  possibly  q^  is  slightly 
negative  beneath  the  Stream  axis  (and  positive  again  farther  north), 
according  to  the  mooring  data.  However,  the  similarity  between  these 
results  Is  probably  fortuitous,  for  the  values  of  q  at  1175  dbar  are 
dominated  by  changing  Tz,  rather  than  by  By  as  in  the  deep  layer  of 
the  numerical  model.  The  dominance  of  By  in  that  layer  is  a  signature  of 
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the  absence  of  thermohal Ine  forcing  In  the  model,  which  could  produce 
meridional  flow  in  the  deep  layer,  thus  requiring  isostrophes  to  depart 
significantly  from  latitude  circles.  The  structure  of  aT/az  (or  ap/az) 
in  the  ocean  --  prescribed  as  density  jumps  in  the  model  --  may  in  fact  be 
determined  to  some  extent  by  thermohaline  processes. 

4.5  Summary 

There  are  a  number  of  ubiquitous  features  which  appear  in  the 
observed  structure  of  the  Gulf  Stream  and  its  analog  in  theoretical  and 
numerical  models.  First  of  all,  there  is  a  tendency  towards  uniform 
potential  vorticity  in  the  shallower  part  of  the  anti -cyclonic  portion  of 
such  jets.  Second,  the  core  of  the  thermocline  jets  is  characterized  by  a 
strong  gradient  in  q,  or  a  potential  vorticity  wall,  along  isothermal 
surfaces.  In  a  way,  even  the  two-layer  inertial  jet  model  contains  this 
feature,  for  where  the  interface  of  the  top  layer  surfaces  there  is  a 
discontinuity  in  potential  vorticity.  Finally,  below  the  thermocline,  the 
observational  and  numerical  data  suggest  that  potential  vorticity  is 
uniform  on  isothermal  (isopycnal)  surfaces. 

The  implications  of  these  features  have  been  discussed  at  length,  but 
are  reiterated  here  for  emphasis.  A  uniform  potential  vorticity  model  is 
at  best  appropriate  for  describing  the  cross-stream  structure  of  the 
anti-cyclonic  part  of  the  jet.  Modeling  the  thermocline  Gulf  Stream  as  a 
potential  vorticity  front  is  better  for  Investigating  large  amplitude 
meandering,  while  linearized  Instability  models  have  dubious  application 
even  though  qy  changes  sign  across  the  Stream.  Finally,  the  difference 
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in  shallow  and  deep  potential  voi  ;ity  structure  may  be  largely 
responsible  for  the  difference  in  the  Lagrangian  flow  patterns  observed  in 
different  parts  of  the  water  column. 


Chapter  5.  Speculative  results  and  interpretations  of  mooring  data 

5.1  Introduction 

In  Chapters  2  and  4  it  was  demonstrated  that  by  using  an  inferred 
horizontal  coordinate,  it  was  possible  to  describe  the  average  velocity 
and  potential  vorticity  of  the  Stream,  or  calculate  fluxes  of  mass, 
momentum  and  kinetic  energy  for  four  separate  events  when  the  Gulf  Stream 
passed  across  the  mooring  site.  Moreover,  the  results  obtained  are  in 
good  agreement  with  past  results,  obtained  primarily  from  hydrographic, 
float,  and  current  profiler  data  all  the  way  across  the  Stream. 

In  this  chapter,  the  energetics  and  dynamics  at  the  mooring  site  are 
examined.  In  particular,  the  year-long  time  series  of  data  ought  to 
provide  new  insight  into  the  energetic  exchanges  between  mean  and  eddy 
flows,  since  for  the  first  time  there  are  concurrent  records  at 
thermocline  and  abyssal  depths.  It  has  proven  difficult  to  address  the 
dynamics  governing  this  complicated  flow,  but  a  kinematical  picture  may 
be  deduced.  Finally,  in  view  of  what  has  been  been  possible  with  the 
G'JSTO  data  set  and  what  questions  remain  unanswered,  suggestions  are 
presented  for  the  directions  future  programs  might  fruitfully  pursue. 

5.2  Energetics  at  the  mooring  site 

In  order  to  discuss  energy  exchanges  between  mean  and  eddy  flow,  it 
is  necessary  to  be  able  to  define  the  two.  Although  it  has  been  possible 
to  describe  an  average  Gulf  Stream  profile  with  horizontal  and  vertical 
structure.  It  is  unlikely  that  there  is  sufficient  data  to  discuss 
deviations  from  that  average.  Yet  if  the  strict  time  average  of  flow  at 


each  depth  is  examined,  the  horizontal  information  pertinent  to  a 
discussion  of  the  role  of  barotropic  instability  is  lost. 


Two  compromises  are  made  in  discussing  the  eddy  kinetic  and 
potential  energy  equations.  First,  they  are  considered  in  the  unrotated 
frame  rather  than  the  rotated  because  the  time-averaging  involved  in 
obtaining  the  equations  introduces  an  interpretative  problem;  in 
addition,  this  approach  is  more  traditional  and  is  more  readily  compared 
with  past  results  as  well  as  the  numerical  results  of  Chapter  3.  Second 
horizontal  resolution  is  limited  to  two  bins,  corresponding  to  < 

13°C  (north  of  the  Stream  axis)  or  >  13°C  (south  of  the  axis). 

—  — 2  — 2 

The  equation  for  eddy  kinetic  energy  K'  =  (l/2)p(u'  +  v'  )  has 

been  discussed  in  some  detail  already;  it  is  obtained  by  adding 
u'x(u-momentum  equation)  *  v'x(v-momentum  equation)  and  then 
time-averaging: 


(TE  +  0  lx  +  *  ly)ir= 


.  p  ( 


u'  +v' 


-)  - 


_  — 7 

u'V (vx+uy)  -  pu,cux 


-  pv  vy  -  .  VHp’  (5- 

where  the  subscript  H  is  used  to  mean  the  horizontal  components  only. 

_  i  — 2  - 

Similarly,  the  equation  for  eddy  potential  energy  r  =  -^  (ga0T*  )/( ©z) 

obtained  by  multiplying  equation  (2-15)  by  qa0T'/  "o  and  time- 
averaginq  is: 


(L-  +  u  i_  +  V  — VP"1" 

1TF  ax  ay' 


r.  a 


■qar 


-firr  tx  + 


Ty)  -  gaowT1-  - 


V..  .  (fc'T-7  gao/2eT) 


(5-2) 
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Recalling  that  -^.V^p'  =  -V^.p'fy  -  p'w^  =  -V^.p'^  +  w'p^  = 

-^.p'foj-  gp'W  =  -?2 -P'^3  +  ga0w'T‘,  (5-1)  and  (5-2)  can  be  added  to  get 
the  equation  for  the  total  eddy  energy: 

(jf  u  ij*  v  -|y)(lT>  -7H.5'  (K'+P* )  -  pu'v'fv/u  )  - 


— (u'T*  T  +  v'T'T  )  -  puv(u,z  -  v,Z)  -  7vp'Jll 


(5-3) 


The  total  eddy  energy  may  be  modified  by  essentially  two  types  of  terms 
that  appear  on  the  RHS  of  equation  (5-3):  1  Exchanges  between  the  mean 

and  eddy  flows,  represented  by  up-  or  down- gradient  momentum  and  heat 
fluxes;  and  2)"radiating"  terms,  which  appear  as  divergences  of 
quantities  depending  only  on  the  eddy  field.  From  (5-1)  and  (5-2)  it  is 
clear  that  the  exchange  between  eddy  kinetic  and  potential  energies  is 
given  by  the  term  +_  g<*  w'T1,  which  appears  with  opposite  sign  in  the 
two  equations. 

To  assess  the  relative  importance  of  terms  on  the  RHS  of  (5-3)  in  a 

gross  rather  than  localized  sense,  one  ought  to  Integrate  over  a  volume, 

as  in  Chapter  3.  Integration  in  z  is  possible  because  of  the  moorinq's 

vertical  resolution;  to  achieve  integration  in  y,  two  separate 

temperature  bins,  corresponding  to  regions  north  and  south  of  the  Stream 

axis,  have  been  used  for  the  time-averaging  process;  integration  in  x 

poses  some  difficulty  and  will  simply  be  ignored,  since  only  a  rough 

qualitative  picture  is  sought.  Furthermore,  it  is  useful  to  recognize 

that  even  in  unrotated  coordinates,  it  is  true  that  v  <<  u  ,  and 

— _  *  y 

— 2 

v'  <<  u'  ,  so  that  the  net  mean-to-eddy  momentum  exchange  is 
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approximately  given  by  -pu'v'u  -  pu'  u  .  To  calculate  the 

y  * 

Reynolds  stress  u'v',  all  the  data  were  divided  into  two  bins  according 
to  whether  T  >  13°C.  or  T  <  13°C  (where  uy  changes  sign).  Then,  for 
each  of  the  new  data  sets,  mean  and  eddy  velocities  were  computed,  as  was 
the  product  u'v'.  These  values  are  given  in  Table  5-1.  The  accompanying 
values  of  uy  were  obtained  by  taking  uy  at  575  dbar  equal  to 
e”1  times  the  maximum  value  attained  analytically  from  (2-B,9).  The 
magnitude  is  assumed  to  decay  with  depth  on  an  e- folding  scale  of  1000 
m.  The  width  of  the  anticyclonic  side  is  taken  as  50  km,  the  distance 
from  T57g  =  13 °C.  to  16.5°C.;  on  the  cyclonic  side,  Ay  r  80  km  (T575 


=  5°C .  to  13°C. ) , 

Previous  long  tine-series  from  measurements  in  the  deep  water  (4000 
m)  beneath  the  Gulf  Stream  (Schmitz,  1977)  have  suggested  that  u'v' 
changes  sign  across  the  axis  (geographical  average)  of  the  Stream  over 
perhaps  2-3  degrees  of  latitude,  such  that  there  is  a  flux  of  eddy  to 
mean  kinetic  energy  (u'v'  >0  south  of  the  Stream);  directly  under  the 
axis,  u'v'  <0  and  has  magnitudes  of  5-15  cm  s  .  The  Reynolds 
stresses  in  the  deep  part  of  the  water  column  here  are  not  terribly 
different  for  the  two  bins,  but  they  do  not  change  sign  across  the  Stream 
axis.  In  addition,  (u'v')y  >  0  across  the  Stream,  the  opposite  sense 
of  Schmitz's  findings.  However,  the  results  in  the  upper  1000  m,  which 
will  make  the  greatest  contribution  to  the  net  momentum  exchange  because 
Uy  is  strongest  there,  are  remarkably  different  from  all  the  deep  water 
values.  For  both  bins  (i.e.,  on  both  "sides"  of  the  Stream),  u'v1  is 
large  and  negative,  and  (u'v')y  >  0,  so  the  more  negative  values 
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occur  on  the  warm  side,  T  >  13*C,  where  also  uy  has  the  greater 
magnitude.  The  net  effect  is  that  of  a  down-gradient  eddy  momentum  flux 
(see  Table  5-1),  which  Implies  a  growth  of  eddy  energy  at  the  expense  of 
the  mean  kinetic  energy  via  this  mechanism.  Furthermore,  vrtille  the 
result  from  deep  measurements  cited  above  suggests  some  symmetry  of 
energetics  across  the  Stream,  the  shallow  measurements  show  there  Is  net 
southward  eddy  transport  of  eastward  momentum  across  the  Stream.  The 
mean  contribution  uv  Indicates  a  similarly  directed  flux,  an  order  of 
magnitude  smaller;  uV  <0  may  be  related  to  the  fact  that  the  average 
flow  was  south  of  east.  Fofonoff  and  Hall  (1983)  found  that  eastward 
momentum  flux  of  the  Gulf  Stream  Is  decreasing  In  this  region;  (u'v'ly 
>  0  is  one  mechanism  that  can  account  for  such  a  decrease  but  was  not 
calculable  in  that  work. 

The  effect  of  the  other  exchange  term  is  hard  to  determine  from  the 

mooring.  Some  historical  data  (Worthington,  1976;  or  Knauss,  1969) 

suggest  that  Gulf  Stream  transport  may  still  be  increasing  at  68°W,  and 
_  — V- 

if  ux  >  0  as  well,  then  -pu'  ux  would  offset  the  effect  of  the 

Reynolds  stresses  In  this  region.  Fofonoff  and  Hall  (1983)  found 

"T- 

ux  <  0  at  this  longitude,  however,  so  that  -pu'  ux  might  enhance 
the  Reynolds  stresses. 

Roth  similarities  and  differences  exist  between  the  data  and 
numerical  results.  On  the  average,  In  the  decelerating  portion  of  the 
numerical  jet,  the  Reynolds  stresses  are  positive  rather  than  negative  In 
the  upper  layer,  about  one- fifth  the  size  of  those  at  the  mooring  site  at 
575  dbar,  and  (u'v')y  >  0.  In  the  deep  layer  uV  changes  sign  across 


the  jet,  but  in  this  layer  (u'v')y  <  0.  However,  because  the 

horizontal  shear  is  so  much  stronger  in  the  top  layer,  there  is  a  net 

transfer  of  mean  to  eddy  kinetic  energy  once  the  budgets  are  integrated 

over  a  volume;  half  of  the  transfer  is  due  to  a  contribution  from  the 
*"7_ 

term  -pu'  ux.  To  compare  absolute  magnitudes,  suppose  that  the 
exchange  estimated  from  Table  5-1  occurs  over  a  width  of  130  km  (Tgy5  = 
5°-16 .5°  C)  and  downstream  for  1000  km.  Then  the  net  transfer  would  be: 

(777  x  50  km  -  252-^  x  80  km)  x  103km  -  19  x  109  J/s, 

g 

which  is  very  close  to  the  value  of  14.0  x  10  J/s  calculated  for  the 

same  mechanism  in  the  decelerating  jet  of  the  numerical  model. 

From  the  data  the  mean  to  eddy  potential  energy  flux  at  the  mooring, 

u'T1  T  +  v'T'  T  ,  can  be  estimated  without  the  necessity  of 
x  y 

separating  the  data  into  bins.  This  flux  is  down-gradient  as  well,  with 
an  integrated  magnitude  less  than  half  as  large  as  the  barotropic 
conversion,  dominated  by  values  from  the  upper  1500  m  or  so  of  the  water 
column  (see  Table  5-2).  As  expected,  v'T'  >  0,  but  all  the  correlation 
coefficients  for  the  heat  flux  calculations  are  small.  The  vertical  eddy 
heat  flux  term  gaQw'T 1  is  positive  (although  again  the  correlation 
coefficients  are  very  small),  and  may  be  estimated  from  values  at  575 
dbar  alone,  since  its  magnitude  falls  off  rapidly  to  negligible 
values.  At  575  db,  w'T'  =  24.979  x  10"3  #C  cm/s;  taking  az  =  725  m. 

Ay  =  130  km: 

fdy  (dz(  ga-WTTr)  =  (130  km)(725  m)(9.81  (10*4  -2L— )  x 

J  J  0  s^  cnr  C 


( 24.979xlO”3*C  cm/s) 


2.310  x  10*  kg  m/s*. 


<db)  (db  >  m)  (cm/s/m)  (10  5°C/m)  (°C  cm/s)  (10'  3 cm/s/m)  (10  eC/m)  (°C  cm/s)  (°C/m) 
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Since  this  term  appears  with  a  minus  sign  In  (5-2),  eddy  potential  energy 
evidently  Is  being  converted  into  eddy  kinetic  energy;  moreover,  It  is 
much  larger  than  the  release  of  mean  to  eddy  potential  energy,  so  that 
the  net  tendency  of  eddy  potential  energy  Is  to  decrease.  The  energy 
pathway--  mean  potential  to  eddy  potential  to  eddy  kinetic  energy-- 
tantalizingly  suggests  the  presence  of  baroclinic  Instability  at  the 
mooring  site.  Caution  is  warranted,  however,  for  all  of  the  calculations 
involved  are  rather  noisy:  the  correlation  coefficients  are  generally 
small  and  not  significantly  different  from  zero  for  the  year  long  data 
record. 

The  net  results  are  not  unlike  the  numerical  model  findings  in  the 
decelerating  jet  and  associated  recirculation  region:  both  types  of 
instabilities  appear  to  be  present,  with  barotropic  instability  possibly 
dominating  within  the  Stream.  The  difference  in  sign  of  the  Reynolds 
stresses  in  the  two  cases  (observational  vs.  numerical)  is  curious  and 
may  be  related  to  the  symmetry  of  the  double  gyre  model.  The  energy 
transfer  of  mean  to  eddy  fields  implies  local  growth  of  eddy  energy,  mean 
or  eddy  advection  of  eddy  energy  away  from  the  mooring  site,  or  radiation 
of  eddy  energy  away  from  the  mooring  site.  The  term  -( a/ay) (v ’ (K'+P* )) 
may  be  estimated  with  the  "two-bin  method"  used  on  the  Reynolds  stresses, 
and  although  it  has  the  desired  sign  to  balance  the  momentum  and  heat 
fluxes,  it  is  at  least  an  order  of  magnitude  smaller  than  those  terms. 
Moreover,  (v'K')y  alone  has  the  opposite  effect  as  in  the  model,  though 
it  is  relatively  much  smaller.  All  of  the  remaining  radiation  type  terms 
are  completely  intractable.  It  is  plausible  that  this  region  is  one  of 


either  temporal  or  spatial  eddy  growth,  and  corresponding  time  and  space 
scales  for  the  Implied  growth  may  be  estimated  from  the  numbers  so  far 
derived,  along  with  an  estimate  of  net  eddy  energy  from  values  In  Table 
2-1.  The  average  along-stream  velocity  u  Is  an  Integrated  value  over 
the  depth.  Then: 


Jdyjdz  j  (p(u’2  ♦  v'2)  +  ^ao^  )  =  8.06  x  lO^kg  m/s2  ; 


-jdy^dz(p  u'V  uy)  -  1 .87  x  10*  kg  m/s3  ; 


Jdy^dlz( ^o) ( vTTT~ Ty  ♦  uTTx)  =  7.36  x  in3  kg  m/s3  ; 


u  =  9  cm/s 

The  above  values  Imply  either  a  growth  rate  r  (calculated  from  |{)^versT6n^ 


_  .  1.87  x  104  ♦  7.36  x  103  .-1  „  00  ..  ,„-7_-l  1  „ 

r  - r* - s  =  3.23  x  10  s  ==>  —  =  36  days 


8.06  x  10J 


or  a  downstream  scale  Lx  for  eddy  energy  growth  (calculated  from 


L  -  of 

x  conversion 


(9  cm/s)(8.06  x  1010  kg  m/s2) 
(1.87  x  104  ♦  7.36  x  103)kg  m/s 3 


=  278  km. 


The  time-averaging  process  over  a  360-day  period  automatically  obscures 
the  Implied  36-day  growth  rate.  Although  the  predicted  scale  for 
downstream  growth  seems  reasonable,  it  should  be  noted  that  the  mooring 
site  Is  In  a  region  that  has  been  Identified  as  a  maximum  In  eddy  kinetic 
and  potential  energies  (Schmitz,  1984;  Richardson,  1983). 
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The  analysis  of  Chapter  3  showed  that  unlimited  spatial  and  temporal 
coverage  of  the  free  jet  does  not  resolve  all  possible  questions 
concerning  Its  energetics.  On  the  other  hand,  diagnosis  of  numerical 
models  needs  to  be  pursued,  since  the  streamfunctlon  and  hence  the 
pressure  work  terms  are  known  everywhere.  However,  the  lack  of  parity 
between  the  accumulation  of  data  and  yield  of  results  suggests  that  a  few 
well  placed  moorings  might  be  as  valuable  as  a  large,  dense  array  for 
addressing  questions  of  the  type  discussed  here. 

5.3  A  kinematic  framework  for  Interpreting  the  flo¥ 

Because  the  equations  for  mean  and  eddy  energies  do  not  depend  on 
small  amplitude  expansions,  different  terms  in  them  may  be  examined  for 
signatures  of  familiar  processes  such  as  barotropic  or  baroclinic 
instabilities.  Analysis  of  dynamical  balances  is  less  tractable, 
however,  because  the  relevant  terms  in  the  vorticity  equation,  for 
example,  involve  so  many  derivatives.  Moreover,  with  Rossby  numbers  of 
about  0.3,  it  may  be  necessary  to  search  for  a  new  dynamical  framework 
that  explains  the  flow.  Although  such  a  dynamical  framework  has  not  been 
fully  developed,  a  kinematic  framework  has  been  explored  that  is 
consistent  with  the  data.  It  is  just  one  interpretation  of  what  is 
occurring  at  the  mooring  site,  and  is  not  necessarily  unique. 

Once  the  data  from  the  mooring  had  been  scrutinized,  a  nurtber  of 
peculiarities  emerged.  The  barotropiclty  of  the  cross-stream  velocity 
field  was  not  anticipated,  nor  were  the  large  magnitudes  of  the  vertical 
velocities.  The  latter  feature  usually  resulted  from  the  difference 


between  the  apparent  translational  velocity  of  the  Stream  (as  Indicated 
by  aT/3t)  and  the  measured  cross- stream  velocities.  Additionally,  there 
Is  the  lack  of  symmetry  In  details  of  the  four  Individual  Gulf  Stream 
crossings.  An  Interpretation  Is  sought  to  relate  these  unusual  features. 

The  scenario  Is  based  on  the  supposition  that  primarily  the  lower 
layer  of  a  two- layer  system  is  observed,  as  suggested  by  the  vertical 
structure  of  vertical  velocity  In  Chapter  2.  In  the  discussion  of  the 
vortlcity  equation.  It  was  found  that  the  vertical  and  cross-stream 
velocities  were  better  correlated  than  v  and  wz;  that  was  ascribed 
to  the  linear  (two- layer)  structure  of  w  and  the  fact  that  wz  would  be 
a  noisy  time-series.  It  was  also  pointed  out  that  the  correlation  can  be 
explained  by  the  following  argument:  vertical  velocities  are  Induced  at 
the  bottom  by  flow  up  or  down  the  bottom  slope;  higher  In  the  water 
column,  water  parcels  moving  vertically  must  also  move  horizontally  -- 
primarily  cross-stream  --  to  remain  on  Isotherms.  Given  that 
cross-stream  velocities  are  nearly  barotroplc,  as  has  been  shown 
empirically,  then  the  vertical  velocities  will  be  greatest  vrfiere  the 
Isotherm  slopes  are  greatest,  that  Is,  in  the  thermocllne  or  at  Its 
equivalent,  the  Interface  in  the  two- layer  model.  The  orientation  of 
Isotherm  slopes  may  be  quite  different  from  the  direction  of  the  bottom 
slope,  whence  the  bottom  vertical  velocity  appears  as  an  Independent 
forcing  mechanism. 

So  far,  of  course,  all  that  has  really  been  said  Is  that  w  * 
u'VZj  =  v  sz^/sy,  where  Zy  Is  the  depth  of  an  Isotherm,  and 
since  v  Is  barotroplc  sz^/sx  Is  assumed  to  be  zero.  The  local 


temperature  change  is  important  too,  however,  and  may  be  represented  as 


Tt  =  -VT*.  where  V  is  then  related  to  the  cross-stream  translational 
velocity,  so  that  w  =  (v-V)3Zj./ay.  It  has  been  noted  that  there  are 
events  for  which  the  apparent  translational  and  cross-stream  velocities 
are  opposite  in  sign;  then  |v  -  vl  >  j v /  and  w  is  larger  in  magnitude. 

On  the  one  hand,  the  EOF  cross- stream  velocity  amplitude  is  highly 
correlated  with  wg  (vertical  velocity  at  4000  m) ,  but  not  at  all  with 
3T/st  at  575  dbar.  One  way  to  interpret  these  results  is  to  assume  that 
the  barotropic  cross-stream  velocity  arises  in  response  to  wg.  The 
size  of  the  vertical  velocities  throughout  the  remainder  of  the  water 
column  depends  on  the  Gulf  Stream's  translation:  in  the  case  of  opposing 
translational  and  measured  velocities,  then  wz  below  the  thermocline 
should  have  the  same  sign  as  wg: 


17  =  vHjfr(S75,  -Iir(4000));  (*.-*).,  >  o;  Iff 

ay  3  v  3  v 


575 


ay 


sy 


4000 


Thus,  if  wg  >0  (whence  v  >  0)  then  T^  <  0  should  imply  wz  >  0;  while 
wg  <  0  (v  <  0)  and  T^  >  0  implies  wz  <  0.  In  fact,  aT/at  at  575  dbar 


and  aw  =  w. 


575 


-  w, 


'4000 


are  negatively  correlated  with  C  =  -.75,  which  is 


significant  at  the  95  percent  confidence  level.  Since  that  result  comes 
from  using  the  entire  time  series,  and  since  vT^  <  0  only  about  half 
the  time,  there  are  evidently  two  different  flow  regimes,  that  can  be 
summarized  as  follows: 


Figure  5.1. 
system  for 
the  bottom, 
ocean, 

1 nterface. 
negative,  1 


Vertical  velocity  amplitude  structures  for  a  two- layer 
cases  discussed  in  text.  Slight  bottom  slope  allows  w  M  at 
Interface  (dashed  line)  represents  thermocline  of  real 
=  vertical  velocity  at  bottom;  Wj  =  vertical  velocity  at 
Case  la  or  II a  occurs  when  given  amplitude  structure  is 
.e. ,  Wg  <  0 . 


Case  I 


Case  I  (Case  la) 


Case  II  (Case  I la ) 


wg  >  0  (Wg  <  0)  >  0  (wB  <  0) 

v  >  0  (v  <  0)  v  >  0  (v  <  0) 

Tt  <  0  (Tt  >  0)  Tt  >  0  (Tt  <  0) 

wz  >0  (wz  <  0)  wz  <0  (wz  >  0) 

Figure  5.1  shows  vrfiat  the  vertical  velocity  structure  would  be  In  the 
two-layer  system  for  these  two  cases;  In  Case  II,  there  are  two 
possibilities,  since  w  at  the  interface  need  not  have  the  same  sign  as  w 
at  the  bottom;  but  the  data  suggest  the  baroclinic  response  (shown  on  the 
right  in  Fig.  5.1  under  Case  II)  is  more  typical.  Johns  and  Watts  (1985) 
present  a  linear  analysis  of  the  temperature  equation,  from  data  just 
downstream  of  Cape  Hatteras,  which  yields  results  analogous  to  Case  I 
described  here;  but  in  that  study,  Case  I  evidently  described  most  of  the 
data,  and  Case  II  was  not  considered  at  all.  Their  Fig.  11  depicts  Case 
I  very  nicely. 

The  individual  events  are  described  rather  well  by  the  various 
cases,  and  all  cases  except  1 1  a  occur  (in  other  words,  Case  II  does  not 
occur  with  negative  bottom  velocities).  The  March  and  early  September 
crossings  are  examples  of  Case  I:  temperatures  are  decreasing,  but 
v  >  0;  examination  of  the  vertical  velocity  time  series  shows  that 
wB  >  0  and  wz  >  0  (vrtiere  wz  is  taken  between  the  thermocline  and  the 
bottom).  The  June  event  is  a  combination  of  two  cases.  Clearly 
Tt  <  0  the  whole  time,  but  v  <  0  for  the  first  11  days  and  v  >  0 
for  the  remainder  of  the  crossing.  Accordingly,  Wg  changes  sign  from 
negative  to  positive  after  10  days;  and  wz  is  generally  less  than  zero 
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throughout  (this  statement  will  be  qualified  below).  Finally,  the  late 
September  event  falls  under  Case  II  with  >  0. 

That  is  just  the  beginning  of  the  story,  however.  In  the  first 
place,  3W/3Z  has  Implications  for  the  vorticlty  balance.  Moreover, 
since  e  for  the  flow  has  been  estimated  to  be  as  large  as  0.3,  3w/az 
may  be  important  to  the  mass  balance.  Finally,  the  above  cases  are 
really  only  a  one-dimensional  description  of  the  flow:  what  is  the 
associated  three-dimensional  picture  Recall  that  in  Section  2.5.2,  it 
was  suggested  that  curvature  of  the  Stream  could  be  important  to  the 
vorticlty  balance.  To  get  at  the  curvature  directly,  consider  vorticity 
in  cylindrical  (rather  than  rotated)  coordinates.  Figure  5.2  gives  a 
definition  sketch  for  the  variables.  Then 


>  1  1  (vrl  I—  1  +  11  1  lii 

**  =  7  sr'  '  ~  r  3X  =  r  3r  *  r  3x 
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(Note  that  now  v  is  long-stream  velocity,  and  may  be  negative  or  positive 

3  1  3 

according  to  the  curvature  of  the  Stream.)  With  v  >>  u,  jp  >>  , 

CV  3  V 

~  -  +  —  and  the  vorticity  equation  becomes 

i  or 

2  2 

3  /v  *  +  u  3V  UV  +  „3V  .V  3V  .  V  3  V  . 

IFF  3r)  r3F‘7  u^  7'ar  rsxsr 


b(u  sin  x 


,  ,  3W 

+  V  cos  x)  =  f  — ■ 


(5-4). 


Now  the  local  change  of  curvature  explicitly  appears,  and  its  size  can  be 
estimated  by  referring  to  maps  constructed  from  satellite  data.  Figure 
5.3  shows  schematically  how  the  curvature  changes  from  May  27  when  r  is 
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Figure  5,3.  Schematic  showing  change  of  curvature  in  Gulf  Stream  from 
May  27,  1983  to  June  1,  1983.  Cross  indicates  mooring  site. 
Path  is  adapted  from  northern  edge  of  front  as  shown  on 
satellite  composites.  Dotted  circle  has  radius  of  about  70  km 
and  approximately  matches  curvature  of  Stream  at  mooring  site 
on  May  2  7  . 
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roughly  70  km,  to  June  1,  when  the  flow  has  straightened  out  so  r  i  s 
essentially  infinite.  During  this  time,  the  along-stream  velocity  at  575 
dbar  v  ~  30  cm/s  whence 

3  /  v  x  „  0  -  (30  cm/s)/(  70  km)  on  „  m-11  --2 

a Vr]  - 5  days -  =  x  10  s  • 

Meanwhile,  |y  between  the  bottom  and  thermocline  is  negative,  and  has 
an  estimated  magnitude: 

f  ||  ~  (.89  x  10-4  s-1)  x  ^42v1'\-3--e— ■■  =  -.78  x  10'Us_2 

Thus,  the  effect  of  changing  curvature  is  more  than  enough  to  balance  the 
squashing  in  the  lower  part  of  the  water  column.  (Notice  that  the  same 
balance  cannot  obtain  above  the  thermocline,  where  w2  must  be  <  0.) 

Proceeding  In  a  similar  but  qualitative  manner  for  the  four 
individual  events  suggests  that  the  observed  flow  patterns  can  be 
accounted  for  by  quasi-fixed  spatial  patterns  like  meanders  being 
advected  past  the  mooring  site,  or  by  patterns  propagating  past  the 
site.  There  is  qualitative  agreement  between  the  calculated  long-stream 
direction  of  flow  and  the  apparent  direction  from  the  satellite  pictures, 
indicating  that  surface  patterns  broadly  reflect  structure  in  the  deeper 
flow  (Fofonoff,  personal  communication).  Figure  5.4  shows  how  this  idea 
of  moving  patterns  Is  consistent  with  all  the  calculations  from  the  data 
for  those  events.  Single  line  arrows  are  selected  dally  long-stream 
directions,  vrfiich  when  placed  end  to  end  suggest  a  spatial  pattern  vrtiich 
could  account  for  flow  direction  at  the  mooring  site  If  the  feature 
passes  over  the  mooring  site  in  the  general  direction  shown  by  the  double 
dashed  line  arrows.  The  X‘s  show  successive  positions  of  the  mooring 


Figure  5.4.  Single  line  arrows  point  In  direction  of  flow  for  selected 
successive  dates  during  each  Gulf  Stream  passage  event;  length  Is 
proportional  to  time  between  successive  arrows.  Double  dashed  lines  are 
velocities  of  meanders  with  shapes  outlined  by  single  arrows,  required  to 
account  for  flow  at  mooring  site.  Successive  qualitative  positions  of 
site  are  Indicated  by  X's.  Along-stream  velocity  In  cylindrical 
coordinates  Is  Indicated  by  v*.  Relevant  Information  on  each  event 
according  to  classification  scheme  discussed  in  text  Is  listed  with  each 
feature. 
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relative  to  the  propagating  features.  With  each  feature  is  a  summary  of 
the  behavior  of  relevant  quantities  during  the  event,  and  its 
classification  according  to  the  above  cases.  The  June  event,  vtoich 
changes  character  halfway  through,  has  already  been  considered  in 
detail.  It  was  noted  that  in  late  May/early  June,  bottom  vertical 
velocities  were  negative,  and  accordingly  cross-stream  velocities  were 
negative  as  well.  However,  temperature  was  locally  increasing  so  that 
evidently  Case  la  is  occurring.  Consistent  with  this  conclusion, 
aw/az  <  0  during  that  time  frame,  and  (as  calculated  above) 

(  yr  )  <  0  as  well.  Between  June  5  and  7,  wg  and  v  change  sign  and 
the  flow  straightens  out  to  a  steady  direction  of  about  90°  true,  while 
aT/at  remains  positive;  this  case  is  like  II  if  aw/az  is  estimated  from 
575  to  4000  dbar,  which  yields  negative  or  small  positive  values. 

However,  aw/az  between  875  dbar  and  the  bottom  is  definitely  positive  for 
the  remainder  of  the  event. 

March  and  early  September  are  good  examples  of  Case  I,  with  the  sign 
of  a/at  (v/r)  consistent  with  the  overall  stretching  between  thermocline 
and  bottom  during  those  events.  As  the  meanders  propagate  or  are 
advected  past  the  mooring  site,  there  is  a  shift  in  each  case  from 
anti-cyclonic  to  cyclonic  flow.  In  later  September,  Case  II  is  observed, 
although  a  clear  indication  of  the  curvature  tendency  for  this  event  is 
lacking.  Going  in  detail  through  the  data,  one  can  find  other  isolated 
examples  of  shorter  duration  that  are  also  consistent  with  the  schematic 
interpretation  presented  and  fall  into  one  of  the  four  cases  enumerated 
above. 


A  distinction  has  been  made  above  between  meanders  "propagating"  and 
"being  advected"  past  the  mooring.  There  is  a  conceptual  difference 
between  the  two  possibilities,  the  latter  occurring  as  a  result  of  an 
ambient  velocity  field  outside  the  Stream  "pushing"  the  meander  along. 

In  this  case,  the  transport  estimates  made  from  the  current  meter  data 
can  be  altered  substantially  by  the  existence  of  the  ambient  field, 
because  without  explicit  knowledge  of  what  the  ambient  velocity  is,  it  is 
impossible  to  separate  it  from  the  velocity  structure  of  the  Stream 
itself.  The  projection  of  the  anbient  velocity  onto  the  along-stream 
direction  can  then  augment  the  transport  estimate  for  an  event  as 
follows:  suppose  the  ambient  field  is  zonal,  while  the  Stream  is 
directed  at  an  angle  a  to  due  East.  Then  the  apparent  transport  will 
exceed  the  actual  Gulf  Stream  transport  (which  would  be  measured  if  one 
could  move  with  the  meander)  by  an  amount 


&T  uan^  COSa  Ay  AZ 

which  for  uamb  =  5  cm/s,  a  =  45°,  Ay  =  100  km,  and  az  =  4000  m  is: 
AT  =  5  x  10~2  ™  x  x  105m  x  4  x  103m  =  14  x  106m3/s , 


which  is  comparable  to  the  transport  differences  between  the  various  Gulf 
Stream  events. 


Finally,  consider  the  importance  of  stretching  to  the  mass  balance. 

In  ordinary  quasi -geos trophic  dynamics,  to  lowest  order 

u  +  v  =  0.  At  the  GUSTO  site,  however,  it  is  possible  that  w 
x  y  z 

affects  the  mass  balance  at  lowest  order.  To  test  this  idea 


quantitatively,  continuity  is  integrated  over  a  cross-section  normal  to 
the  Stream,  from  575  to  4000  dbar: 


151 


C  I5”  77  **  dz  *  f75<JN-;s>  il  ‘  f "  <w575-  w4000>  * 
ys  ‘'4000  J4000  ys 


=  0 


(5-5) 

If  It  Is  assumed  that  the  Stream  maintains  a  ~teady  width  (an  implicit 
assumption  throughout  the  analysis  thus  far),  then  v^  =  v$  and: 

77  =  (>,575-  w4000)dy-  M  -  4000  “  dy  d2  (5'6) 

Thus,  the  transport  calculated  for  the  "lower  layer"  can  change  in  the 
downstream  direction  if  there  is  squashing  or  stretching  in  that  part  of 
the  water  column.  The  RHS  of  (5-6)  has  been  estimated  for  the  March  and 

June  events.  For  March,  when  wz  was  basically  positive,  the  RHS  has  a 

2  2 
value  of  -59.7  m  /s.  For  June,  wz  <  0  and  RHS  *  49.8  m  /s.  (Notice 

that  the  widths  of  the  two  events  from  Table  2-3  are  nearly  the  same.) 

c  o 

Estimated  transport  for  the  two  events  differs  by  32  x  10  m  /s , 

6  3 

about  half  of  which  occurs  below  575  db.  With  aM  =  16  x  10  m  /$,  a 
downstream  distance  ax  can  be  estimated,  over  which  squashing  or 
stretching  must  act  to  produce  the  observed  transport  difference: 


AX 


16  x  10^  m^/s 
(49.8  to  59.7)  m/s 


268  to  321  km. 


Over  a  length  scale  of  about  300  km,  a  change  in  transport  below  575  db 
can  occur  that  Is  comparable  to  the  observed  differences  between  the 
March  and  June  events.  This  length  scale  is  intriguingly  similar  to  the 
downstream  spatial  growth  scale  estimated  in  5.2.  Since  velocities  have 
merely  been  extrapolated  to  the  surface  to  obtain  the  total  transports, 
they  reflect  the  changes  observed  below  575  dbar.  However,  in  the 


situation  described  above,  in  the  "upper  layer"  wz  generally  should 
have  the  opposite  sign  as  in  the  lower  layer,  and  a  compensating  change 
in  transport  ought  to  be  observed  in  that  layer  if  instruments  were  there 
to  measure  it. 

Alternatively,  the  assumption  of  a  fixed-width  Stream  may  sometimes 
be  violated.  Then,  balancing  the  last  two  terms  on  the  LHS  of  (5-5) 
would  give  (assuming  v^,  vs  are  approximately  barotropic): 

(vM  -  vs)  ~  0(50  m2/s)  ==>  vM  -  v$  ~  0(1.5  cm/s), 
so  that  in  the  presence  of  stretching  (squashing),  the  Stream  would  be 
narrowing  (widening)  at  at  rate  of  1.5  cm/s  or  about  1.5  km/day. 

Clearly,  with  a  single  mooring  that  depends  on  an  assumption  of  fixed 
width  for  horizontal  information,  such  a  possibility  cannot  be  tested. 

Although  the  scheme  presented  here  can  acccount  for  some  of  the 
peculiarities  observed  in  the  data,  it  is  incomplete  in  the  sense  that  it 
is  not  predictive.  Even  if  the  behavior  of  the  bottom  vertical  velocity 
is  known,  for  example,  it  is  unknown  whether  a  Case  I  or  II  type  of  event 
is  occurring.  More  serious,  perhaps,  is  not  knowing  what  gives  rise  to 
wB  in  the  first  place:  is  it  indeed  independent  of  the  response  in  the 
rest  of  the  water  column,  or  is  it  somehow  an  integral  part  of  that 
response  Closer  examination  of  the  velocity  time-series  at  4000  db 
suggests  that  the  flow  there  is  actually  more  complicated  than  the  cross- 
and  long-stream  EOF  decomposition  for  the  water  column  would  indicate: 
at  times  the  cross-stream  velocities  there  are  bigger  than  the 
long-stream  component,  and  up  to  twice  as  large  as  cross-stream 
velocities  in  the  remainder  of  the  water  column.  Since  GUSTO  is  the 


first  example  of  concurrent  deep  and  thermocllne  Gulf  Stream  velocity 
time-series.  It  has  raised  many  Interesting  questions 

5.4.  Indication  for  future  directions 

Now  that  it  Is  technically  feasible  to  collect  long  time-series  of 
current  and  temperature  measurements  throughout  the  water  column  In  the 
Gulf  Stream,  our  understanding  of  that  current  should  advance 
dramatically.  The  GUSTO  mooring  was  the  first  successful  deployment  of 
such  a  mooring,  and  the  results  from  the  data  collected  point  to  the 
directions  that  future  Investigations  might  take.  Analysis  of  the  GUSTO 
data  appears  to  justify  the  identification  of  the  Gulf  Stream  primarily 
as  a  discrete  feature  with  a  well-defined  velocity  structure.  Thus,  a 
few  well-placed  moorings  can  provide  considerable  coverage  of  the  Stream, 
if  use  is  made  of  temperature  as  a  horizontal  coordinate,  and  If  a 
reasonable  definition  of  flow  direction  is  applied.  Using  the 
decomposition  of  velocities  into  their  along-  and  cross-stream  components 
shows  that  the  vertical  structure  at  the  GUSTO  mooring  site  Is  accounted 
for  by  a  baroclinlc  along-stream  model  and  a  decoupled,  barotroplc 
cross-stream  mode.  The  Inferred  horizontal  Information  may  be  used  to 
estimate  mass,  momentum,  and  kinetic  energy  transports  of  the  Gulf 
Stream.  It  Is  also  possible  to  construct  a  horizontal  and  vertical 
profile  of  the  average  Stream  velocity  structure,  from  which  an  average 
potential  vortlclty  section  may  be  constructed.  The  current  meter  data 
have  also  been  used  to  deduce  the  existence  of  strong  vertical  velocities 
In  the  Stream,  with  maximum  rms  values  at  thermocllne  levels,  and  a 
vertical  structure  resembling  the  first  baroclinlc  mode. 
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Certain  types  of  Information  are  required  to  clear  up  points  that 
are  left  ambiguous  by  the  GUSTO  data  set.  Data  from  the  thermocllne  to 
the  surface  are  needed  to  determine  vrtiether  the  velocity  structure  Is 
sharper  on  the  cyclonic  or  anti-cyclonic  side.  Also,  the  Rossby  nuiriber 
has  been  estimated  to  be  0.2  or  0.3  at  thermocllne  levels,  but  could 
possibly  be  larger  at  shallower  levels  where  the  relative  vortlcity  -uy 
may  be  quite  strong.  Finally,  results  from  the  GUSTO  data  have  suggested 
that  a  two- layer  system  might  adequately  model  Stream  dynamics; 
information  on  the  velocity  structure  above  the  thermocline  is  necessary 
to  determine  whether  or  not  two  layers  would  be  sufficient. 

Yorticlty  balances  were  difficult  to  determine  at  the  mooring  site, 
but  the  vertical  velocity  calculations  suggested  that  stretching  was  very 
important.  Rudimentary  comparison  with  changes  In  the  Stream  path 
curvature  showed  that  they  were  probably  sufficient  to  balance  the 
stretching,  but  careful  and  detailed  Inspection  of  satellite  data  in 
conjunction  with  current  meter  data  will  be  required  to  answer  this 
question  more  fully.  The  multiple  horizontal  derivatives  Involved  in  the 
vortlcity  equation  require  greater  spatial  resolution  than  a  single 
mooring  can  provide.  More  than  one  mooring  deployed  in  the  cross-stream 
direction  would  be  useful  In  addressing  a  number  of  other  Issues  as 
well.  It  was  found  in  Section  5.3  that  the  energy  budgets  at  the  mooring 
site  could  not  be  very  accurately  determined;  the  analysis  of  barotroplc 
energy  exchanges  between  eddy  and  mean  flow  especially  would  have 
benefited  from  greater  cross-stream  resolution.  Transports  could  be 
monitored  more  accurately  with  more  moorings  across  the  Stream,  since  It 
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would  take  less  time  for  the  entire  Stream  to  be  profiled  In  a  single 
passage  over  the  sites;  then  It  might  be  possible  to  determine  whether 
transport  really  changes  as  dramatically  In  time  as  the  GUSTO  data 
suggest,  or  vrfiether  apparent  changes  are  due  to  an  ambient  velocity  field 
advectlng  features  past  the  mooring  site.  It  would  also  be  possible  to 
test  the  extent  to  which  the  Stream  maintains  a  constant  width. 

The  most  intractable  aspect  of  observational  and  theoretical  Gulf 
Stream  analyses  Is  that  of  downstream  changes  In  the  flow  structure, 
because  they  are  so  slight  that  their  signal  Is  swamped  by  variations  in 
time  and  In  the  cross-stream  direction.  The  diagnostic  energetic 
analysis  of  Chapter  3  yielded  several  results  that  might  be  tested 
observatlonally,  however.  For  example,  in  the  jet  itself,  the  terms 
Involved  In  maintaining  the  mean  kinetic  energy  budget,  particularly  In 
the  decelerating  portion  of  the  jet,  were  greater  by  an  order  of 
magnitude  than  those  In  the  eddy  kinetic  energy  budgets.  Analysis  of 
data  from  moorings  separated  by  as  much  as  1500  km  In  the  downstream 
direction  might  be  able  to  Identify  whether  that  Is  Indeed  true. 

(Fofonoff  and  Hall  (1983)  tried  to  address  the  point  but  had  to  take 
relatively  synoptic  data  as  representative  of  the  mean,  so  that  eddy 
energies  had  to  be  Ignored.) 

One  of  the  major  thrusts  of  this  work  has  been  to  justify  the 
treatment  of  the  Gulf  Stream  as  a  well-defined  flow  structure  In  the 
ocean.  This  approach  gives  rise  to  the  most  basic  philosophical  question 
that  Investigators  must  address  In  the  future:  when(1f  ever)  Is  It 
appropriate  to  use  Eulerlan  averages  In  examining  regions  of  Gulf  Stream 


flow,  as  was  done  In  the  numerical  analysis  Our  definitions  of  the 
general  or  time-averaged  ocean  circulation  may  have  to  be  refined,  that 
we  may  distinguish  between  an  observable,  Eulerlan  average  crlculatlon  In 
the  ocean,  and  the  existence  of  a  boundary  current  with  an  average 
structure  but  variable  position  and  orientation,  such  that  It  affects  the 
interior  general  circulation  in  order  one  fashion.  It  seems  that 
observational  tools  are  sufficiently  advanced  to  resolve  these 
philosophically  different  approaches  In  the  near  future. 
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Appendix  A.l 


It  is  to  be  shown  that 


fx<“kV2V  *  Ty<\  \  »2  V  ’ 

vrfiere  k  Is  the  layer  index  and  KE^  =  ^  (V^2  ♦  For  neatness,  in  the 

following  the  subscript  k  Is  omitted.  Thus: 

i_(u  <hH)  *  |p.  (v  H2Y-)  = 

>  Ux>  *  V»  *  V  * u  **«w  *  lu1f>x  vyy  * 
v  +  fxxy  ♦  (v  f  >y  V-xx  . 

fc"  t  txx)  ♦  |y(v  Vfyy)  ♦  |y(U  V^y>  -  (Uf)yVxy  . 

^(»^xy)  •(«^xV'“,^y*('V»  • 

ax'  x'  Ty1  y'  y  xy  y  xy  x  xy 

C  A  ]  C  B  ] 

V  f  f  ♦  u  ^  f  ♦  u  +vV'5'  +  V  f  a 

x  1  xy  x  yy  *x  yy  y  'xx  y  'xx 

C  -A  ]  <-vtyfyy)  C  -B  ]  (-u*x*xx) 

|_(V6  .  Vv)  *  fy<-«  .  »«)  -  a  -  v  - 

V  i_(i^2)  .  u  —  (i'f2)  = 

jy7  y'  ax7Tx' 
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Appendix  A. 2 


It  is  to  be  shown  that: 

J*\+l/2’hk+l/2^  =  J^k»hk+l/2^  =  J^k+l,hk+l/2^ 


k,nk+l/2  1  =  Ul  1 k+1 ,nk+l/2 
fc 

k  T  "k+1  k+1/2  ’  gk+l/21  'f|c+1  Tk‘ 


u,  Hk+1  ^k  +  Hk^k+1  w  12 _ /Tt,  ^  * 

where  'k+i/2  ~  - H — V  H  i -  *  "*  =  - 1  *  “u+i  -  t  i,  /  • 


Thus,  using  the  well-known  properties  of  the  Jacobian: 

>  ,.Hk+l  \  +  Hk  \+l  fo  ,4,  _  t>  \\  _ 

n  ji  Vi/2*hk+i/2}  = J{ — nr^" - *  ^ J  " 

f. 


Hk  +  ^k+1  *  9k+l/2  k+1  k 


9k+i/2°V  HknT  CJ(Hk+1^k’  ^k+l)  +  J(Hk^kn*-^k)]  " 

9k+l/2(V  ^CHk+lJ(tk’tk+l)  +  HkJ(^k’^k+l)] 
r-£“J(tk*Vi) 

9k+l/2  K  K  1 

H)  J(tk,hk+1/2)  =  J(%.STTTT^k+l-  *k»  - 


III)  J(  ^k+l'hk+l/2*  =  J(  ^k+1  ’gk+1/2(  ^k+1"  ^k^ 


J(^kn*^k)  -  r^Jf^k'Vi1' 

gk+l/2  k  1  k  9k+l/2  K 


as 


References 


Bryden,  H.L.  (1976)  Horizontal  advection  of  temperature  for  low-frequency 
motions.  Deep-Sea  Res.,  23,  1165-1174. 

Bryden,  H.L.  (1979)  Poleward  heat  flux  and  conversion  of  available 
potential  energy  in  Drake  Passage.  J.  Mar.  Res.,  37U),  1-22. 

Bryden,  H.L.  (1980)  Geostrophic  vorticlty  balance  in  midocean. 

J.  Geophys.  Res.,  85  (C5),  2825-2828. 

Bryden,  H.L.  (1982)  Sources  of  eddy  energy  in  the  Gulf  Stream 
recirculation  region.  J.  Mar.  Res.,  40  (4),  1047-1068. 

Chow,  J.C.  and  W.R.  Holland  (1985)  A  multi-layer  quasi -geostrophic 
numerical  model  of  ocean  circulation.  In  preparation. 

Dewar,  W.K.,  and  J.M.  Bane  (1984)  The  subsurface  energetics  of  the  Gulf 
Stream  near  the  Charleston  Bump.  J.  Phys.  Oceanog.,  accepted. 

Fofonoff,  N.P.  (1954)  Steady  flow  In  a  frictionless  homogeneous  ocean. 

J.  Mar.  Res..  13,  254-262. 

Fofonoff,  N.P.  (1962)  Dynamics  of  ocean  currents.  The  Sea;  Ideas  and 
Observations  on  Progress  In  the  Study  of  the  Seas,  Vol.  1.  Physical 
Oceanography .  M.N.  Hill,  Ed.,  Wlley-Interscience,  323-395. 

Fofonoff,  N.P.  (1981)  The  Gulf  Stream  System.  In:  Evolution  of  Physical 
Oceanography,  B.A.  Warren  and  C.  Wunsch,  eds..  The  HIT  Press, 
Cambridge,  Massachusetts,  112-139. 

Fofonoff,  N.P.  and  M.M.  Hall  (1983)  Estimates  of  mass,  momentum,  and 
kinetic  energy  fluxes  of  the  Gulf  Stream.  J.  Phys.  Oceanog.,  13 
(10),  1868-1877. 


161 


Fugllster,  F.C.  (1963)  Gulf  Stream  '60.  Progress  In  Oceanography,  Vol . 

1,  Pergamon  Press,  265-383. 

Haldvogel,  D.B.,  and  W.R.  Holland  (1978)  The  stability  of  ocean  currents 
In  eddy-resolving  general  circulation  models.  J.  Phys.  Oceanog.,  8^ 
393-413. 

Halkln,  D.  and  T.  Rossby  (1985)  The  structure  and  transport  of  the  Gulf 
Stream  at  73*  W.  J.  Phys.  Oceanog.,  submitted. 

Hall,  M.M.  and  H.L.  Bryden  (1985)  Profiling  the  Gulf  Stream  with  a 
current  meter  mooring.  Geophys.  Res.  Letters,  accepted. 

Halllwell,  G.R.,  Jr.,  and  C.N.K.  Mooers  (1983)  Meanders  of  the  Gulf 

Stream  downstream  from  Cape  Hatteras  1975-1978.  J.  Phys.  Oceanog., 
13,  1275-1292. 

Harrison,  D.E.  (1979)  Eddies  and  the  general  circulation  of  numerical 
model  gyres:  an  energetic  perspective.  Rev.  Geophys.  and  Space 
Phys.,  17  (5),  969-979. 

Harrison,  D.E.  and  A.R.  Robinson  (1978)  Energy  analysis  of  open  regions 
of  turbulent  flows  --  mean  eddy  energetics  of  a  numerical  ocean 
circulation  experiment.  Dyn.  Atmos.  Oceans,  2,  185-211. 

Holland,  W.R.  (1978)  The  role  of  mesoscale  eddies  In  the  general 
circulation  of  the  ocean  —  numerical  experiments  using  a 
wind-driven  quasi -geos trophic  model.  J.  Phys.  Oceanog.,  8,  363-392 

Holland,  W.R.  (1985)  Simulation  of  mesoscale  ocean  variability  In 
midlatitude  gyres.  In:  Advances  In  Geophysics,  Vol.  xx. 

Holland,  W.R.  and  D.B.  Haldvogel  (1980)  A  parameter  study  of  the  mixed 
Instability  of  Idealized  ocean  currents.  Dyn.  Atmos.  Oceans,  £, 
185-215. 


162 


Holland,  VI. R.  and  L.B.  Lin  (1975)  On  the  generation  of  mesoscale  eddies 
and  their  contribution  to  the  oceanic  general  circulation.  I.  A 
preliminary  numerical  experiment.  J.  Phys.  Oceanog.,  5^,  642-657. 

Holland,  W.R.,  and  P.B.  Rhlnes  (1980)  An  example  of  eddy-induced  ocean 
circulation.  J.  Phys.  Oceanog.,  10  (7),  1010-1031. 

Holland,  W.R.,  and  W.J.  Schmitz  (1985)  The  penetration  scale  of  the  Gulf 
Stream.  J.  Phys.  Oceanog.,  submitted. 

Johns,  W.E.  (1985)  Dynamics  and  structure  of  Gulf  Stream  meanders 

northeast  of  Cape  Hatteras,  N.C.  Ph.  D.  Thesis,  University  of  Rhode 
Island,  232  pp. 

Johns,  W.E.  and  D.R.  Watts  (1985)  Gulf  Stream  meanders:  Observations  on 
the  deep  currents.  J.  Geophys.  Res.,  accepted. 

Knauss,  J.A.  (1969)  A  note  on  the  transport  of  the  Gulf  Stream. 

Frederick  C.  Fugllster  Sixtieth  Anniversary  Volume,  Deep-Sea  Res. , 

16  (Suppl.), 117-123. 

Luyten,  J.R.  (1977)  Scales  of  motion  In  the  deep  Gulf  Stream  and  across 
the  Continental  Rise.  J.  Mar.  Res. ,  35  (1),  49-74. 

Luyten,  J.R. ,  J.  Pedlosky  and  H.  Stommel  (1983)  The  ventilated 
thermocllne.  J.  Phys.  Oceanog.,  13  (2),  292-309. 

Luyten,  J.R.,  and  A.R.  Robinson  (1974)  Transient  Gulf  Stream  meandering. 
Part  II:  Analysis  via  a  quasi -geostrophlc  time-dependent  model. 

J.  Phys.  Oceanog.,  4  (20),  256-269. 

McCartney,  M.S.  (1982)  The  subtropical  recirculation  of  Mode  Waters. 

J.  Mar.  Res.,  40,  Suppl.,  427-464. 


V'.v 


.v.'-.v; 


•  • 

1777777 


\V.-4v:v 
■w  aV-V« 


.N  . 


163 


McDowell,  S.,  P.  Rhines,  and  T.  Keffer  (1982)  North  Atlantic  potential 
vortlclty  and  Its  relation  to  the  general  circulation.  J.  Phys. 
Oceanog.,  12  (12),  1417-1436. 

Munk,  W.H.  (1950)  On  the  wind-driven  ocean  circulation.  Journal  of 
Meteorol ogy ,  7_,  79-93. 

Orlanskl,  I.  (1969)  The  Influence  of  bottom  topography  on  the  stability 
of  jets  In  a  baroclinlc  fluid.  J.  Atmos.  Scl.,  26,  1216-1232. 

Owens,  W.B.  (1984)  A  synoptic  and  statistical  description  of  the  Gulf 
Stream  and  subtropical  gyre  using  SOFAR  floats.  J.  Phys.  Oceanog., 
14,  (1),  104-113. 

Pedlosky,  J.  (1970)  Finite-amplitude  baroclinlc  waves.  J.  Atmos.  Sci., 

27  (1),  15-30. 

Pedlosky,  J.  (1979)  Geophysical  Fluid  Dynamics.  Sprlnger-Verlag,  624  pp. 


Pedlosky,  J.  (1982)  Finite-amplitude  baroclinlc  waves  at  minimum  critical 
shear.  J.  Atmos.  Scl.,  39  (3),  555-562. 

Pratt,  L.J.  and  M.E.  Stern  (1985)  Dynamics  of  vortlclty  fronts. 

J.  Fluid  Hech.,  submitted. 

Raymer,  M.E.,  A.  Spencer,  and  H.L.  Bryden  (1985)  Temperature  correction 
for  mooring  motion  on  GUSTO  current  meter  time  series.  W.H.O.I. 
Technical  Report,  85. 

Rhines,  P.B.  and  W.R.  Young  (1982)  Homogenization  of  potential  vortlclty 
In  planetary  gyres.  J.  Fluid  Mech. ,  122,  347-368. 

Richardson,  P.L.  (1983a)  A  vertical  section  of  eddy  kinetic  energy 
through  the  Gulf  Stream  System.  J.  Geophys.  Res.,  88  (C4), 
2705-2709. 

Richardson,  P.L.  (1983b)  Eddy  kinetic  energy  In  the  North  Atlantic  from 
surface  drifters.  J.  Geophys.  Res.,  88  (C7),  4355-4367. 


.1 


-I 

m 


Robinson,  A.R.,  J.R.  Luyten,  and  G.  Flier  1  (1975)  On  the  theory  of  thin 
rotating  jets:  A  quasl-geostrophlc  time  dependent  model. 

Geophysical  Fluid  Dynamics,  Vol.  6. 

Robinson,  A.R.  and  P.P.  Nlller  (1967)  The  theory  of  free  Inertial 
currents.  I.  Path  and  structure.  Tellus,  19,  269-291. 

Schmitz,  W.J.,  Jr.  (1977)  On  the  deep  general  circulation  In  the  Western 
North  Atlantic.  J .  Mar .  Res ■ ,  3!5  (1),  ?l-28. 

Schmitz,  W.J.,  Jr.  (1980)  Weakly  depth- dependent  segments  of  the  North 
Atlantic  circulation.  J.  Mar.  Res. ,  38  (1),  111-133. 

Schmitz,  W.J.,  Jr.  (1984)  Abyssal  eddy  kinetic  energy  in  the  North 
Atlantic.  J.  Mar.  Res.,  42,  509-536. 

Stommel,  H.  (1948)  The  westward  Intensification  of  wind-driven  ocean 
currents.  Trans.  Amer.  Geophys.  Union,  29^  202-206. 

Stommel,  H.  (1950)  The  Gulf  Stream:  A  brief  history  of  the  Ideas 
concerning  Its  cause.  The  Scientific  Monthly,  70  (4),  242-253. 

Stommel,  H.  (1965)  The  Gulf  Stream.  University  of  California  Press, 
erkeley,  248  pp. 

Talley,  L.D.  (1982)  Instabilities  and  radiation  of  thin,  barocllnlc 
jets.  Ph.  0.  Thesis,  Hass.  Inst,  of  Tech. /Woods  Hole  Oceanog. 

Inst.,  233  pp. 

Warren,  B.A.  (1963)  Topographic  Influences  on  the  path  of  the  Gulf 
Stream.  Tellus,  15 ,  167-183. 

Warren,  B.A.  and  G.H.  Volkmann  (1968)  Measurement  of  volume  transport  of 
the  Gulf  Stream  south  of  New  England.  J .  Mar .  Res . ,  26  (2),  110-126 


Watts,  D.R.  (1983)  Gulf  Stream  Variability.  In:  Eddies  In  Marine 
Science,  A.R.  Robinson,  ed.,  Sprlnger-Verlag,  114-144. 

Watts,  D.R.  and  W.E.  Johns  (1982)  Gulf  Stream  meanders:  observations  on 
propagation  and  growth.  J.  Geophys.  Res.,  87  (C12),  9467-9476. 
Worthington,  L.V.  (1976)  On  the  North  Atlantic  Circulation.  The  Johns 


Hopkins  Oceanographic  Studies  6  110  pp. 

it,  D.G.  (1981)  Barocllnlc  Instability  in  Drake  Passage.  J.  Phys 


i 


DOCUMENT  LIBRARY 
April  9,  1985 

DISTRIBUTION  LIST  FOR  TECHNICAL  REPORT  EXCHANGE 


Institute  of  Marine  Sciences  Library 
University  of  Alaska 
O'Neill  Building 
905  Kovukuk  Ave..  North 
Fairbanks.  AK 

Attn:  Stella  Sanchez-Wade 
Documents  Section 
Scripps  Institution  of  Oceanography 
Librarv.  Mail  Code  C-075C 
La  Jolla,  CA  92093 

Hancock  Library  of  Biology  &  Oceanography 

Alan  Hancock  Laboratory 

University  of  Southern  California 

University  Park 

Los  Angeles.  CA  90089-0371 

Gifts  &  Exchanges 
Library 

Bedford  Institute  of  Oceanography 
P.O.  Box  1006 

Dartmouth.  NS,  B2Y  4A2,  CANADA 

Office  of  the  International 
Ice  Patrol 

c  o  Coast  Guard  R  &  D  Center 
Avery  Point 
Groton,  CT  06340 

Library 

Physical  Oceanographic  Laboratory 
Nova  University 
8000  N.  Ocean  Drive 
Dania.  FL  33304 

NOAA.  EDIS  Miami  Library  Center 
4301  Rickenbacker  Causeway 
Miami.  FL  33149 

Library 

Skidaway  Institute  of  Oceanography 
P.O.  Box  13687 
Savannah,  GA  31416 

Institute  of  Geophysics 
University  of  Hawaii 
Library  Room  252 
2525  Correa  Road 
Honolulu,  HI  96822 


MIT  Libraries 

Serial  Journal  Room  14E-210 
Cambridge,  MA  02139 

Director,  Ralph  M.  Parsons  Laboratory 

Room  48-311 

MIT 

Cambridge,  MA  02139 

Marine  Resources  Information  Center 

Bldg.  E38-320 

MIT 

Cambridge,  MA  02139 
Library 

Lamont-Doherty  Geological  Observatory 
Colombia  University 
Palisades,  NY  10964 

Library 

Serials  Department 
Oregon  State  University 
Corvallis,  OR  97331 

Pell  Marine  Science  Library 
University  of  Rhode  Island 
Narragansett  Bay  Campus 
Narragansett,  RI  02882 

Working  Collection 
Texas  A&M  University 
Dept,  of  Oceanography 
College  Station,  TX  77843 

Library 

Virginia  Institute  of  Marine  Science 
Gloucester  Point,  VA  23062 

Fisheries-Oceanography  Library 
151  Oceanography  Teaching  Bldg. 
University  of  Washington 
Seattle.  WA  98195 

Library 

R.S.M.A.S. 

University  of  Miami 

4600  Rickenbacker  Causeway 

Miami,  FL  33149 


Library 

Chesapeake  Bay  Institute 
4800  Atwell  Road 
Shady  Side.  MD  20876 


I.  RfPOPT  NO. 


WHOI-85-16 


REPORT  DOCUMENTATION 
PAGE 


4.  Titl.  and  SuMftla 

Horizontal  and  Vertical  Structure  of  Velocity,  Potential  Vorticity  and 
Energy  in  the  Gulf  Stream 


7.  Authors) 

Melinda  M.  Hall _ 


9.  Pirformlni  Organization  Name  and  Address 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Massachusetts  02543 


S.  Recipient's  Accession  i 


12.  Sponsoring  Organization  Name  and  Addrwt 

National  Science  Foundation 
and 

Office  of  Naval  Research 


ft.  Itapart  Dal* 

February  1985 


I.  Parforming  Organization  Rapt.  No. 

WHOI-85-16 


10.  Praiact/TaalimBrti  Unit  No. 


II.  Conlroct(C)  or  Grant(G)  No. 

(ciOCE-8208746 
N00014-82-C-0019, 
NR  083-004 


19.  Type  of  Rapoit  4  Period  Covered 

Ph.D. 


IS.  Supplementary  Not  as 

This  thesis  should  be  cited  as:  Melinda  M.  Hall,  1985,  Horizontal  and  Vertical  Structure  of  Velocity, 
Potential  Vorticity  and  Energy  in  the  Gulf  Stream.  Ph.D.  Thesis.  MIT/WHOI,  WHOI-85-16. 


IS.  AMtract  (Limit:  200  wonts)  ->«>-- 

A  current  meter  mooring  was  deployed  for  one  year  at  the^  mean  position  of  the  Gulf  Stream 
at  68^ W.  Projecting  the  velocity  onto  along-  and  cross-stream  components  shows  that  the  vertical 
structure  of  velocity  at  the  site  is  described  by  a  baroclinic  along-stream  and  decoupled,  barotropic 
cross-stream  component.  The  well-defined  relationship  between  thermocline  temperature  and 
cross-stream  distance  allows  the  use  of  the  former  as  a  horizontal  coordinate  and  enables  estimates 
of  mass,  momentum  and  kinetic  energy  fluxes  of  the  Gulf  Stream  for  four  events  when  it  passes 
over  the  mooring  site.  Alternatively,  an  average  vertical  and  horizontal  velocity  section  of  the 
Stream  may  be  constructed,  and  using  that,  a  potential  vorticity  section.  There  is  a  strong  jump 
from  low  to  high  values  of  thermocline  potential  vorticity  crossing  the  Stream  from  south  to  north. 
Time-series  of  vertical  velocities  have  been  calculated  and  their  rms  sizes  are  large  (0.05  cm/s). 
The  Rossby  number  below  the  thermocline  is  estimated  to  be  $  0.1,  but  in  the  thermocline  it  is  about 
0.3.  Diagnostic  calculations  of  the  kinetic  energy  budgets  in  a  quasi-geostrophic  eddy-resolving 
general  circulation  model  suggest  that  downstream  convervgence  of  mean  kinetic  energy  in  the 
Gulf  Stream  is  primarily  balanced  by  mean  ageostrophic  flow  against  a  pressure  gradient. 


/ 


17.  Document  Analyst*  a.  Descriptors 


1.  Gulf  Stream  structure 

2.  Energy  budgets  of  free  jets 

3.  potential  vorticity  in  Gulf  Stream 


b.  Idantlfiare/Opan-Endad  Tarms 


e.  COSAT!  Flald/Oroup 


IS.  Availability  Stataman: 

I  Approved  for  publication;  distribution  unlimited. 


IS.  Sacurtly  Claaa  (Tbl*  A* port) 

UNCLASSIFIED 

SO.  Sacurtly  Claaa  (TMa  !’•«*) 


(Saa  ANSI-23*.  IS) 


Saa  Instruct lani  an  I 


Zl.  Ho.  at  Paco* 

165 _ 

72.  Prica 

OPTIONAL  POma  Z7Z  (4-77) 
(Formally  NTIS-3S) 
Oapartmant  at  Cemmatca 


END 

FILMED 

7-85 

DTIC 


